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V.  CW  LASER  ROCKET  ENGINE 


1.  INTRODUCTION  AND  ENGINE  CONFIGURATIONS 

Two  basic  engine  configurations  have  evolved  for  CW  laser  propulsion 
The  character  of  the  designs  depends  significantly  upon  the  manner  in  which 
the  laser  radiation  is  introduced  into  the  absorption  chamber.  In  one  case 
the  laser  beam  enters  the  chamber  passing  through  the  rocket  exhaust  plume 
without  having  its  energy  absorbed  in  the  plume  or  nozzle  (see  Section  V.  1). 
This  concept  comprises  the  single  port  engine.  In  the  other  configuration 
the  beam  enters  the  chamber  through  a window  (solid  or  aerodynamic ) with 
all  of  the  absorption  occurring  in  the  chamber.  The  nozzle  and  plume  part 
of  this  so-called  two-port  engine  rocket  are  not  different  from  conventional 
chemical  rocketry. 

a.  Single-Port  Engine;  Bleaching 

In  Figure  33  we  see  the  schematic  of  a single-port  CW  laser  propul- 
sion engine.  The  laser  beam  enters  through  the  nozzle  which  also  serves 
as  a collector  for  the  beam.  The  propellant  used  with  this  scheme  is  fed 
from  a tank  into  the  absorption  chamber  where  it  absorbs  laser  radiation 
via  molecular  or  inverse  bremstrahlung  absorption.  Some  quant itative 
estimates  of  laser  absorption  waves  and  the  absorption  zone  structure  are 
given  in  Section  V.  3. 

A propellant  that  does  not  absorb  in  the  plume  or  nozzle  but  does 
absorb  in  the  chamber  to  raise  the  stagnation  enthalpy  is  required.  Further- 
more, a nearly  step  function  behavior  of  the  absorption  coefficient  between 

1 1 
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the  chamber  and  the  nozzle  is  preferred  In  general,  the  temperature  will 
not  have  a step  function  variation  with  engine  location  and  bleaching  will  not 
occur.  The  temperature  variation  will  depend  upon  the  nozzle  geometry 
and  on  the  kinetic  processes  (dissociation  and  recombination).  A slow 
temperature  variation,  for  example,  would  keep  the  absorption  coefficient 
high  in  the  nozzle.  This  would  only  permit  a fraction  of  the  radiation  to 
reach  the  chamber  and  hence  alter  the  entire  absorption  process  upon  which 
the  engine  was  originally  designed.  Critical  experiments  for  the  single- 
port configuration  are  proposed  in  Section  V.  5. 
b.  Two-Port  Engine 

In  Figure  34  there  is  a schematic  of  a two-port  cw  laser  propulsion 
engine.  The  laser  beam  enters  through  a side  port  onto  a focusing  mirror 
which  deflects  the  beam  through  an  aerodynamic  window.  The  propellant 
and  the  aerodynamic  window  source  flow  medium  are  the  same.  Here,  as 
in  the  single-port  device,  selected  seedants  can  be  ado.'d  to  enhance  or 
suitably  modify  the  absorptive  characteristics  of  the  flow.  Note  again  the 
need  for  an  absorption  coefficient  that  absorbs  the  radiation  in  the  chamber 
under  high  pressure  conditions.  In  the  next  section  the  thermodynamics  of 
various  candidate  propellants  is  discussed. 

The  two  key  areas  of  research  needed  for  the  single-port  engine  are 
also  needed  for  the  tw'o-port  scheme.  In  addition  there  is  a third  important 
area  of  concern  - the  aerodynamic  window  flow.  A supersonic  aerodynamic 
window  flow  can  be  designed  to  operate  at  a practically  infinite  pressure 
ratio  between  the  high  pressure  chamber  and  the  low  pressure  ambient. 

The  precise  flow  characteristics  in  the  vicinity  of  the  aerodynamic  window 
are  rather  complex;  further  downstream  approaching  the  focus  of  the  laser 
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beam  and  the  absorption  chamber,  there  is  a supersonic  flow  which  must 
be  decelerated  to  subsonic  velocities  via  a normal  shock.  The  focus  of  the 
beam  (absorption  zone)  is  in  the  subsonic  portion  of  the  chamber  flow.  The 
energetics  of  the  essentially  constant  pressure  heating  is  then  completely 
analogous  to  the  processes  occurring  in  chemical  propulsion  systems. 

Since  the  absorption  processes  will  induce  laser  absorption  waves  in  the 
chamber  that  will  travel  toward  the  beam,  the  aerodynamic-window  flow 
configuration  will  be  perturbed. 

It  is  obvious  that  problems  remain  to  be  solved  for  the  two-port 
configuration.  We  have  examined  some  of  the  matters  described  above. 

In  Section  V.4  we  describe  the  gasdynamic  matching  analysis  for  the  two- 
port  system  with  the  purpose  of  obtaining  scaling  maps  that  size  experi- 
mental engines  for  various  power  levels.  In  Section  V.  5 we  describe  a 
burner  experiment  that  can  potentially  address  many  of  the  questions  we 
have  posed  in  the  present  discussions.  We  must  add  that  most  of  our 
attention  has  been  paid  to  the  two-port  engine.  This  reflects  our  view  that 
the  increased  flexibility  of  this  configuration,  as  contrasted  to  the  inflexible 
need  for  absorption  coefficient  tailoring  with  the  single-port  engine,  has 
more  promise  for  early  success.  We  must  add  that  whatever  is  learned 
regarding  the  physics  of  the  absorption  process  will  be  applicable  to  both 
CW  laser  propulsion  engine  configurations. 
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2.  THERMODYNAMICS 


The  thermodynamics  of  addition  of  energy  to  the  flowing  medium  in 
the  absorption  chamber  will  be  discussed.  There  are  a number  of  connected 
areas  of  interest.  Firstly,  we  define  the  enthalpy  increase  required  to  pro- 
vide a certain  specific  impulse  (for  a specified  mission).  A laser  system 
must  be  capable  of  delivering  the  enthalpy  increase  to  the  propellant  of 
interest.  Secondly,  several  propellants  will  be  described  and  assessments 
regarding  their  applicability  will  be  made.  High  temperature  water  appears 
to  be  a particularly  interesting  propellant.  Thirdly,  certain  practical 
aspects  of  the  thermodynamic  efficiency  of  propellants  as  well  as  loss  mech- 
anisms will  be  discussed, 
a . Specific  Impulse 

An  important  measure  of  a propulsion  system  is  its  specific  impulse. 
Under  uniform  exhaust  flow  conditions,  i.  e.  , all  the  fluid  particles  are 
leaving  the  rocket  nozzle  exit  plane  at  the  same  velocity,  the  specific  im- 
pulse is  defined''  as 

Ig  = U^/g  (seconds)  (277) 

where  U-,  is  the  constant  exit  speed  and  g is  the  acceleration  of  gravity. 
Since  the  chamber  speed  is  typically  a fraction  (~l/lO)  of  the  sound  speed, 
the  kinetic  energy  per  unit  mass  in  the  chamber  will  be  small  compared  to 
the  thermal  energy  introduced  via  laser  radiation.  Hence  the  energy  equa- 
tion between  the  absorption  chamber  and  the  nozzle  exit  is 

U?  = n/“2  (h5"-”h7)  (278) 

The  subscript  numbering  system  refers  to  the  stations  in  a configured 
two-port  laser  propulsion  engine  (see  Figure  52  and  Table  4). 


where  h,.  is  the  chamber  enthalpy  and  hy  is  the  exhaust  nozzle  enthalpy 
The  expansion  of  the  flow  in  the  nozzle  will  reduce  the  thermal  energy 
continually.  Hence,  to  maximize  the  exit  velocity  is  set  to  zero.  Then 
the  maximum  specific  impulse  I is  related  to  the  enthalpy  change  by: 


Ah5  = 


(279) 


Note  that  h^  is  replaced  by  Ah^  since  enthalpy  differences  from  some 
common  reference  state  are  always  considered.  Relation  (279)  is  plotted 
in  Figure  35.  Chemical  rockets  have  propulsion  systems  that  are  theo- 
retically capable  of  about  500  seconds  as  shown  in  Table  3.  The  chemical 
regime  is  shaded  on  the  figure.  Laser  propulsion  rockets  are  not  confined 
to  the  chemical  energy  of  reaction  between  fuel  and  oxidizer  as  a limiting 
factor  for  Ahg.  Conceivably,  it  should  be  possible  to  attain  specific  im- 
pulses in  the  thousands  of  seconds.  Mission  analysis  studies  (see  Section  II) 
have  shown,  however,  that  specific  impulses  between  500  and  1000  seconds 
are  of  considerable  interest.  Laser  propulsion  is  needed  because  there  are 
no  existing  systems  that  operate  in  this  range. 

Based  upon  Figure  35  it  is  evident  that  one  requires  an  enthalpy 
increase  of  at  least  104  J/gm  to  obtain  a specific  impulse  of  ~500  seconds. 
This  statement  reflects  a perfectly  efficient  process.  If  one  factors  in  the 
potential  losses,  then  one  must  supply  something  like  30  or  40%  excess 
Ah^  to  achieve  the  specific  impulses  indicated.  Note  that  LH-, -LO,  c hemical 
rockets  are  ~60%  efficient.  Loss  terms  will  be  itemized  in  a subsequent 
subsection;  they  are  dominated  by  the  heat  transfer  to  the  walls  in  the  ab- 
sorption chamber, 
b.  Propellants 

The  specific  impulse  variation  with  Ah^  shown  in  Figure  35  defines 
a criterion  for  the  selection  of  a propellant.  We  must  provide  at  least  104 
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j/gm  to  a laser  propulsion  rocket  absorption  chamber  to  exceed  specific- 
impulses  achievable  by  chemical  propulsive  systems. 

This  statement  is  significant  from  a number  of  points  of  view. 

(i)  Any  propellant  conceivable  will  have  as  much  enthalpy  as  desired 
if  one  is  willing  to  operate  over  large  temperature  and  pressure  ranges. 

(ii)  Large  (temperature  and/or  pressure)  excursions  will  not  be  a 
favorable  circumstance  from  the  point  of  view  of  the  absorption  process 
since  the  absorption  coefficients  (for  molecular  or  Inverse  Brenistrahlung) 
depend  critically  upon  the  temperature  and  pressure. 

(iii)  The  extraction  (or  efficient  use)  of  the  available  enthalpy 

is  an  additional  important  factor  influencing  propellant  choice.  Practical 

rocket  nozzle  expansion  bells  have  area  ratios  that  are  perhaps  between  10 

and  100;  typically  the  ratio  is  near  40.  Improvement  in  the  state-of-the-art 

3 

may  perhaps  permit  an  area  ratio  of  up  to  1 0 or  even  higher,  but  there  will 
eventually  be  a practical  limitation. 

(iv)  Enthalpy  extraction  from  the  heated  gas  may  be  restricted  due 

to  dissociation  and  vibrational  specific  heat  lag.  Thermodynamic  properties 
of  propellants  are  most  usually  displayed  in  Mollier  charts,  in  which  en- 
thalpy is  plotted  against  entropy, 
c.  H-.Q  Mollier  Chart 

In  Figure  36  a Mollier  Chart  for  H^O  is  illustrated.  Enthalpy  and 
entropy  are  the  ordinate  and  abscissa  on  the  Mollier  Chart;  the  thermo- 
dynamic variables  illustrated  in  Figure  36  are  the  pressure,  p in  atmo- 
spheres, the  temperature,  T in  °K  and  the  compressibility  factor  Z 
which  is  defined  by  the  gas  law  p = Zp  Rq  T where  p is  the  density  and 

R is  the  universal  gas  constant.  Z is  a measure  of  the  dissociation  of 
o * 
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the  gas.  The  chart  was  generated  utilizing  the  AERL  equilibrium  chemistry 
program.  Note  immediately  that  there  is  sufficient  available  enthalpy  to 
achieve  an  excess  of  500  sec  specific  impulse.  Obviously,  there  is  sufficient 
enthalpy  available  for  all  candidate  propellants,  however  for  11^,0  it  is 
available  at  reasonably  low  temperatures.  In  the  next  section  one  observes 
that  a 60%  efficient  process  requires  temperatures  around  5000°K  for 
interesting  pressures. 

Particularly  interesting  boundaries  to  put  on  the  Mollier  chart  for 
H-,0  are  the  absorption  coefficient  contours  for  1%  molar  Al^O^  in  H^O, 
since  Al^O^  dissociates  to  AlO  which  is  a very  strong  absorber  of  10.  6 pm 
radiation.  These  boundaries  are  shown  on  the  skeleton  of  the  Mollier  chart 
in  Figure  37.  Note,  in  particular,  the  nesting  of  large  values  of  the  ab- 
sorption coefficient  in  the  high  pressure  region  of  the  chart,  which  is  for- 
tunate since  initial  high  pressure  facilitates  enthalpy  extraction, 
d.  Efficiency 

We  now  evaluate  the  efficiency  of  the  chamber  absorption  (heating! 
process  for  candidate  propellants.  The  efficiency  e is  defined  as 


€ = 


Ah&  - Ah 

aTT^ 


7 


(2801 


where  Ah-,  is  the  nozzle  exit  enthalpy.  If  the  exit  enthalpy  is  zero,  this 
means  that  the  exhaust  products  in  the  nozzle  are  at  the  same  state  that 
they  were  prior  to  the  absorption  of  laser  radiation  - hence  the  process 
is  perfectly  efficient.  Let  vis  now  relate  the  efficiency  to  the  specific  im- 
pulse, by  vising  Eq.  (278).  Thvis, 


2 . 2 

g I 

££_ 

2 Al-lr 


= 1 - 


Ah. 

Ahc 


(281) 


22 


ENTHALPY  ( joule*  /gram) 


OOl  OOOl 


T • 57 00*K 
5500 
5000 
4500 


4000 

5500 

3000 


A (0  ABSORPTION  COUNTOURS 


2500 


ENTHALPY  REFERENCED 
TO  ZERO  AT  298°K 


2000 


- — 1500 
1000 

T • 600* K 


00  700  8 00 

ENTROPY  (cal  /gram  °K) 


h- 


From  Eq.  (281)  wo  observe  that  specifying  I and  t determines  Ah, 
and  Ah-,.  Then  starting  from  a base  level,  we  can  obtain  values  of  the 
exit  pressure  p„,  for  given  values  of  the  absorption  chamber  pressure  pr  . 
The  pressure  ratio  can  further  be  related  to  the  area  ratio  of  the  nozzle 
The  variation  of  'Pj.  and  A^/A^  (A ^ is  the  nozzle  throat  area)  with  Mach 

number  as  the  parameter  for  constant  \ flows  are  given  below: 


On  Figure  38  we  show  the  variations  cf  p_  p_  and  A„/A,  for  various 

7 5 i o 

values  of  y . Even  though  the  nozzle  expansion  process  is  one  of  variable 
\ , one  can  still  obtain  valuable  information  from  a constant  \ calculation. 
Two  values  of  the  area  ratio  are  of  interest  as  they  relate  to  the  pressure 
ratio.  An  area  ratio  of  10  represents  a practical  (nonphysical)  upper 
bound  to  that  ratio.  Rocket  nozzles  just  become  too  long.  (Note  the  wide 
variation  of  p^/p.,  depending  upon  y . ) Another  interesting  value  is  for 
an  area  ratio  of  40  which  represents  a typical  high  performance  rocket 
expansion  bell  area  ratio. 

e.  Efficiency  of  Water  Rocket 

We  apply  the  preceding  ideas  to  the  water  rocket,  i.  e.  , the  11,0 
Mollier  chart  described  in  the  previous  section.  A specific  impulse  of 
700  sec  is  selected.  This  alters  the  Ahj.  depending  upon  the  efficiency; 
however,  the  difference  between  Ahj.  and  Ah-,  is  constant,  to  provide 
the  700  sec  of  specific  impulse.  The  cases  of  30%,  60%  and  °0%  efficiency 
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Figure  38  Kxpansion  Pressure  Ratio  Related  to  Area  Ratio 
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are  shown  on  the  skeletal  H20  Mollier  chart  in  Figure  39.  The  cycle  base- 
line has  been  taken  as  1000°K.  Note  that  the  90%  efficiency  case  utilizes 
nearly  all  of  its  available  enthalpy.  Also,  the  thermodynamic  process 
between  chamber  and  nozzle  exit  has  been  taken  to  be  an  isentropic  one. 

This  is  usually  a good  approximation.  For  each  efficiency  we  have  taken 
four  chamber  pressures  that  drop  to  four  different  exit  pressures  through 
the  isentropic  expansion  process.  These  cases  translate  to  Figure  40 
which  is  an  efficiency  map  that  has  been  generated  for  a constant  y = 1.3. 
This  value  of  y is  considered  to  be  an  average  value  for  H^O  as  it  trav- 
erses between  the  chamber  and  the  exit  nozzle.  We  note  that  at  a chamber 
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pressure  of  about  100  atm  and  a nozzle  bell  area  ratio  between  10  and  10 
the  efficiency  is  55  to  75%.  This  is  comparable  to  efficiencies  achieved 
in  LH,  - LO-,  rocket  engines.  These  engines,  however,  cannot  obtain  a 
specific  impulse  of  700  sec  as  we  would  be  obtaining  with  a laser  propulsion 
scheme.  It  is  to  be  noted  that  even  though  decreasing  the  chamber  pressure 
decreases  the  efficiency,  one  may  achieve  interesting  specific  impulse  levels 
at  low  temperatures.  Specific  impulse  as  function  of  temperature,  derived 
from  the  Mollier  chart,  is  shown  in  Figure  41. 
f.  Losses 

The  propellant  enthalpy  is  reduced  by  heat  transfer  losses  to  the 
confining  walls  of  the  engine.  A good  discussion  of  heat  transfer  in  the 

(4) 

rocket  chamber  environment  is  given  by  Summerfield. 


(4)  Summerfield,  M.  , "The  Liquid  Propellant  Rocket  Engine  in 

Princeton  Aeronautical  Paperbacks,"  Number  1,  Liquid  Propellant 
Rockets,  Princeton  University  Press  (I960). 
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He  finds  that  the  convective  heat  transfer  is  given  by  the  expression 


% 


Nu  k 
x 

L 


(T  - T ) 
' x w 


(284) 


where  sub  x refers  to  the  station  of  interest,  T is  the  wall  temperature 


and  the  Nusselt  Number  is  given  by 


Nux  = a <Rex)4/5  (Pr)l/3 


(285) 


The  quantity  a is  a constant,  the  Reynolds  Number  is  based  upon  local 
conditions  and  the  relevant  length  L^.  Pr  is  the  Prandtl  Number  and  k 
is  the  thermal  conductivity  of  the  propellant.  Consequently,  the  convective 
heat  transfer  is 

q = K (L  f1/5  (T  )4/5  (P  )4//S  (286) 

MC  X X X 

with  K a constant.  The  transport  properties  (viscosity  and  thermal  con- 
ductivity) were  approximated  to  be  linearly  increasing  functions  of  the 
temperature  T . Relation  (286)  implies  that  the  convective  heat  transfer 
can  be  significantly  reduced  by  decreasing  the  pressure.  We  utilize  this 
behavior  in  a subsequent  section  to  scale  a low-pressure  experiment. 


The  radiative  heat  transfer  losses  are  given  by 

q = e ct  T 4 
'r  x 


(287) 


where  the  emissivity  can  approach  unity  for  gases  under  high  temperature 
conditions.  Typical  chemical  rocket  chambers  experience  75%  of  their 
losses  via  convective  heat  transfer  and  25%  via  radiative  heat  transfer, 
but  they  operate  above  atmospheric  pressure.  Subatmospher ic  pressure 
operation  would  reduce  qc  and  therefore  make  q^  dominant.  Effective 
chamber  design  might  be  employed  to  capture  some  of  the  radiant  heat 
loss . 
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Other  losses  to  the  system  arise  due  to  heating  and  flow  nonuniform- 
ities. Heating  nonuniformities  arise  due  to  nonuniform  absorption  of  the 
laser  radiation  leading  to  "hot  spots"  that  must  "mix"  with  relatively  cold 
propellant.  Such  mixing  will  more  than  likely  be  turbulent  and  it  will  be 
entropy  producing.  These  thermal  effects  must  be  considered  in  addition 
to  the  nonuniform  velocity  effects  produced  by  the  boundary  layer  and  heat 
transfer  losses  during  the  nozzle  expansion  processes. 

All  of  the  effects  discussed  above  lead  to  a reduction  of  the  equi- 
librium thermodynamic  efficiencies  described  in  the  preceding  subsections. 


T'”' 
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3.  FLUID  DYNAMICS 


The  absorption  of  laser  radiation  by  a fluid  whether  stationary  or 
in  motion  leads  to  the  dynamic  change  of  that  fluid  medium.  In  general, 
laser  absorption  fluid  mechanics  is  very  complex;  phenomenon  may  be 
steady  or  unsteady,  explosive  or  simply  due  to  volumetric  heating. 

The  physics  of  any  particular  scheme  is  strongly  governed  by  the 
absorption  process.  Laser-supported  waves,  whether  they  are  detonation 
waves  or  deflagration  waves,  are  closely  linked  to  the  absorption  process 
Furthermore,  the  structure  of  the  steady  wave  phenomenon  is  likewise 
dominated  by  the  absorption  of  laser  radiation.  A gas  may  absorb  laser 
radiation  at  a specific  wavelength  via  molecular  absorption  in  some 
temperature  range.  With  increasing  temperature  the  absorption  co- 
efficient typically  decreases  as  the  absorbing  specie  dissociates  until 
inverse  bremstrahlung  leads  to  a rapid  increase  of  absorption.  The 
transition  between  the  two  types  of  absorption  will  have  an  important  effect 
on  the  type  of  absorption  wave  produced. 

In  the  sections  that  follow,  one -dimens ional  laser  absorption  wave 

phenomenon,  i.  e.  , the  complete*  Hugoniot  curve,  and  the  structure  of 

quasi -one -dimensional  absorption  zones,  are  calculated. 

a.  Hugoniot  Curve  for  One -Dimensional  Laser -Suppo rted  Waves 
(5) 

Raize  r obtained  a relation  for  the  laser -suppo  rted  detonation 

wave  velocity  under  the  "zero  counter  pressure"  limit.  His  expression  is 
apparently  only  valid  in  the  upper  (detonation)  branch  of  the  Hugoniot  curve 

jfi 

The  word  complete  refers  to  the  calculation  of  both  the  upper  (detonation) 
and  lower  (deflagration)  branches  of  the  Hugoniot  curve  for  one -dimens ional 
flow  with  energy  addition  through  laser  radiation.  Raizer  (Ref.  5)  com- 
puted the  detonation  wave  branch  in  his  classic  paper. 

(5)  Raizer,  Yu.  P.  , (1965)  "Heating  of  a Gas  by  a Powerful  Light  Pulse,  " 
Soviet  JETP,  2\_,  No.  5,  1009. 
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relating  pressure  to  density.  If  one  attempts  to  evaluate  his  relation  on 

the  lower  branch  of  the  Hugoniot,  one  obtains  a square  root  singularity. 

(5 ) 

This  difficulty  will  be  highlighted  when  we  reexamine  the  Raizer  analysis 
in  order  to  obtain  an  equation  which  admits  a solution  on  the  deflagration 
b ranch. 

The  one -dimensional  gasdynamic  equations  with  laser  energy  ab- 
sorption read  as  follows  for  steady  flow  (where  ( ) indicates  differentia- 

tion). p is  the  density,  u is  the  velocity,  I is  the  radiation  intensity  (W/cm"). 
and  a is  the  absorption  coefficient. 

(pu)x  = 0 ( 288 ) 

(pu)  ux  = -px  <28q> 

+ 1=0  (2‘H1) 

x 

laser  radiation  parallel 
to  flow 

(291) 

laser  radiation  antiparallel 
to  flow 

(292) 


[(Pu)  (h  + u2/2)] 


- 1 


I = 6 a I where  6 = 

x 


+ 1 


h = h (p.  p) 


P = Z (T.p)  p RqT 


(293) 


a = a (p.  T) 


(2'M) 
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The  implicit  approximation  made  in  writing  down  these  equations 
is  that  we  are  excluding  diffusive  effects.  Diffusive  effects  are  also 
neglected  when  deriving  the  shock  adiabatic  and  the  Hugoniot  curve  when 
chemical  energy  is  released  into  the  flow. 

Suppose  there  is  a discontinuity  zone  somewhere  in  the  flow.  This 
is  depicted  in  Figure  42.  Conditions  on  either  side  of  the  zone  are  ex- 
amined. The  frame  of  reference  considered  makes  the  discontinuity  zone 
(DZ)  fixed  in  space.  Across  DZ  from  Eq.  (288) 


Pu  = constant  = m 


(295) 


P)ul  =P2U2 


(296) 


Furthermore,  utilizing  Eq.  (295)  in  Eq.  (289)  provides 


Pj  + nnij  = p^  + mu 


(297) 


across  DZ.  Finally,  the  energy  Eq.  (290)  integrated  and  evaluated  awa\ 
from  the  DZ  gives 


m (hj  + Uj2/2)  + Ij  = m (h ^ + U2V2>  + 


(298) 


The  light  intensity  on  one  side  of  the  DZ  is  known  (i.  e.  . either  Ij  or  I^); 
the  intensity  on  the  other  side  of  the  DZ  depends  upon  the  absorption  pro- 
cesses In  turn  the  absorption  physics  (absorption  coefficient)  is  coupled 
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to  the  dynamics  (flowfield)  as  well  as  the  thermodynamics.  We  now  con- 
sider Eq.  (291).  For  the  case  with  6 = -1  and  taking  the  laser  flux  Iq  = Ij 


whe  re 


(299) 


(300) 


and  the  integral  is  taken  between  region  [ 1]  and  [ Z] . 

For  the  antiparallel  situation  (6  = 1),  - Iq  and 

h - I,  - ‘o  (' 


(301) 


If  we  consider  the  parallel  situation  (6  = -1)  with  the  direction  of  the  laser 
radiation  the  same  as  the  flow  direction,  then  Eq.  (298)  reads 

h2  + U2Z/Z  = hl  + Ul2/2  + 7^  n°2) 

for  the  complete  absorption  of  laser  radiation  in  DZ.  Region  1 is  the  un- 
perturbed state  with  p^.  p^,  hj  known.  The  wave  velocity  (of  DZ)  is 
Uj  = W.  The  Hugoniot  conditions  then  read 

P2u2  = pj  W = m (303) 


36 


p2  + m = pj  + m W 


(304) 


h2  + «22/2  = hj  + W2/2  + I0/p,  W 


(305) 


In  order  to  proceed  further,  the  perfect  gas  relation  with  constant  y is 
employed. 

h = ...X.  E. 
h v - 1 P 


(306) 


Raizer  ^ used  the  same  relationship  in  his  analysis. 

It  is  more  convenient  to  manipulate  the  equations  alter  defining 
new  variables.  Let 


" a pi/p; 


(307) 


Equation  (304)  reads 


P = P2/Pi 


(P-1) 


(300) 


Substitution  of  Eq.  (300)  into  the  square  of  Eq.  (205)  after  using  Eq  (306) 
provides  the  following  relationship  between  P-1  = t,  and  (q-l)/r|  o: 
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£ (r-1  + (j/2)c  + C ZTo  (r  -1  + a/2) 


(31 


2 2 2 
+ 5 rc  a + hq  <7=0 


where 


r = Y/(Y-n  and  H02  = ro2  pi/pi3 


(31 


Equation  (310)  is  a cubic  equation  for  C,  (pressure)  as  a function  of  the  o 

(5) 

(density).  This  relation  is  more  general  than  the  Raizer  - derived  re - 

t 

lation.  The  Raizer  formula  is  obtained  from  our  basic  equations  by  substi- 

J 

tuting  Eq.  (309)  into  Eq.  (306)  (without  squaring);  he  obtains 


Raizer  omitted  the  first  term  in  the  numerator  by  making  the  "zero  counter 
pressure"  approximation.  We  include  this  term  for  completeness.  The 
difficulty  in  Eq.  (312)  is  the  square  root  term  (1-q)^^.  On  the  detonation 
branch  with  n < 1 (Pj  < P^'  an  acceptable  value  for  the  pressure  can  be 
obtained;  however,  on  the  deflagration  branch  n > 1 (P  j > would  be 

forced  to  take  accept  an  imaginary  pressure  if  we  used  Eq.  (312).  The 
use  of  our  Eq.  (3101  circumvents  this  difficulty  and  permits  the  calculation 
of  the  deflagration  branch. 

Prior  to  describing  solutions  for  the  generalized  Hugoniot  curve  we 
must  specify  the  physical  limits  of  the  parameters. 


(i)  The  density  is  a positive  finite  quantity:  0 < q < *> 

(ii)  The  pressure  is  a positive  finite  quantity:  0 < P < °o 

( iii)  A solution  for  Eq.  (3101  is  given  by  £ = 0,  a = 0 independent 

2 

of  the  value  of  H . This  point  corresponds  to  W =0  since 
Pl  = p^,  and  pj  = p^  for  this  situation. 

(iv)  The  Hugoniot  curve  exists  in  the  second  and  fourth  quadrants 
of  the  coordinate  system  formed  by  the  intersection  of  the 

£ = 0 line  and  q = 1 line.  This  conclusion  is  reached  from 
Eq.  (309)  which  states  that  all  solution  points  on  the  Hugoniot 
curve  must  have  negative  slopes  in  the  £ -q  plane. 

The  forbidden  regions  are  shown  shadded  in  Figure  43  along  with 
solutions  curves  obtained  by  solving  Eq.  (310).  The  Hugoniot  curves  on 
the  deflagration  branch  lead  to  negative  pressures  at  values  of  q that 
depend  upon  the  laser  radiation  intensity.  Note  how  rapidly  negative  pres- 
sures are  realized  for  H = 10  whereas  a gentle  approach  at  nearly  con- 
stant pressure  is  found  for  Ho  = 1.  The  variation  of  the  value  of  q whence 
the  pressure  first  becomes  negative,  i.  e.  , d{.r  can  be  found  vs  Hq  from 
Eq.  (310)  by  setting  £ = -1.  We  obtain  a parabolic  equation  for  a = 

(q  _ l)/q  . This  algebraic  equation  can  be  readily  solved.  In  particular 
for  y = 1.  40, 


acr  = 5/6+  Hq2/18  - \j (Hq4  + 30  Hq2)/324 


[313) 


A plot  of  q vs  Hq  is  shown  in  Figure  44.  Note  that  the  slowest  variation 

of  n with  H occurs  in  the  neighborhood  of  H 1 . Furthermore,  as 
1 c r o o 

H becomes  excessively  large  q -*  i and  as  H becomes  very  small 
o cr  o 


cr 


6. 
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An  important  concern  for  laser  propulsion  is  the  wave  propagation 

2 

velocity  W.  Since  \p^/pj  = a^  , the  undisturbed  speed  of  sound  squared, 
Eq.  (309)  becomes 

\ W2/a,2  = (P  - 1)/(1  - q)  (314) 


We  plot  this  velocity  ratio  for  our  solution  curves  as  shown  in  Figure  45. 
From  Eq.  (314),  we  note  that  W will  be  sonic  (equal  to  ap  when  (P  - 1)/ 

(1  - q ) is  y . This  line  is  shown  on  the  same  figure.  Note  that  there  are 
no  subsonic  solutions  for  large  values  of  the  parameter  Hq  (=  10).  For 
Hq  = 1 there  is  a small  regime  of  subsonic  waves  prior  to  the  negative 

pressure  limit.  As  one  further  reduces  the  power  number  (H  = l/lO), 
there  is  an  extensive  region  of  subsonic  wave  solutions.  For  the  range 
of  q 's  shown,  the  wave  speed  becomes  as  small  as  one-half  aj. 
b.  Absorption  Zone  Gasdynamics 

Consider  the  basic  gasdynamic  equations  with  laser  energy  addition. 
These  equations  read: 


TlP"s'=  0 

( 115) 

dp  , du  „ 

(316) 

!♦»■£)*  -x 

d (IA) 
dx 

(317) 

^ ■ - * • 

(318) 

h = h (p,  p) 

(319) 

T = T (h,  p) 

(32  0) 

H = h + } u2 

(321) 

where  A is  the  area  variation  of  the  channel. 
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In  order  to  integrate  these  equations,  we  transform  the  spatial 
coordinate  x to  the  optical  depth  coordinate  j using  the  transformation 

dT  s a dx 

The  preceding  equations  then  read: 


-XT  (p  u A)  = 0 


du 

dT 


/ dh  , du  \ 

+ u if] 

d(IA)  _ 


1 d (LA) 
A 


dr 


dr 


IA 


(322) 

(32  3) 
(324) 

( 325) 
(326) 


Equations  (319),  (320)  and  (321)  remain  unchanged.  All  the  dependenl 

variables  shall  be  presumed  known  at  the  original  station  Xj  or  r | The 

reference  quantities  will  be  denoted  by  Uj,  p,  , Tj,  h^ , Hj  = h ^ 

2 

+ u.  /2,  A^  and  Ij.  Note  that  the  independent  variable  T is  dimension- 
less. If  we  now  choose  to  define  dimensionless  variables  utilizing  our 
reference  quantities  and  denote  them  by  ( ) = ( )/(  )j,  then  the  equations 

read: 

d 


dr 


(p'  u ' A ’ ) = 0 


*2l 

dr 


P « 


du 

dT 


= 0 


p*  u'*‘Sf  = - p,  d <■*  a>!) 

H dT  1 


~Xt~ 


II  = 


('  ■ ^t)  " * (-wr) u 


2 


d = -r  a 

dT 


(32  7) 
(328) 

( 32^) 

(3  <0) 

(331) 
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'.r  ' ■ — r* 3 — " 7 1 1 


wmm 


It 


where  a compressibility  number 


a power  ratio 


„ 2 / 2 

C 1 ^ Pl/Pl  U1 


Pj  — Ij  A j/P]  Uj  Aj  Hj 


(3321 


( 333) 


* * A*  , 

p u A =1 


and  the  ratio  (Uj  /2Hj)  contained  in  Eq.  (350)  characterizes  the  laser  energy 
addition  problem. 

Equation  (327)  simply  integrates  to 

(3  34) 

The  transformation  to  the  optical  depth  coordinate  T makes  it  easy 
to  integrate  Eq.  (331).  The  solution  reads 

* * (T,  -T) 

I*  A*  = e 1 (335) 

The  use  of  Eqs.  (334)  and  (335)  in  Eq.  (329)  permits  the  integration  of  the 
equation.  The  solution  for  the  stagnation  enthalpy  becomes 

) (336) 


* T 1 'T 
H = 1 + Pj  (1  - e 


We  also  have  a thermodynamic  relationship  (319)  which  may  in  general  be 
a complex  function 

h = h (p  , p ) (337) 

It  now  remains  to  solve  the  momentum  equation  relating  the  variation  of 
pressure  to  the  velocity  variation.  We  begin  by  multiplying  Eq.  (328)  by 


A and  utilizing  Eq.  (334)  to  find 


du  * dA  d (A  p ) 

■ ■ ■■■  — n ■ ■ " — - 1 


c. 


dr 


dT 


dT 


( 338) 


We  can  integrate  Eq.  (338)  to  the  implicit  form 

1 Sic  \ i 

U = F(T)  - F(Tj)  + -J-  Uj  -A  P + Aj  Pj 

C1  c, 

45 


( 339) 


where 


F(T)  - F(t  ,)  = 


* , . dA  , 

P (r)  dr 


(340) 


The  system  of  equations  can  be  solved  algebraiclly  for  some  idealized  situ- 
ations. In  general,  one  must  resort  to  a numerical  calculation  in  order  to 
solve  the  system. 

c.  Constant  Pressure  Approximation 

From  Eq.  (328)  note  that  if  the  flow  is  nearly  incompressible  (Cj^  » 1), 
p « = constant.  Upon  making  this  approximation,  one  can  solve  Eqs.  (334), 

(336)  and  (337)  simultaneously  or  we  can  obtain  a single  equation  for  the 
density  hi  Equation  (336)  and  the  definition  of  H by  Eq.  (330)  implies 

(>-^) h%  W u*2*  i+p>  (■  -*T,’T) 


Equation  (334)  gives 


(’  ‘ 2»l)  " + (™l)  „<2’ 


= 1 + P 


p*  A* 


On  the  initial  condition  isobar  (pj  = 1)  in  relation  (337)  we  have 


h’:  = hj  ' (p^*  = 1,  p*).  Hence 

2 


(‘  ‘ 2«i) 


»*!*  (Pi*  = 1.  P*>  + 


(Jill')  __i 

\2Hi/  *2 


= 1 + P, 


Tl-^ 


p*  A* 


For  a perfect  gas. 


or  in  dimensionless  form 


h = -X-  £. 

y - 1 p 


h*  = p*/p* 


(341) 


(342) 


(34  3) 


and  for  the  pj  ' = 1 isobar, 


* >i' 

h i = P 


,.-l 


.,.-1 


Therefore  using  Eq.  (343)  in  Eq.  (341)  and  defining  R = p , we  get  a 


quadratic  relation  for  R , i.  e.  , 

a R*2  + b R + c = 0 

where 


(34 


. pJ)  _l 
\2IV  A* 


u. 


b = 1 


2H . 


- c = 1 + P 


The  solution  is  therefore 


R*  = 


♦ J 


b - 4 a c 


2a 


(34 


We  further  note  that  the  ratio  of  Uj  /2Hj  will  be  small  for  our  low  speed 
high  pressure  environment.  Consequently  a » 0,  b « 1.  Then  Eq.  (344) 
simply  gives 

R*  = 1 + P1  ^1  - eTl  ^ (34 


or 


* 

P = 


1 


(34 


1 + P,  (l  - e 1 ) 

The  selection  of  the  perfect  gas  law  made  it  easy  to  solve  Eq.  (342) 


If  we  do  not  make  this  approximation,  but  let  a and  b be  zero  and  one 


respectively  (equivalent  to  zero  velocity),  then  the  equation  for  the  density 


(348) 


hj  (1,  p ) = 1 + Pj  (1  - e 


Tl-T 


A curve  fit  for  (1,  p ) e.  g.  , from  a Mollier  diagram  could  be  used  to 

obtain  an  explicit  analytical  relationship  for  p . The  inversion  of  Eq.  (348) 
gives,  in  general, 

T,  -T 


* * 

p = p 


(349) 


* j^l  + Pj  (1  - eTl  T)j 

whence  p is  a function  of  the  argument  in  brackets. 

At  constant  pressure  we  can  obtain  an  approximate  form  for  the 

J{C 

enthalpy  variation  with  density,  i.  e.  , hj  (1,  p ) under  nonideal  conditions. 
We  no  this  for  high-temperature  water  - both  as  an  example  and  as  a poten- 
tially important  practical  case.  An  examination  of  the  H^O  Mollier  chart 
shows  that  at  a high  pressure  (50  atm),  the  variation  of  enthalpy  is  approxi- 
mately 


h,  (1,  p ) = (p  ) 


,=V5/4 


(350) 


in  the  temperature  range  of  3500°K  to  5500°K.  Expression  (350)  differs 
from  Eq.  (343)  because  of  dissociation  effects.  Calculations  for  the  ideal 
gas  Eq.  (343)  and  the  real  gas  Eq.  (350)  will  be  described  in  the  subsequent 
section  on  the  molecular  absorption  of  AlO. 
d.  Molecular  Absorption  of  AlO 

The  absorption  coefficient  of  1%  Al^O^  molar  in  H^O  at  a pressure 
of  50  atm  peaks  at  an  interesting  (from  a mission  point  of  view)  temperature 
and  value  as  calculated  by  Douglas -Hamilton.  ^ The  curves  obtained  may 
be  fitted  quite  well  to  be  parabolas  in  the  semi-log  plot  of  a (cm  S vs 
T(°K).  A curve  fit  of  the  50  atm  case  was  made  (c.  f.  Figure  46)  and  the 


(6)  Douglas-Hamilton,  D.  H.  . AERL  private  communication  (1976). 
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resulting  relation  reads 


In  Q = - 9.  24  x 10"7  (5000  - T)2  - 0.  916 


(351) 


The  gas  law  may  be  written: 

T [°K]  = °-  2*9  P fatm| 

Z p [gm/ cm3) 

The  compressibility  factor  is  obtained  from  the  Mollier  chart  as: 

Z = 0.  02 1 3/p2/ 3 

for  p = 50  atm.  Hence,  the  temperature  is  related  to  the  density  as 

T = 514. l/p1/3 


(352) 


( 35  3) 


(354) 


Substitution  of  Eq.  (354)  into  Eq.  (351)  then  provides  an  explicit  relationship 
between  the  absorption  coefficient  and  the  dimensionless  density  p , i.  o.  , 


a = exp 


- 9.  24  x 1 


0‘7  ^5000 


3500 


) 


0.  016 


; 35  5) 


where  the  reference  state  (Tj)  has  been  chosen  to  be  3500°K.  (This  state 

corresponds  approximately  to  the  combustion  temperature  of  a stochiometr ic 

mixture  of  gaseous  H7  and  gaseous  at  50  atm,  c.  f.  Section  V.  5.  ) 

The  molecular  absorption  coefficient  (355)  was  used  in  conjunction 

with  the  gasdynamic  equations  in  order  to  determine  the  structure  of  the 

absorption  zone.  Both  an  ideal  equation  of  state  (h*  ~ p ')  and  molecular 

-5 /4 

H^O  Mollier  chart  - derived  equation  of  state  (h*  ~ p*  " ' ) was  utilized  to 

integrate  the  equations.  A particular  case  of  constant  area  (A*  = 1)  was 
selected.  Hence  the  solution  for  the  intensity  is  merely  a decaying  expo- 
nential in  the  r plane  as  given  by  Eq.  (335).  This  is  plotted  in  Figure  47; 
note  that  i/lj  is  1%  of  its  initial  value  at  T » 4.6. 

The  calculations  presented  below  correspond  to  the  dimensionless 
compressibility  number,  Cj  = 8.  66  and  power  ratio  Pj  = 2.9,  23.7.  Since 


Figure  47  Density  and  Intensity  Variation  in  Optical  Depth  Plane 
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Cj  = y /M,  where  M is  the  Mach  number,  increasingly  massive  amounts 

of  heat  at  low  subsonic  Mach  numbers  are  added  as  increases.  The 

density  is  also  shown  in  Figure  47.  Note  that  most  of  the  variation  with  t 

ceases  after  r > 2.  Also  note  the  sensitivity  to  ideal  (h*  ~ n*  ~S  as 

-5/4 

opposed  to  molecular  (h*  ~ p*  ' ) gas  behavior.  In  Figure  48  we  illustrate 

the  temperature  variations  under  the  same  conditions.  Since  the  peak  ab- 
sorption at  50  atm  is  near  5000°K  (see  Figure  46)  it  will  become  evident 
that  absorption  behavior  in  the  zone  is  strongly  influenced  by  the  tempera- 
ture variations.  One  observes  that  when  the  heating  is  very  strong  (P ^ » 1), 
the  temperature  exceeds  5000°K  very  rapidly  leading  to  a turn-off  of 
absorption  as  shown  in  Figure  49.  It  is  interesting  to  note  the  change  in 
character  between  the  Pj  ~ 2.  9 case  and  the  Pj  = 23.  7 case.  Physically 
one  would  prefer  to  observe  the  variation  of  the  absorption  coefficient  with 
x (cm).  The  result  of  transforming  our  calculations  from  the  t -plane  to 
the  physical  plane  are  shown  in  Figure  50.  At  the  high  intensity  levels,  the 
absorption  of  laser  radiation  is  virtually  complete  after  a few  centimeters. 
On  the  contrary,  tens  of  centimeters  are  needed  to  absorb  the  radiation  at 
the  low  intensity  levels.  This  apparent  sensitivity  of  q to  P^  has  signifi- 
cant implications  in  the  design  of  an  experimental  engine  because  the 
"physical"  absorption  length  may  become  excessively  large, 
e.  Absorption  Stability 

The  stability  of  laser -heated  flows  is  an  important  concern  for  any 

(7) 

absorption  scheme.  Wu  and  Pirri  examined  the  "local"  stability  of  a 
laser  heated  flow  in  a geometry  representative  of  a typical  rocket  motor 


(7)  Wu,  P.  and  Pirri,  A.  N.  , AIAA  J.  14,  No.  3,  390(1976). 
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Figure  50  Absorption  Coefficient  in  Physical  Plane 


with  an  absorption  chamber  followed  by  a throat  and  then  an  expansion 
nozzle.  rhe  absorption  coefficient  was  modeled  as  K pn  I'1"  where  Kt 
is  a constant  and  the  Inverse  Bremsstrahlung  values  n = 2.  m - - '■>  1 
are  taken. 

They  obtained  curves  of  neutral  stability  tor  a spe.  itu  nozzle  con- 
figuration and  laser  power  input.  Ilieir  unstable  regions  .ere  found  to 
be  upstream  of  the  throat  in  the  subsonic  portion  of  the  flow.  Plie  analysis 
presented  cannot,  however,  answer  the  question  of  how  an  unstable  dis 
turbance  behaves  as  it  travels  to  a new  location  in  the  flow.  Nevertheless, 
one  is  cautioned  to  examine  the  stability  of  laser  absorption  zone  heating 
prior  to  designing  a CW  engine. 

The  favorable  aspects  of  the  absorption  behavior  of  AlO  to  10.  b pm 
radiation  (Figure  SO)  were  described  previously;  but  the  stability  aspects 
of  the  heated  zone  have  not  been  treated.  It  is  not  possible  to  apply  the  Wu 
and  Pirri^'  stability  calculations  to  the  molecular  absorption  situation. 
Since  the  stability  of  the  laser-heated  zone  as  treated  by  \V u and  Pirn^ 
is  dependent  upon  so  many  aspects  of  the  physical  problem  such  as  absorp- 
tion mechanism,  channel  area  variation,  real  gas  phenomena,  beam  direc- 
tion (relative  to  flow  direction)  etc.  , it  would  appear  to  be  fruitful  to  addres 
the  stability  question  from  a more  global  point  of  view.  Iliat  is,  one  could 
hope  to  answer  the  question:  What  type  of  absorption  coefficient  behavior 

in  the  flowing  environment  would  lead  to  complete  stability'*  Such  an  invest i 

(8) 

gation  will  be  patterned  after  the  work  of  Jacob  and  Mani.  ' 


(8)  Jacob,  J.  and  Mani,  S.  "Thermal  Instability  in  High  Power  Laser 
Discharges"  (11>7S). 


Ik  l 


I 


■RB 


4.  SCALING  STUDIES  FOR  TWO-PORT  CW  ENGINE 

In  the  past  experimental  CW  laser  propulsion  engines  have  been  very 
small  due  in  part  to  the  laek  of  sufficient  laser  power  available.  We  will 
confirm  these  features  of  laser  rockets  in  detail  for  various  rockets  charac- 
terized by  a constant  ratio  of  specific  heats, 
a.  Nozzle  Throat  Size 

The  cross-sectional  area  of  the  sonic  nozzle  throat  following  the 
laser  absorption  section  is  given  by  the  following  expression: 


Y5 


t 1 


A6  = 


- 1 


( 3 5 1> ) 


where  sub-6  refers  to  the  nozzle  throat  and  sub-5  refers  to  the  absorption 
chamber:  rh  is  the  mass  flow,  Tq  is  the  stagnation  temperature,  pQ  is  tin* 
stagnation  pressure,  R^  is  the  gas  constant,  W is  the  molecular  weight, 
and  y is  the  ratio  of  specific  heats. 

Under  ideal  conditions  the  maximum  attainable  exit  velocity  is  given 


as: 


2 2 V5  Ro  ... 

U7  - y5  - 1 W5  1 5 


(357) 


The  specific  impulse  is  simply  given  by  u^/g,  hence,  a mission-derived 
value  for  I defines  the  right  hand  side  of  Eq.  (357).  A relation  betweer 
T,.(0K)  and  Igp(sec)  is  given  below  where  W^  is  given  in  gm/gm-mole. 


Ts(  K)  = 5.  78  x 10 


W)*v;. 


(358) 


substitution  of  Eq.  (358)  into  Eq.  (356)  provides: 

m 8 I 

a / = — nY,) 
6 Pq5  5 


(359) 


where 


W) 


V5+  1 

v5  - hi/2 


(360) 


r(v5)  varies  monotonically  from  0.  339  (y,.  = 1.1)  to  0.616  (y^  = 1 . 67). 

The  primary  variation  of  in  Eq.  (359)  is  therefore  due  to  rh  lSp/PQ5- 

If  we  select  an  I of  1000  sec,  then 
sp 

A,  (cm2)  = 0.967 
6'  p05(atm) 

This  relation  is  plotted  in  Figure  51  with  the  nozzle  throat  area  as  the 
parameter.  In  addition,  we  have  added  the  laser  power  scale  on  the  left 
hand  side  of  the  figure  assuming  all  of  the  laser  power  goes  into  heating 
the  gas,  i.  e.  , 


■ **  - 5 — y 

g I 

p sp 

where  IP  is  the  laser  power.  We  have  gone  through  this  relatively  simple 
exercise  to  highlight  the  observation  that  low  pressure  operation  is  favorable 
from  the  point  of  view  of  increasing  the  geometry  of  candidate  CW  engine 
configurations.  For  example,  a laser  power  of  10  kW  absorbed  into  a flow- 
ing gas  at  0.  2 gm/sec  and  10  * atm  pressure  would  have  a throat  area  of 
2 

7.  5 cm  or  a throat  radius  of  1 . 5 cm.  The  remainder  of  the  scaled-up 
engine  would  be  a nice  laboratory  device,  i.  e.  , a nozzle  area  ratio  of  100 
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Figure  51  Throat  Area  for  1 = 1000  sec 
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would  provide  an  exit  radius  of  15  cm.  The  energy  addition  in  the  chamber 
would  yield  similar  sizes.  Subatmospher ic  energy  addition  would  be  practi- 
cal at  altitude  and  of  course,  in  exoatmospheric  applications.  Subatmo- 
spheric  experiments  would  imply  operation  in  a partially  evacuated  chamber 
if  all  of  the  other  aspects  of  the  simulation  appear  feasible.  Iwo  key 
aspects  are  the  rocket  heat  transfer  and  the  absorption  of  the  laser  rad- 
iation. 

(4) 

Based  upon  experience  with  rocket  motors  Summerfield  one  can 
readily  show  that  the  convective  heat  transfer  rate  varies  as  follows  with 
pressure: 

q ~ p5  l'^b 

For  a reduction  of  chamber  pressure  from  10  atm  to  10  atm  or  a factor 
of  f00  reduction,  the  convective  heat  transfer  is  reduced  by  nearly  two 
orders  of  magnitude.  Hence,  the  heat  transfer  problem  is  diminished  at 
low  pressures. 

The  situation  with  the  absorption  phenomena  is  not  as  favorable:  in 
fact,  most  gases  have  decreased  absorption  coefficients  at  low  pressures. 

A typical  molecular  absorber  of  10.  6 pm  radiation  (AlO)  behaves  in  this 
manner.  Inverse  Bremsstrahlung  absorption  is  also  decreased  at  low 
pressures.  This  may  imply  the  use  of  very  long  absorption  chambers  in 
low  pressure  rockets  - but  this  is  not  a critical  concern  for  the  proof  of 
principle  simulation  suggested  in  Section  V.  5. 

b.  Two-Port  Rocket  Scaling  Considerations 

The  analysis  of  the  two-port  rocket  considers  nine  different  station 
locations.  These  stations  are  defined  in  the  schematic  two-port  engine  in 
Figure  52  and  described  in  Table  4.  One  can  calculate  the  relative  sizes 
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Figure  52  Schematic  Two-Port  CW  Laser  Engine 


TABLE  4.  STATIONS  IN  TWO-PORT  LASER  ENGINE  ROCKET 


of  such  an  engine  once  a specific  mission  requirement  has  been  established. 

A specific  impulse  of  1000  sec  is  chosen.  This  corresponds  to  an  interest- 
ing mission  as  defined  in  the  trajectories  section.  The  maps  developed  will 
provide  indications  of  the  physical  size  of  two-port  engines  as  a function  of 
laser  power  level. 

The  maps  are  derived  by  considering  the  one  dimensional  gasdynamic 
heat  addition  problem  using  gas  tables.  Constant-y  calculations  are  per- 
formed even  though  the  ratio  of  specific  heats  varies  considerably  for  the 
thermodynamic  states  of  interest.  (See  Appendix  C). 

As  configured  in  Figure  52,  the  engine  consists  of  three  key  elements 
that  provide  for  the  sizing  of  the  engine;  aerodynamic  window,  absorption 
chamber,  and  rocket  nozzle.  The  nozzle  throat  sets  the  "level"  of  the 
rocket  geometry  - and  thereby  is  a key  aspect  of  the  scaling.  The  elements 
are  discussed  in  turn. 

c.  Aerodynamic  Window 

In  Figure  53  we  see  an  exploded  view  of  a simply  wedge  supersonic 
aerodynamic  window.  The  interior  wedge  angle  Q is  chosen  such  that  the 
flow  does  not  turn  upstream  (antiparallel)  to  the  incoming  laser  radiation. 
This  aerodynamic  window  must  operate  under  pressure  conditions  that  vary 
from  1 atm  at  sea  level  to  essentially  vacuum  when  operating  exoatmospheri- 
cally.  The  maximum  flow  deflection  angle  u>max-C0(degrees)  is  given  by  an 
infinite  pressure  ratio  (for  the  vacuum  case).  The  variation  of  -u>  and 

6 is  given  in  Figure  54  for  various  values  of  and  the  y of  the  window  flow 
gas.  Note  that  we  have  indicated  the  minimum  required  to  achieve  the 
desired  turning  of  the  flow  for  each  y - this  corresponds  to  0 - 0°.  The 
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significant  variation  of  the  minimum  window  flow  Mach  number  vari- 

ation with  y is  dramatic  and  is  shown  in  Figure  55.  This  figure  implies 
that  a very  low  Mach  number  can  be  utilized  for  the  flow  of  a monatomic 
gas  like  for  example.  Low  has  immediate  implications  with  respect 
to  the  geometry  between  station  1 and  station  2.  These  details  are  con- 
sidered below.  Shown  on  Figure  56  is  the  area  ratio  between  station  1 and 
station  2 for  various  values  of  and  y.  The  use  of  monatomic  gases 
(y  = 1 . 67),  would  imply  the  use  of  an  ~ 2 flow  and  an  area  ratio  of  ~ 2. 
A diatomic  gas  such  as  would  imply  the  use  of  an  M^  ~ 3 with  an  area 
ratio  of  ~ 5.  Furthermore,  a polyatomic  gas  such  as  H^O  (with  a low- 
temperature  y = 4/3)  would  require  a window  a 4 and  an  area  ratio  of 
> 10.  If  the  window  throat  is  completely  "open,  " then  the  relationship 
between  A^/Aj  and  dj/d^  is  given  by 


2 = 1 'V1  ' 


(364) 


This  is  plotted  on  the  right  hand  side  of  Figure  56.  In  lieu  of  a completely 

"open"  ring  of  area  Aj  with  height  dj,  one  could  alternatively  think  of  the 

window  throat  as  made  up  of  N small  holes  of  radius  6j.  The  window  throat 

2 

area  A^  is  then  equivalent  to  Ntt6j. 

The  relationship  to  the  area  ratio  then  becomes 


N^a 


1/2 


1 


1 


(365) 


2 2 

This  relationship  is  also  plotted  on  Figure  56.  We  will  utilize  these  quantities 
shortly  when  we  specify  the  geometry  of  a particular  engine.  Previous  rough 
scaling  efforts  have  shown  that  the  smallest  dimension  of  the  rocket  config- 
uration shown  will  be  the  window  throat  dimension  - d^  or  6j . 
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The  window  area  which  permits  the  entrance  of  the  laser  beam  into 

the  engine  is  defined  by  the  radius  d^  as  shown  on  the  schematic  Figure  53. 

This  area  is  sized  by  the  breakdown  limits  of  the  laser  radiation.  A good 

3 2 

number  for  "dirty  air"  breakdown  is  10  kW/cm  . The  breakdown  limit 
defines  a minimum  d^  for  a particular  power  level.  Breakdown  limitations 
reveal  that  a practical  engine  would  be  meters  in  diameters  whereas  a 
laboratory  device  would  be  centimeter  size, 
d.  Absorption  Chamber 

The  absorption  chamber  is  between  the  rocket  nozzle  and  the  aero- 
dynamic window.  Figure  57  gives  the  area  ratio  A^/A^  = A^/A^  as  a 
function  of  the  window  Mach  number  M^  and  6.  Two  other  inputs  needed  to 
generate  this  figure  are  the  absorption  chamber  initial  Mach  number,  taken 
as  M^  = 1/10,  and  the  ratio  of  stagnation  temperatures  before  and  after 
absorption,  i.  e.  , T^/Tq^.  Note  how  insensitive  this  ratio  A^/A^  is  to 
M^  and  y A value  of  between  2/10  and  6/10  is  reasonable. 

The  region  between  station  3,  the  subsonic  flow  downstream  of 
normal  shock,  and  station  4,  the  preabsorption  M^  = 1 / 1 0 flow,  is  a sub- 
sonic diffuser  with  a slight  inclination  (rs  6°)  to  avoid  flow  separation  and 
the  concomittant  losses.  The  subsonic  diffuser  separates  the  two  constant 
area  sections:  supersonic  diffuser  section  and  the  absorption  section. 

The  addition  of  heat  in  a constant  area  duct  neglecting  friction  (Rayleigh 
line)  leads  to  the  determination  of  the  ratio  A5/A^  as  shown  in  Figure  58. 
Note  the  slow  variation  with  y for  the  various  levels  of  T^/Tq^.  Further- 
more, the  area  ratio  between  the  absorption  chamber  and  the  nozzle  throat 

is  < 3 for  all  values  of  Tr.c/Tr. . > 4. 
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Figure  57  Subsonic  Diffuser  Dimensions 
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e.  Expansion  Nozzle 

It  is  evident  from  Figures  57  and  58  that  the  nozzle  throat  scales 
the  geometry  of  the  absorption  region  through  the  nozzle  throat  area  A^ . 

An  expression  for  A^  was  developed  previously  (Eq.  (359))  which  developed 
a map  (Figure  51)  that  highlighted  the  advantages  of  low  pressure  operation. 
The  parameters  that  fix  A^  are  Igp,  IP,  pQ5,  and  the  propellant.  The 
propellant  influence  is  subtle  since  it  appears  in  the  ratio  of  specific  heats 
as  well  as  in  our  assumption  (used  in  Eq.  (362))  that  all  of  the  laser  radi- 
ation is  absorbed  in  the  medium. 

The  pressure  difference  between  the  absorption  chamber  and  nozzle 
exit  plane  will  accelerate  the  propellant  gases  as  the  nozzle  exit  area  A?  is 
increased  from  A^.  The  variation  over  a large  range  of  pressure  ratios 
and  y are  shown  in  Figure  38.  A practical  range  of  area  ratio  A-,/A^  is 
between  10  and  10  ; this  implies  pressure  ratios  between  104  and  10°. 
Furthermore,  a significant  variation  of  A^/A^  with  y at  the  largest  area 
ratios  is  to  be  noted.  A further  connection  between  two  key  elements  of 
the  engine  is  provided  by  combining  Figure  56  and  58  to  provide  the  ratio 
of  throat  areas  A^/A,.  This  ratio  is  given  in  Figure  59.  It  varies  signifi- 
cantly with  y and  M^. 

f.  Scaling  Maps 

The  considerations  of  the  preceding  section  indicate  that  the  specifi- 
cation of  a "general  engine  concept"  would  necessitate  the  construction  of 
a multidimensional  map.  In  order  to  reduce  the  number  of  variables,  three 
interesting  cases  are  chosen.  These  cases  are  defined  in  Table  5.  Some 
remarks  relating  to  the  choice  of  various  parameters  are  now  given.  The 
three  cases  are  distinguished  primarily  by  the  selection  of  y;  based  upon 
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Figure  55  the  choice  of  y prescribes  a minimum  M^.  We  have  selected 
such  an  The  stagnation  temperature  ratio  allows  for  a significant 

range  of  temperature;  specifically,  if  preheating  is  necessary  to  avoid 
condensation,  T^j.  can  reach  10,  000°K.  The  area  ratios  in  Table  5 are 
then  taken  from  previously  developed  figures. 

All  three  cases  are  considered  to  develop  a specific  impulse  of  1000 
sec,  but  at  two  different  chamber  pressures:  p^^  =10  atm  and  p^,.  =1/10 
atm.  One  map  has  been  drawn  for  each  of  the  three  cases;  since  there  is 
a relatively  small  difference  between  the  cases  Tq,./T  ^ = 8 and  16,  we 
display  only  the  Tq,-/Tq^  = 8 case.  The  10  atm  situation  is  described:  the 
1 / 1 0 atm  on  the  maps  are  simply  shifted  to  the  right  by  two  orders  of  magni- 
tude. The  scaling  maps  are  given  in  Figures  60,  61  and  62.  Case  1 (y  = 

1.  20)  is  interesting  from  the  standpoint  that  the  window  throat  area  A^  is 
rather  small  compared  to  the  other  dimensions.  Case  2 (y  = 1.40)  brings 
the  dimensions  closer  together.  Finally,  Case  3 (y  = 1.67)  indicates  that 

all  the  dimensions  in  the  rocket  are  less  than  order  of  magnitude  apart. 

2 

We  also  show  the  minimum  window  area  A_  = rrd  which  results 

2 w 

(y  2 

from  the  use  of  10  W/cm  as  the  breakdown  limit. 
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Figure  62  Two-Port  Scaling  Map  (Case  3) 
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5.  experiments  for  c.w  laser  propulsion 

In  order  to  propel  a vehicle  efficiently  utilizing  laser  radiation  the 
absorption  mechanism  must  be  understood.  This  section  defines  a number 
of  experiments  that  could  potentially  aid  in  establishing  the  practicality  of 
CW  laser  propulsion. 

Absorption  phenomena  (I)  are  most  fundamental.  There  are  physical 
properties  of  gases  that  must  be  known  - to  as  yet  an  unspecified  tolerance  - 
in  order  to  design  a two- port  engine.  It  is  obvious  as  observed  in  Section 
V.3  that  the  theoretical  modeling  of  the  absorption  region  depends  strongly 
upon  the  absorption  coefficient  behavior.  Whether  the  absorption  zone  is 
one  diameter  in  width  or  a hundred  meters  in  width  is  an  important  concern. 
The  important  differences  between  molecular  and  Inverse  Bremsstrahlung 
are  also  very  important  to  chamber  design.  And  in  view  of  the  fact  that  a 

v 

can  have  a behavior  that  differs  in  character  for  an  individual  propellant  as 
well  as  for  a propellant  with  a modifying  seedant,  experimental  measure- 
ments of  the  absorption  coefficient  are  much  needed.  (The  following  section 
describes  Type  I experiments.) 

Absorption  gas  dynamics  (II)  involves  Type  I absorption  phenomena 
in  a flowing  environment.  The  key  question  from  a practical  point  of  view 
is  can  one  control  the  la ser- supported  waves;  a secondary  question  relates 
to  the  behavior  of  la  ser- supported  waves  as  one  modifies  the  absorption 
phenomenon  through  the  "null"  between  molecular  absorption  and  absorp- 
tion by  Inverse  Bremsstrahlung. 

The  two- port  propulsion  system  (III)  involves  Type  I and  Type  II 
phenomenon  in  addition  to  the  aerodynamic  window  flow  region  and  the 


nozzle  expansion  flow  that  are  respectively  fore  and  a ft  of  the  absorption 
chamber.  An  experiment  of  Type  III  involves  most  of  the  elements  of  a 
real  rocket  system  and  is  therefore  the  most  complex  type  of  experiment. 

The  selection  of  a Type  III  experiment  obviously  depends  upon  the  information 
of  Type  I and  Type  II  measurements. 

A number  of  Type  I absorption  experiments  are  described  in  the 
next  section.  This  discussion  is  followed  by  the  description  of  a high- 
pressure  burner  experiment  that  consists  of  perturbing  and  modifying  an 
H^O  flame  produced  by  the  stochiometric  mixing  of  gaseous  with  gaseous 
O^.  This  is  essentially  a Type  II  experiment  with  Type  I and  Type  III  cases 
also  possible;  the  limitation  of  the  burner  experiment  is  the  restriction  lo 
H^O  or  seeded  H^O  as  the  propellant.  The  last  section  outlines  what  is 
primarily  a Type  III  experiment:  low-pressure  two- port  simulation, 
a.  Absorption  Measurements 

The  design  of  the  cw  laser  propulsion  engine  depends  critically  upon 
the  absorptive  characteristics  of  candidate  propellants.  In  this  section  we 
describe  shock-tube  experiments  that  can  be  utilized  to  make  absorption 
measurements.  For  this  experiment  one  is  singularly  interested  in  laser 
beam  transmission  through  the  propellant  medium  which  will  typically  be  at 
high  temperature  and  high  pressure.  The  shock  tube  is  ideally  suited  for 
the  investigation  because  (1)  the  thermodynamic  conditions  (T,p)  can  be 
controlled  and  (2)  wide  range  of  propellants  (fuel  and  seed)  can  be  examined. 
Some  initial  assessments  have  been  made  of  utilizing  shock  tube  to  measure 
the  absorption  characteristics  of  H^O  (and  seeded  H^O). 

The  shock  tube  configuration  limits  the  temperature  range  of  H.,0 
to  approximately  3500°K  in  the  pressure  range  of  2 - 10  atm.  This 
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operating  range  is  shown  on  Figure  63  where  the  rurves  at  3,s,  and  8 torr 
correspond  to  the  initial  pressure  of  II  ,0  at  2'*8°K.  Specific  calculations  for 
the  3 torr  case  are  shown  on  Table  6. 

The  proposed  shock  tube  configuration  and  diagnostics  are  shown  on 
Figure  64.  Typical  shock  tube  conditions  of  330°K,  4 atm  1 1 ,0  can  be  ob- 
tained with  a 30  atm  11  , driver  pressure  and  saturated  water  vapor  t- S1  Cl 
providing  a driven  II  ,0  pressure  of  I torr.  Higher  temperature  and  pres- 
sures can  be  achieved  using  an  inert  gas  filler,  e.g.,  Argon. 

The  diagnostics  for  the  experiment  are  shown  on  Figure  t>4.  A 
Schlieren  setup  would  be  utilized  to  examine  the  boundary  layer  and  observe 
the  position  of  the  shock  wave.  The  boundary  layer  characteristics  of  in- 
terest are  the  thickness  and  the  thermal  characteristics.  The  primary 
diagnostics  would  be  10. 6 pm  transmission  measurements  employing  a 
CO,  laser.  l’he  windows  would  be  made  of  / nSe  (Irtran  4).  Temperature 
probes  along  the  tube  and  at  the  end  of  the  tube  would  be  employed  to  deter- 
mine the  shock  velocity.  Itoth  pure  II  ,0  and  seeded  11  ,0  would  be  examined. 
The  physical  shock  tube  would  be  approximately  8 feet  long  and  I 1 / 2 in. 
in  diameter.  The  windows  would  be  one  centimeter  from  the  end  wall. 

Note  that  the  possibility  of  absorption  of  the  laser  radiation  in  the 
shock  tube  boundary  layer  was  investigated  and  found  to  be  insignificant. 

b-  U'Btlzl  'reasure  I'nrner  experiment 

A potentially  attractive  means  of  laser  propulsion  is  by  means  id'  a 
two-port  rocket  configuration  utilizing  a CW  laser.  However,  scaling 
studies  have  shown  that  an  experimental  two- port  system  would  be  extremely 
small  unless  very  high  power  lasers  were  used.  In  order  to  avoid  the  prob- 
lems associated  with  such  small  systems,  a CW  simulation  experiment  was 
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Subscripts  refer  to  typical  shock  tube  configuration 
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Figure  64  Shock  Tube  Configuration  and  Diagnostics  (Illustrated  for  H^O) 
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conceived.  In  principle  this  system  uses  the  combustion  energy  of  and 
O-,  to  heat  the  resultant  H^O  to  temperatures  high  enough  to  be  close  to  the 
laser-heated  rocket  range.  Consequently  the  optical  properties  of  the  heated 
gas  can  be  investigated  without  using  a very  hig  power  laser.  We  are  pri- 
marily interested  in  determining  the  absorption  characteristics  at  high 
temperatures,  the  interaction  between  la ser- absorption  waves  and  the  flow 
field,  the  operation  of  an  aerodynamic  window  for  the  laser  absorption- flow 
problem;  and  the  coupled  radiation- gas  dynamics  phenomenon  in  two- port 
rockets . 

In  order  to  obtain  enthalpies  representative  of  laser  propulsion  ap- 
plications, we  will  use  the  High  Enthalpy  Absorption  Test  - Cell  (HEAT- 
Cell ] which  is  an  experiment  that  utilizes  the  laser  to  perturb  an  existing 
high  temperature  environment  to  a larger  enthalpy  environment.  Since  one 
of  our  candidate  propellants  is  H^O,  it  is  appropriate  to  use  a gll^/gO, 
motor  since  the  combustion  products  are  predominantly  H,Ot  ^ ' see  Figure 
65.  The  mixture  of  gaseous  H^  and  gaseous  O,  at  super  atmospheric  pres- 
sures provides  combustion  at  high  chamber  temperatures.  The  adiabatic- 

flame  (chamber)  temperature  for  a pressure  of  34.01  atm  is  shown  in 

(9) 

Figure  66  vs  the  O^/U^  mixture  ratio  (by  weight).'  The  stochiometnc 
mixture  ratio  of  8 gives  the  highest  temperature  of  approximately  3500°K. 

In  this  preliminary  analysis  we  will  consider  the  burner  flame  to  be  a homog- 
eneous one  from  the  temperature  point  of  view.  The  actual  structure  of  the 
flame  fronts  will  be  studied  subsequently.  We’will  consider  this  state  to 

(9)  Wilkins,  R.  L.  , "Theoretical  Evaluation  of  Chemical  Propellants,  " 
Prentice- Hall,  Inc.,  Englewood  Cliffs,  N . J . (1963)  p.  181. 
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EQUILIBRIUM  (2H2  -1-02) 


Figure  65  Mole  Fraction  of  Primary  Constituents  in  High 
Temperature  Water 
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be  the  initial  condition  in  our  test  cell.  This  high  temperature  - high 
pressure  environment  is  a severe  one  from  the  materials  point  of  view.  A 
small  160  lb  thrust  W^/O^  rocket  described  by  Ferriso*10*  et  al.  , utilizes 
water  cooling  to  preserve  the  rocket  chamber  integrity.  This  rocket  is 
shown  in  Figure  67.  A variation  of  this  rocket  can  be  used  for  the  investiga- 
tions discussed  above.  The  modifications  required  include: 

(i)  Buffer  gas  flow  to  insulate  the  test  cell  walls. 

(ii)  Laser  aerodynamic  window  designed  to  match  burner  flow. 

(iii)  Seedant  introduction  to  enhance  absorption. 

(iv)  Probe  ports. 

The  advantages  of  utilizing  a rocket  motor  are  indicated  in  Figure  68. 

First  the  combustion  process  provides  a high  enthalpy  initial  condition.  We 
can  see  this  on  the  skeleton  of  an  H^O  Mollier  Chart.  The  basic  combus- 
tion state  is  isolated  at  10  atm  and  3500°K — thereby  providing  an  increase 

4 . 

in  enthalpy  over  room  temperature  HzO  of  about  1.6  x 10  j/gm.  Absorp- 
tion characteristics  of  steam  at  these  temperatures  are  of  interest  by  them- 
selves. However,  we  would  like  to  obtain  higher  temperatures  to  simulate 
realistic  laser  propulsion  missions.  Hence,  to  obtain  an  interesting  temp- 
erature  of  4500°K,  the  laser  absorption  of  only  2.7  x 10  j/gm  would  be 
needed  as  shown  on  Figure  68.  The  Ferriso  rocket  mass  flow  of  0.  2 
kg/sec  implies  the  use  of  a 5.4  MW  laser  power  source  to  achieve  the 
indicated  Ah.  This  total  Ah  change  from  room  temperature  conditions 
would  yield  a total  Ah  - 4.  3 x 104  j/gm.  It  corresponds  to  an  ideal  specific 

JjfO)  Ferriso,  C.,  et  al.  , "High  Temperature  Infrared  Emission  and 
Absorption  Studies,"  General  Dynamics/Astronautics.  Scientific 
Report  AE61-0910,  (1961). 
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Figure  67  160#  Thrust  Rocket  Burner 
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Figure  68  Laser  Perturbed  State  on  H-O  Mollier  Chart 


impulse  of  about  920  seconds,  (see  Figure  41)  and  if  realistic  nozzle  ex- 

2 

pansion  pressure  ratios  (10  ) are  considered  one  would  still  obtain  a useful 
Ig  of  something  like  550  seconds.  The  perturbation  to  4500°K  at  constant 
pressure  is  an  interesting  one,  but  it  requires  a high  power  laser  to  achieve 
the  high  temperature  state.  The  skeletal  Mollier  chart  (Figure  68)  leads 
one  to  conclude  that  high  pressure  operation  is  favorable  for  laser  propul- 
sion applications  because  the  enthalpy  change  Ah  required  to  increase  the 
chamber  temperature  from  3500°K  to  4500°K  decreases  as  the  pressure 
inc  reases . 

The  various  HEAT- Cell  simulation  possibilities  are  outlined  in 
Table  7.  These  possible  experiments  are  illustrated  in  Figure  69. 

Experiment  1 is  a combustion  only  experiment.  The  main  purpose 
of  this  experimental  configuration  is  to  measure'  the'  absorption  properties 
of  water  at  high  temperatures  and  pressures.  Pressure  levels  up  to  50 
atm  would  be  investigated  and  peak  temperatures  around  "t500°K  would  be 
achieved.  A low  power  (10  W)  probe  laser  transverse'  te>  the-  burner  flame 
would  be  utilized  to  measure  the  absorptance  of  the  medium.  The  combus- 
tion only  experiment  would  thus  provide  = a ^ (p,  1 ) for  H^O,  In  addition 
to  H^O,  one  could  study  the  properties  of  seeded  H^O  by  introducing  a separ- 
ate inlet  port  or  by  adding  the*  se'edant  to  one1  of  the-  burner  propellant  s . 

Water  cooling  was  quite  satisfactory  for  the  General  Dynamics  rocke-t  motor 
discussed  above.  Water  cooling  could  be  used  and  the  increase  in  tempera- 
ture  of  the  water  could  be  measured  to  gain  some  understanding  of  wall  he'at 
transfer. 

Experiment  2 would  add  CW  laser  radiation  to  the  "combustion- 
generated inferno"  at  a relatively  low  power  level  (e.g.,  10  KW  from  the 

91 
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EXPERIMENT  1 

COMBUSTION  ONLY 


EXPERIMENT  2 


COMBUSTION  AND  LASER 


Figure  69  Heat-Cell  Experimental  Options 
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EXPERIMENT  3 

COMBUSTION  AND  LASER 


EXPERIMENT  4 


LASER  ONLY 


Figure  69  Continued 
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AERL,  Inc.  I1PL).  This  perturbation  of  the  rocket  scheme  would  provide 
valuable  information  beyond  that  pained  in  Experiment  1.  First,  the  absorp- 
tion of  laser  radiation  increases  the  chamber  temperature  level  so  the  probe 
laser  could  extend  the  range  of  absorption  coefficient  information.  Second, 
the  perturbing  beam  could  be  sent  in  parallel  to  the  burner  flow  as  shown 
in  Figure  69  or  it  could  be  sent  in  through  the  plume  (anti-parallel).  Axial 
absorption  would  be  monitored.  The  transition  from  molecular  to  inverse 
Bremsstrahlung  absorption  would  be  investigated.  The  stability  of  laser 
heating  problem  would  also  be  assessed.  The  heat  transfer  problem  would 
be  more  acute  here;  however,  water  cooling  would  still  maintain  the  integrity 
of  the  rocket  motor. 

Experiment  3 differs  from  Experiment  2 in  that  the  increase  in  per- 
turbing laser  power  necessitates  the  use  of  an  aerodynamic  window  to  get 
the  beam  into  the  chamber.  Employing  a laser  and  an  aerodynamic  window 
in  concert  with  a rocket  motor  situation  requires  a clever  design  for  the  in- 
let region.  The  solution  of  this  problem  then  provides  a more  interesting 
experiment  than  that  of  Experiment  2.  The  absorption  coefficient  can  be 
mapped  out  further  and  the  mating  of  the  aerowindow  with  "absorption" 
chamber  would  be  tested.  The  increase  in  laser  power  would  further  in- 
crease the  heat  transfer  to  the  rocket  chamber  walls.  The  window  flow  of 
hot  steam  as  well  as  the  rocket  flame  would  to  some  extent  protect  the  walls 
from  the  laser-absorbed  enthalpies. 

Experiment  4 begins  to  resemble  the  prototype.  The  rocket  motor 
aspect  of  the  HEAT-Cell  is  now  completely  eliminated.  All  the  enthalpy 
needed  to  raise  the  high  pressure  water  to  4500°K  now  comes  from  the  laser 
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(which  is  not  a perturbing  laser).  Further  absorption  coefficient  measure- 
ments can  thus  be  made,  and  transition  from  molecular  to  Inverse  Brems- 
strahlung  can  be  studied.  Propulsive  thrust  and  efficiency  can  be  measured 
and  the  viability  of  two- port  laser  propulsion  can  be  assessed. 

Note  that  the  single  port  concept  can  be  investigated  in  Experiment  3 
and  Experiment  4 by  turning  the  aerodynamic  window  flow  off  and  only  flow- 
ing propellant  into  the  rocket  chamber.  It  is  rather  obvious  that  the  series 
of  experiments  described  are  very  flexible;  additional  experiments  could  be 
mated  to  the  HEAT-Cell. 

c.  Low-Pressure  Two-Port  Simulation 

The  advantages  of  operating  a CW  engine  at  low  (absorption  chamber) 
pressures  have  been  mentioned  previously.  In  Figure  70  we  illustrate  a 
scaled  drawing  of  a CW  engine  operating  at  10  ^ atm  pressure.  The  charac- 
teristics of  the  engine  were  determined  from  the  scaling  maps  of  Section 
V.4.  One  of  our  purposes  in  designing  the  engine  was  to  make  it  as  large 
as  possible  in  all  its  elements.  For  a specific  impulse  of  1000  sec.  Figure  51 

provides  the  throat  area  when  a laser  power  level  is  chosen.  If  we  select 

2 

10  KW  as  the  laser  power,  then  A^  varies  between  8.3  and  14.9  cm  for  y 
between  y = 1.1  and  1.67  respectively.  Aerodynamic  window  calculal ions 
(in  particular  Figure  59)  indicate  that  one  obtains  the  largest  sizes  at  low 
and  y = 1-67.  The  choice  of  = 2 selects  case  3 (Table  5)  which  is 
illustrated  in  Figure  62.  The  length  (14.6  cm)  of  the  subsonic  diffuser  is 
obtained  by  making  the  diffuser  diverge  at  7°. 

The  length  of  the  absorption  chamber  has  somewhat  arbitrarily  been 
selected  to  be  equal  to  its  diameter.  Rockets  motors  typically  have  square 
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dimensions.  The  length  of  this  chamber  is  the  weakest  part  of  the  design; 
more  specifically,  since  absorption  at  low  pressures  may  require  large  ab- 
sorption lengths,  the  absorption  chamber  may  have  to  be  "stretched"  out. 
The  use  of  a short-  abso rption-  length  seedant  would  keep  the  chamber  geom- 
etry square.  An  area  ratio  between  the  throat  (A^)  and  the  nozzle  exit  was 
chosen  to  be  1 00 . This  is  a state-of-the-art  number.  Foray  =1.67  gas, 
an  expansion  pressure  ratio  (p5/P7)  of  1 . 5 x 1 0 is  obtained  (see  Figure  38). 
The  length  of  the  nozzle  then  becomes  73.2  cm  and  as  shown  on  the  scaled 
drawing,  the  nozzle  represents  the  major  portion  of  the  engine  length 
(1 . 24  m) . The  system  mass  flow  is  0.2  gm/ sec . 

It  now  becomes  evident  that  the  low-pressure  engine  is  well  suited 
to  laboratory  experiments.  Most  significantly,  at  low  pressures  available 
high-power  lasers,  e.g.,  AERL,  Inc's  HPL  would  provide  the  requisite 
laser  radiation.  The  mass  flows  are  not  excessive  and  the  heat  transfer 
rates  will  be  significantly  reduced  from  the  high-pressure  simulation.  The 
experimentation  would  focus  upon  the  use  of  a seeded-monatomic- like  gas 
(y  = 1.67).  Since  the  operation  is  at  sub-atmospheric  pressure  levels,  the 
experiment  would  be  performed  in  an  evacuated  environment.  The  ambient 

6 - 3 

pressure  level  would  have  to  be  less  th an  6 x 10  atm  (5.1  x 10  torr). 
The  primary  diagnostic  for  the  experiment  would  be  the  measurement  of 
thrust  — and  hence  specific  impulse. 

The  goal  of  the  experiment  would  be  "proof  of  principle".  Scaling 
to  high  chamber  pressures  might,  however,  be  difficult  due  to  different 
absorption  behavior  even  if  the  low-pressure  simulation  were  successful. 
We  would  then  view  the  low-pressure  engine  as  useful  for  exoatmospheric 
applications . 


98 


I 


VI.  BEAM  TRANSMISSION 


1.  INTRODUCTION 

We  wish  to  find  the  conditions  under  which  a beam  of  power  in  the 
9 

order  of  1 0 W can  be  transmitted  through  the  atmosphere.  We  will  con- 
fine the  discussion  to  a CO^  laser  beam  at  10.6  pm,  making  an  angle  with 
the  vertical  of  60°  or  less. 

The  only  significant  absorbers  present  in  the  atmosphere  are  CO^ 
and  H^O . The  first  of  these  absorbs  by  the  inverse  of  the  CO^  laser  tran- 
sition, CO2  (100)  + he  -*  C02  (001).  This  transition  is  capable  of  being 
saturated  or  bleached,  and  when  the  molecule  is  bleached  absorption  be- 
comes effectively  negligible.  The  potential  for  bleaching  out  the  CO^  ab- 
sorption component^'  • ^1  extremely  important  for  propagation  of  pulsed 
CO^  laser  beams,  especially  those  directed  near  to  the  vertical,  as  in  the 
case  of  laser  propulsion  applications.  We  will  derive  the  conditions  under 
which  this  effect  is  significant  below.  We  note  that  CO^  is  present  in  the 
atmosphere  at  a molar  concentration  of  330  x 10  ^ . The  second  absorber 

is  H^O;  there  is  some  indication  that  absorption  takes  place  in  the  contin- 

113) 

uum  of  a weakly  bound  H^O  dimer.  ' ' The  extremely  rapid  vibrational 


(11)  Wood,  A.D.,  Camac,  M.,  and  Gerry,  A.E.,  Appl.  Optics  10,  1877 
(1971). 

(12)  Avizonis,  P.  V . , Butts,  R.,  and  Hegge,  B.  , Appl.  Optics  14,  1911 
(1975). 

(13)  Adiks,  T.  , Arefyev,  V.,  and  Dianovklokov,  V.,  Sov.  JQE  5,  481  (1975). 
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readjustments  of  the  11,0  molecule  imply  that  the  molecule  cannot  be  bleached, 
anti  the  absorption  it  produces  must  be  added  to  that  of  the  CO,  molecule.  The 
11,0  content  of  the  atmosphere  is  extremely  variable,  and  in  order  to  estimate 
transmission  we  have  adopted  the  values  for  the  U.S.  standard  atmosphere 
given  by  Valley,'  which  correspond  to  1%  molar  H^O  at  sea  level. 

The  absorption  coefficients  obtained  will  then  be  used  to  calculate 
the  degree  of  distortion  (thermal  blooming)  produced  in  the  beam;  the  phase 
corrections  required  to  reduce  this  phenomenon  and  the  maximum  correction 
obtainable  will  then  be  computed.  These  two  calculations  have  been  performed 

(15) 

on  the  Lincoln  Laboratory  propagation  code,  kindly  arranged  by  J.  Herrmann. 
We  proceed  to  a more  detailed  discussion. 

CO,  ADSORPTION  BLEACHING 

The  CO,  molecule  absorbs  10.6  pm  radiation  on  the  resonant 
(100)  -*•  (001)  transition.  The  molecular  vibrational  energy  level  diagram 
is  shown  in  Figure  71,  in  which  the  major  energy  transfer  processes  are 
indicated.  The  photon  absorption  transfers  the  molecule  to  the  (001)  state. 

There  are  four  possible  loss  modes  from  this  state;  transfer  to  a nitrogen 
molecule  in  the  reaction 

CO2(001)  + (0)  - CO.,  (000)  + N , ( 1 ) ; (166) 

collisional  deactivation  to  the  bending  mode  in  the  reaction 

C02  (001)  + M -*  CO2(110)  + M,  (167) 


i 


(14)  Valley,  S.  L.  , Atmospheric  Temperature,  Density,  Pressure  and 

Moisture.  Handbook  of  Geophysics  and  Space  Environments,  AFCRL, 
Office  of  Aerospace  Res.,  U.S.A.F.  1965. 

(I1')  Hermann,  J.,  private  communication. 
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ASSEMPTIONS 


1.  SYMMETRIC  STRETCH  & BENDING  IN  EQUILIBRIUM 

2.  DURING  PULSE  t » b 


3. 

001  LOSS: 

VV  TRANSFER  TO  Nj 

< 

< 

VV  TRANSFER  TO  030 

rdi 
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010  LOSS: 

VT  TRANSFER  TO  KINETIC 

Td2 

5. 

ATMOSPHERIC  CO  CONTENT:  330  ppm  MOLAR 

6. 

ONLY  C02 

ABSORPTION  IS  CONSIDERED 

7.  ONLY  N2.  02  AND  H20  EFFECT  ON  Td1  AND  T d2  CONSIDERED 
H 4256 


Figure  71  CO^  " Kinetic  Model 
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where  M is  N^,  02  or  H20;  stimulated  emission  by  another  photon  in  the 
laser  transition,  the  reverse  of  which  is  the  absorption  process 

CO,  (100)  + hv  - CO2(001);  O68' 

or  direct  radiative  transition  to  the  ground  state.  The  latter  process  is  very 
much  slower  than  the  others  and  may  be  neglected.  The  most  rapid  process 
is  transfer  to  nitrogen;  for  the  conditions  here  the  nitrogen  acts  as  a vast 
thermal  bath  into  which  the  absorbed  energy  is  deposited.  The  N->  vibra- 
tional energy  either  deactivates  via  collisions  with  H20  molecules  (Tq  ~ 30  ms) 
or  by  transferring  the  energy  back  to  a C02  molecule,  which  then  collisionally 
deactivates  to  the  bending  mode,  and  thence  to  the  ground  state  (Tc  ~ 10  ms). 

The  characteristic  times  above  are  for  sea  level,  and  the  H20  deactivation 
decreases  much  more  rapidly  with  height  than  that  due  to  C02  transfer  and 
deactivation;  we  therefore  neglect  deactivation  of  N2  by  H2G.  The  character- 
istic times  may  be  obtained  from  the  rate  constants , ( l6)  and  we  find  that  proc- 
ess  (366)  is  ~ 100  times  faster  than  process  (367),  and  thus  dominates  the  loss 
processes  from  CC>2  (001).  The  time  for  process  (366)  is  labelled  Tyv,  and 

that  for  (367)  is  t these  are  shown  as  functions  of  altitude  in  Figure  72.  i j 

Simultaneously,  depletion  of  the  (100)  state  by  photon  absorption  re-  r 

suits  in  its  repopulation  via  collisional  transfer  from  the  bending  mode,  which  ^ j 

is  extremely  rapid,  and  this  leads  to  repopulation  of  the  bending  mode  from  j 

I j 

the  ground  state.  This  process  is  strongly  catalyzed  by  H20,  and  the  char-  ; j 

acteristic  time  for  (010)  - (000)  transfer,  labelled  Td2  in  Figure  71  and 
shown  as  function  of  altitude  in  Figure  72,  is  dependent  on  H20  concentration. 

Thus  the  fourth  process  important  for  bleaching  is: 

CO,  (000)  + M -*  CO?  (010)  + M.  (369) 

(16)  Lewis,  P.  F.  and  Trainor,  D.W.,  AERL  AMP  422  (Nov.  1974). 
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We  now  proceed  to  write  the  kinetic  equations  for  these  processes,  first 
listing  the  assumptions  made  in  the  derivation, 
a . Assumptions 

1.  Atmospheric  constituents  are  O^.  CO^  and  H^O.  Collisional 

deactivation  by  process  (367)  where  or  are  the  third  body  have  the 

identical  rate  = 3x10  ^ cm^/sec  for  T ~ 300°K,  and  for  H^O 

k , . = 4 x 10"13.  We  ignore  the  small  variation  of  k with  temperature, 
d 1 a 1 

Collisional  VT  transfer  by  process  (369)  with  or  O.,  as  third  body  have 
the  identical  reaction  coefficient  k^  = 3 x 10  ^ e ^60/Tg  cm^y'seC(  while 
for  H00  k,-,  = 2 x 10-11  e‘9()0/Tg>  where  Tg  is  gas  kinetic  temperature. 

C02  - N2  VV  transfer  (process  (366))  has  the  rate  constant  kvv  = 5.8  x 10 
cm3/ sec.  All  the  above  rates  are  obtained  from  Lewis  and  Trainor.*1^ 

2.  The  atmospheric  temperature,  density  and  humidity  for  the  stand- 

(14) 

ard  U.S.  atmosphere  given  by  Valley'  will  be  used. 

3.  The  CO^  vibrational  modes  are  harmonic,  and  while  Vj  and 
are  connected  by  the  rapid  exchange  of  vibrational  information  between  the 
two  modes,  the  series  is  connected  with  the  other  modes  only  by  proc- 
esses (367)  and  (368). 

4.  For  repeated  laser  pulses,  typically  250  pps  of  ~ 10  ps  pulse 
duration,  the  CO^  (001)  population  fraction  during  the  pulse  will  always  be 
much  greater  than  the  (1)  population  (N.,  bath  assumption),  so  that  for 
the  present  discussion  the  reverse  rate  of  process  (366)  will  be  negligibly 
slow  compared  to  the  forward  rate. 

5.  Gas  temperature  change  will  affect  the  molecular  kinetics  insig- 
nificantly, and  constant  temperature  may  be  assumed  during  the  absorption 
process . 
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b.  Kinetic  Equations 

The  rate  equation  for  the  CO,  (001)  population  may  be  written: 
£[CO2(001)]  = - kvv{[N2  (0)][CO2  (001)]  - [n2  ( 1 )]  [CO,  (000)]( 

- [0O2  (001)1/Tdl  + ft  ■ 

where  molar  concentrations  per  cubic  centimeter  are  shown  in  square 

brackets,  ^ is  flux  density  at  10.6  pm,  he  is  photon  energy  at  10.6  pm, 

and  q is  the  absorption  coefficient.  The  above  equation  may  be  simplified 

by  defining  r = [CO,  (001)]/ [CO,  (000)]  and  b = [n  (1)]/[N  (0)1  . Then 

-3380/T3  -3395/Tn 

r = e where  T(  is  the  v,  vibrational  temperature , and  b = e , 

where  T^,  is  N,  vibrational  temperature.  Due  to  the  photon  pumping 

T,  » T^j,  and  hence  r ^ b;  we  neglect  b.  Using  the  definition  of  Tyy  we 


ir  f 1 1 ! 

* = ' r[Tvv  + Tdl  J 


he  4 Nt 


(370) 


where  (I1  is  CO,  molar  fraction  and  N_  is  total  gas  number  density.  Ground 

X 

state  depletion  may  be  neglected  since  r « 1.  At  this  point  we  note  that  a is 


given  by 


where 


a = F (s  - r) 


-J  (J+  Df)r/T& 


/ Mc  > Q er  »c  Nt  (2J  + 1)  e * 

Zv  k n0  4 TZ  1 T 3//Z  N.  o~(l  + M /M.)1/2 

R > . J J c J 


(372) 


(17)  Douglas -Hamilton,  D.  H.  and  Lowder,  R.S.,  AERL  CO,  Laser  Kinetics 
Handbook,  1974. 


r 


. 


for  small  r and  s,  and  the  quantities  have  the  following  meaning  and  values 


for  P (20): 


CO^  mol  wt 
mol  wt 

Boltzman  const. 


Avogadro  const. 

Wavelength 

Upperstate  Rotational 
Quantum  Number 

Rotational  Spacing 

Gas  temperature 

Number  density  of 
species  j (N^  + O^) 

Collisional  Cross 
Section  for  j (N^  + O^) 


1 . 38054  x 10 


6. 023  x 10' 


10.59  jim 
19  for  P (20) 


0.  53°K 


Nj  ~ ^N  NT  + ^o  NT  = NT 


- 1 5 2 

8.7x10  13  cm 


Inverse  Einstein  transition  5.38  sec 
probability 

2 

and  s = CO?  (010)  / CO-,  (000)  , so  s is  the  lower  state  population  frac- 


Substituting  the  appropriate  values  for  air  in  Eq.  (372)  yields  the 
numerical  value  for  the  P (20)  transition: 


7.39  -201''Tf 


(373) 


For  constant  s,  the  characteristic  time  for  r to  reach  the  steady  state  value 

r , from  Eq.  (370),  will  be: 
s s 

tc  = tvv  r* + ^/^vi  ' 
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whore 


and 


VV 


= [rv'v  + Tdl] 


- 1 


<^V  = hv  <i>c  NT/(FTVV) 

The  term  if,  is  then  a saturation  flux  for  vibrational  transfer  to  N^,  and 
defines  the  flux  at  which  photon  pumping  into  CO,  (001)  becomes  as  rapid 
as  vibrational  transfer  out  of  C.O^,  (001).  The  numerical  value  of  if,  ^ in  the 
atmosphere  is  shown  in  Figure  7 1.  Numerical  values  for  the  characteristic 
time  are  shown  in  Figure  74  as  function  of  altitude,  and  we  see  that  for  fluxes 
greater  than  10^  W/cm^  the  characteristic  time  is  always  less  than  I p sec. 

We  also  see  from  Eq.  ( 170)  that,  using  the  definition  of  </)  v , the 
s teady-state  value  of  (001 ) population  fraction  will  be 

2 

s 

s s l + </>  y /<f> 

so  that  the  absorption  coefficient  when  the  upper  state  has  reached  steady 
state  will  be 


( 174) 


1 + (ft  / if)  v 

_2  . _ -"«°/T 

where  u - F ( r - s ) is  the  unperturbed  value  of  n , with  r e 
-960/T 

and  s = e 

We  conclude  that  the  value  of  n given  by  Fq.  (17  5)  will  be  reached 
in  the  times  given  in  Figure  74,  i.e.  , very  rapidly  compared  with  pulse 
durations  of  10-20  psec.  We  also  note  that  depletion  of  CO,  (100)  will  re- 
sult in  a further  decrease  in  the  absorption  coefficient  on  a longer  time 
scale.  The  kinetic  equation  for  CO ^ (100),  allowing  for  the  double 


' i 


L. 


107 


LOCAL  VALUE  OF 
SATURATION  FLUX,  W CM 


r 


Figure  73  CO^  kinetic  Saturation  Fluxes  vs  Altitude 

CO^  (001)  -»  (1)  Corresponds  to<f>^  in  Test,  and 

CO.  (000)  -*  (100)  Corresponds  to  A 
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degeneracy  of  the  bending  mode  and  assuming  identical  vibrational  tempera  - 

-960/ T 

tures  for  v j and  v^,  again  ignoring  ground-state  depletion  since  e & 

« 1,  becomes: 


? £*£  _ 2 a 

dt  ~ ‘ hv  <|'c  Nt 


2 (s  - s) 


(376) 


We  observe  that  since  Tyy  « (see  Figure  72),  the  CO^  (001)  steady 

state  given  by  Eq.  (374)  will  be  reached  very  rapidly  compared  to  the  time 

taken  for  CO-,  (100)  to  reach  steady  state.  We  are  therefore  justified  in 

taking  r = r , substituting  in  Eq.  (376)  and  integrating  the  equation  to  ob- 
s s 

tain  CO 2 (010)  and  CO^  (100)  population  fractions  as  function  of  time.  De- 
fining (j>  = hv  N^/(  F Tj^) , which  is  the  flux  for  which  CO^  (100)  depop- 

ulation becomes  significant,  and  defining  the  following  quantities: 

V*.  - 3T^/Tdi 


A = 11  + 4 A s 


s = e 


-960/ Tg 


2 A s 4-  1 - A 

2 A s + 1 + A 


the  solution  of  Eq.  ( 376)  becomes: 


T + T75T 


.14-1). 

. . z -2 

4 A s 


w fjf  ae 

i _ct 

1 - ae 


( 377) 
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Here  again  the  important  point  is  that  the  steady-state  value  (for  t - ® ) is 


a 

nr 


4A  s 


(A  - ir 

(i  + <i> /</>)-  ’ 


which  is  reached  in  a characteristic  time 


T = T.,/A 

c d2' 


(378) 


(379) 


We  conclude  therefore  that  the  absorption  coefficient  will  rapidly 
drop,  in  times  given  in  Figure  74,  to  the  values  given  by  Eq.  (375).  It  will 
continue  to  drop,  in  times  given  by  Eq.  (379),  to  the  values  given  by  Eq 
(378).  We  therefore  plot  the  Upper  State  Bottleneck  (Eq.  (375)),  Lower 
State  Bottleneck  (Eq.  (378)),  and  characteristic  time  for  moving  from  the 
first  condition  to  the  second  (Eq.  (379))  in  Figure  75,  as  functions  of  alti- 
tude. We  notice  that  the  important  contribution  comes  from  the  upper  state 
bottleneck,  and  that  subsequent  lower-state  depletion  decreases  the  absorp- 
tion coefficient  by  factors  from  2 to  4.  At  3 km  altitude  the  transition  be- 
tween these  absorption  modes  takes  only  1 jj  sec;  but  at  10  km  it  takes  10  p sec , 
and  at  20  km  about  100  psec.  We  conclude  that  it  will  be  conservative  to 
assume  that  only  the  upper  state  bottleneck  will  be  important  for  the  laser 

propulsion  case.  We  will  proceed  therefore  to  calculate  the  absorption  in- 

/R 

a dz , ignoring  the  time  dependence  and  taking 
absorptance  as  given  by  Eq.  (375). 

c . Absorption  Integral  and  Beam  Distortion 

The  Distortion  Number  N,  which  is  a measure  of  the  beam  phase 
distortion  produced  by  thermal  blooming,  is  defined  as 

k 


IP  a dz 


J D(V  + R ci>) 


l l 1 


N 


( 380) 


NORMALIZED  ABSORPTION  COEFFICIENT,  a/a 


where  k is  the  wave  number,  IP  and  D the  beam  power  and  diameter,  x 
is  the  initial  height  and  R is  the  range,  V the  wind  velocity  transverse  to 
the  beam,  Roi)  is  the  averaged  slew  velocity  and  J is  constant.  We  note 
that  a enters  only  as  an  integral  over  the  range.  For  vertical  propagation, 
absorption  of  CO,  P (20)  radiation  is  unimportant  after  about  20  km,  and 
the  beam  will  not  become  further  distorted  thereafter.  We  will  therefore 
calculate  the  absorptance  integral  from  height  x to  20  km,  based  (as  dis- 
cussed above)  on  the  upper-state  bottleneck.  Using  the  standard  atmosphere 
values  we  obtain  the  absorbed  fraction  of  the  beam  shown  in  Figure  7b,  in 
which  the  effects  of  11,0  as  an  absorber  and  as  a vibrational  deactivant  have 
been  included.  For  comparison  the  bed  absorption  at  height  11  is  shown  in 
Figure  76a  as  function  of  flux  and  time  into  pulse.  Thus  starting  from  a 
height  of  3 km,  with  neglible  flux  about  4 3%  of  the  beam  is  absorbed  ami 
57%  transmitted  (this  ratio  is  less  favorable  in  more  humid  atmospheres). 

A flux  of  10^  W/cm*”,  however,  which  should  still  be  above  the  aerosol 
breakdown  limit,  would  only  be  12%  absorbed;  88%  transmitted.  We  would 
expect  the  distortion  number  to  be  a factor  of  4 lower  than  in  the  unbleached 
case.  Consequently  bleaching  has  a significant  effect  on  beam  transmission. 
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Figure  76a  Bleaching  Atmospheric  CO^ 

Local  absorption  coefficent  a normalized  to  local  unbleached 
value  a,  as  a function  of  time  r into  pulse  at  constant  flux 
level  </>  and  at  constant  altitude.  Altitude  H is  in  kilometers. 
This  absorption  coefficient  must  be  added  to  that  due  to  local 


water  vapor  content. 

(a)  <t>  = 105  W/cm2 
(c)  $>  = 3 x 10^  W/cm2 


(b)  <J>  = 106  W/cm2 
(d)  <t>  = 1 0^  W / cm2 


3.  PHASE  DISTORTION  AND  PHASE  CORRECTION 


We  are  indebted  to  J.  Herrmann  at  Lincoln  Laboratory^®*  for 
performing  calculations  on  the  Lincoln  Laboratory  atmospheric  phase  distor- 
tion code  in  which  the  amount  of  blooming  produced  was  estimated,  and  the 
phase  corr  ection  (and  the  degree  to  which  the  beam  can  be  corrected)  were 
calculated.  In  the  computations  so  far  completed  bleaching  was  not  taken 
into  account;  we  can  however  estimate  its  effect  by  using  the  corrected  value 

of  / adz  (see  Figure  76).  Blooming  and  the  requisite  phase  corrections 
x 

calculated  for  various  diameter  and  slew  rates  will  be  summarized  below. 

The  following  calculations  apply  to  cw  and  mp  lasers  pointing  upward 
through  the  atmosphere.  Scaling  laws  for  such  a complicated  scenario  are 
not  available;  therefore,  we  had  to  run  specific  cases  in  order  to  explore 
the  influence  of  at  least  some  of  the  relevant  parameters. 

All  propagation  calculations  were  performed  with  the  following 
parameters  fixed: 

Beam  shape  = truncated  Gaussian  with  l/e  power  diameter 
equal  to  2 m. 

Initial  altitude  of  laser  = 3 km, 

Wind  velocity  = 10  m/sec  of  uniform  direction  at  all  altitudes. 

Two  wavelengths  were  used:  10.6  n ( CC and  4.8  fj  (CO),  for  calcula- 
tions with  diffraction- limited  beams.  For  each  of  the  two  wavelengths,  an 
appropriate  altitude- dependent  absorption  coefficient  was  used,  correspond- 
ing to  mid- latitude  summer  conditions.  For  the  CO^  the  effect  of  kinetic 
cooling  is  included. 

(18)  J.  Herrmann,  D.P.  Greenwood,  P.  Kafalas.  Anticipated  Propagation 
Conditions  for  High  Energy  Lasers  at  White  Sands  Missile  Range,  Con- 
tract DAAH01-75-C-  1272. 

(19)  J.  Herrmann,  Private  Communication. 
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The  parameters  range,  slewing  rate,  elevation  angle,  and  the  over- 
lap number  (for  multi-pulse  mode)  were  partially  explored. 

The  results  are  presented  as  relative  area  increase  as  a function  of 
average  power,  where  the  area  of  a beam  is  defined  as  the  area  containing 
63.2%  of  the  total  power.  In  some  cases,  a diameter  of  an  aperture  re- 
ceiving this  fraction  of  power  is  used. 

Effects  due  to  aerosols  and  turbulence  are  not  included. 

For  some  cases,  a predictive  phase  correction  of  thermal  blooming 
was  performed. 

We  give  a summary  of  the  principal  results. 

a . Range  Dependence 

The  relative  area  increase  is  only  a weak  function  of  the  range,  as 
shown  in  Figure  77,  where  we  plot  the  area  increase  as  a function  of  power 
for  two  ranges,  for  cw  and  multipulse  (MP)  beams,  with  pulse  overlap  num- 
ber N ^ = 2.  (N^  is  number  of  laser  pulses  passing  through  given  volume  of 

gas  at  laser  exit.)  The  results  show  that  the  effective  atmospheric  bloom- 
ing lens  is  a thin  lens  located  near  the  transmitter.  We  performed  there- 
fore most  of  the  calculations  for  the  200-km  range. 

b.  Slewing  Rate  Dependence 

The  relative  area  increase  is  shown  in  Figure  78  for  two  slewing 
rates  for  a cw  and  MP  ( N ^ = 8)  beam.  Ihe  area  increase  is  roughly  pro- 
portional to  the  slewing  rate;  therefore,  we  limited  our  calculations  to  a 
slewing  rate  of  40  mrad/sec. 

c . Elevation  Angle  Dependence 

For  cw  beams  without  slewing,  the  effective  lens  at  an  elevation 
angle  of  30°  is  twice  the  effective  lens  at  90°.  llut  the  presence  of  slewing 


Figure  77 
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Effect  of  Range  on  Multipulse  CO.,  Thermal  Blooming 
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Figure  78 


Effect  of  Slew-Rate  for  CW  and  MP  CO,  Laser  Beams 
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reduces  this  effect,  as  shown  in  Figure  79  for  a 40-mrad/sec  slewing  rate. 
The  elevation  angle  dependence  is  smaller  for  the  multi-pulse  case  because 
the  overlap  distance  is  independent  of  the  elevation  angle. 


d . Kinetic  Cooling 

Only  for  the  cw-CO^  laser  do  kinetic  cooling  effects  appear.  We 
compared  the  same  conditions  with  and  without  kinetic  cooling  (by  reducing 
the  relaxation  time  to  zero).  The  pres  . -nee  of  kinetic  cooling  reduces  the 
thermal  blooming  by  a factor  of  two.  This  is  shown  in  Figure  80.  The  MP 
lasers  are  not  affected  by  kinetic  cooling  because  the  overlap  does  not  reach 
into  the  region  (6  - 10-km  altitude)  where  the  relaxation  time  becomes  ef- 
fective. The  effects  of  CO^  bleaching  have  not  been  taken  into  account  here. 

e . Overlap  Number 

The  MP  case  has  one  parameter  more  than  the  cw  case,  for  which 
we  use  the  overlap  number.  The  overlap  number  Nq  is  defined  as  the  num- 

< 

ber  of  pulses  per  flow  time  at  the  transmitter.  In  Figure  81  we  show  the 
area  increase  as  a function  of  average  power  for  three  overlap  numbers. 

t 

An  overlap  number  going  to  infinity  approaches  cw  conditions.  In  our  case, 
this  transition  is  complicated  due  to  the  appearance  of  kinetic  cooling  for 

f 

this  limit.  ; 

i 

f . Blooming  for  CO  j 

Thermal  blooming  for  CO  wavelengths  and  absorptions  is  shown  in 

I 

Figure  82  for  cw  and  two  MP  cases.  The  integrated  absorption  coefficient 
smaller  for  CO  than  for  CO^  (0.242  vs  0.541),  but  the  absorption  is  more 
concentrated  near  the  transmitter.  Slewing  will  therefore  reduce  CO^ 
blooming  more  than  CO  blooming.  The  reference  area  is  smaller  for  CO; 
therefore,  we  get  a larger  relative  blooming  effect  for  CO  than  for  CO^. 
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Figure  70  Elevation  Angle  Dependence  of  Thermal  Blooming  for 
CO.,;  Multipulse  and  CW 
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Thermal  Rlooming  Dependence  on  Overlap  Number  for 
CO2  Multipulse 
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g . Comparison  of  CO.,  a nd  CO 

In  order  to  compart*  tht*  thermal  blooming  for  CO^  and  CO,  we  plot 
in  Figure  83  the  diameter  (of  the  63.2%  flux  area)  in  absolute  units  as  a 
function  of  power.  The  difference  in  the  blooming  effects  is  then  properly 
compa  red . 

h . Phase  Correction 

We  can  perform  different  phase  corrections  far  thermal  blooming. 
Our  present  cw  code  can  be  used  with  the  simple  predictive  correction, 
which  does  not  permit  a proper  inclusion  of  kinetic  cooling.  We  performed 
therefore  the  cw  phase  correction  without  kinetic  cooling  effects.  Our  MP 
code  can  be  used  with  the  iterative,  predictive  phase  correction.  Adaptive 
corrections  were  not  made. 

i.  Search  for  the  Phase  Corrections  (CO.,  - cw) 

Phase  corrections  for  the  CO.,  - cw  laser  were  performed  for  two 
altitudes,  10  and  200  km,  and  they  are  presented  in  Figure  84,  where  we 
plot  the  relative  area  increase  as  a function  of  the  scale  parameter.  The 
correction  for  10  km  is  not  as  good  as  it  is  for  200  km,  because  the  effec- 
tive thermal- blooming  lens  is  less  thin  than  for  the  200  km  cast'  (but  it  is 
still  thin) . 

j.  Phase  Correction  for  CO„,  - cw 

In  Figure  85  we  plot  the  relative  area  increase  as  a function  of  power 
for  the  CO,  - cw  beam  for  two  altitudes,  uncorrected  and  corrected. 

k.  Comparison  of  Corrections 

The  correction  for  the  10  km  range  is  less  complete  than  for  the 


200  km  range,  but  the  reference  beam  is  smaller.  We  plot,  therefore,  the 
diameter  in  absolute  units  for  both  cases  in  Figure  86.  If  the  criterion 
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Figure  85  Phase  Correction  for  Thermal  Blooming  in  CW  CO, 
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of  usefulness  is  a given  target  size,  then  the  lesser  is  a given  target 
size,  then  the  lesser  correction  of  10  km  is  not  important. 

1 . Correction  for  MP 

A few  corrections  were  performed  for  the  CO,-MP  case,  shown  in 
Figure  87.  Multi-pulse  cases  can  be  corrected  to  a much  higher  power 
level  than  cw  cases.  The  correction  is  very  good  for  the  lower  overlap 
number.  The  phase  excursion  required  for  the  case  P = 300  (Nq  = 2)  is 
about  100  radians  on  the  wavefront,  corresponding  to  about  100  microns 
mirror  displacement. 

We  conclude  from  the  above  survey  that  for  low  pulse  numbers 
multipulse  is  preferable  to  cw  for  transmission  through  the  atmosphere. 

We  will  consider  only  CO^  10.6  pm  radiation,  since  CO  propagation  is  in- 
ferior, and  it  is  difficult  to  see  how  the  use  of  chemical  lasers  at  shorter 
wavelengths  would  not  be  considerably  more  expensive  than  CO,  lasers. 

Since  it  appears  that  a stacked  array  of  modular  units  of  the  CO^ 
laser  type  will  be  the  most  promising  method  of  beam  generation,  the  use 
of  a large  combining  chandelier  mirror  provides  a natural  way  to  impose 
phase  corrections,  since  the  mirror  could  be  faceted  with  each  facet  separ- 
ately controllable  and  each  laser  module  directing  its  radiation  at  the  cor- 
responding facet.  The  modules  would  be  phase- locked  by  a pulse  from  a 
low-power  pulse- initiating  laser,  so  that  phase  corrections  to  the  combined 
beam  could  be  provided  by  independently  adjusting  each  facet  on  the  chande- 
lier mirror.  An  alternate  scheme  would  be  to  provide  phase  modulation  of 
the  low  power  pulse  initiating  laser,  which  would  produce  the  required  phase 
correction  in  the  high  power  amplified  beam.  This  would  have  the  advantage 
of  manipulating  the  beam  while  at  low  power. 


130 


We  are  therefore  interested  in  the  propagation  and  phase  correction 
of  a large  pulsed  beam,  between  10  and  30  meters  in  diameter.  Accordingly 
calculations  have  been  performed  on  the  Lincoln  Laboratory  (Herrmann- 
Bradley)  code  and  the  results  are  shown  in  Figure  88,  in  which  beam  area 
at  1000  Km  is  plotted  against  laser  power,  for  various  mirror  diameters 
and  slew  rates.  At  these  sizes  the  overlap  number  is  typically  Nq  « 100-300, 


for  a 100  pps  laser,  and  the  resultant  beam  is  effectively  cw  from  a propaga- 
tion standpoint.  Nevertheless,  the  large  mirror  diameter  provides  a beam 
whose  distortions  can  be  corrected  so  well  that  at  1000  Km  range  (the  maxi- 
mum Brennschluss  distance  for  laser  propulsion  to  geosynchronous  orbit) 
the  energy  within  a 2 meter  spot  will  be  1 09  watt  provided  the  effective  slew 
rate  is  2.5  mrad/sec.  We  should  point  out  here  that  for  a launch  from,  for 
example,  Pike' s Peak,  the  average  winds  are  so  high  that  the  effective  slew 
rate  is  much  greater  than  this  value  throughout  the  absorbing  region  of  the 
atmosphere.  In  addition  the  effect  of  CC>2  bleaching  gives  a further  advan- 
tage by  reducing  the  distortion  number  by  a factor  of  4,  as  discussed  above. 

The  results  have  been  summarized  by  scaling  the  distortion  number 
from  Eq.  (380)  as  function  of  bleaching  (dependent  via  / a dx  on  the  laser 
flux  density)  and  mirror  diameter.  The  distortion  number  may  be  written 


where  y is  the  adiabatic  exponent  for  air  and  Cs  is  the  sound  speed,  IP  is 
the  laser  power,  V the  wind  velocity  (including  slew  rate,  D the  mirror 
diameter  and  X the  optical  wavelength.  Thus  for  constant  power  and 
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Figure  88  Phase  Corrected  Area  of  Laser  Beam  at  Target  1000  km 
Distant,  vs  Laser  Power  for  Various  Mirror  Diameter 
and  Slew  Rates.  At  these  diameters  CW  and  MP  lasers 
propagation  is  indistinguishable. 


wavelength  the  distortion  number  scales  as 

r R 

N CC  J a dz/D.  (381) 

X 

The  abscissa  of  Eq.  (381)  scales^®'  as  and  since  the  effect  of  bleach- 

ing is  to  reduce  adz  by  a factor  of  two,  it  is  evident  that  when  bleaching 

3/2 

is  considered  the  beam  power  at  target  will  be  2 ' that  when  bleaching  is 

ignored.  For  the  case  of  the  20  meter  diameter  mirror  we  therefore  con- 

9 

elude  that  a phase-corrected  10  watt  beam  can  be  propagated  into  a 2 meter 
spot  at  1000  km.  The  propagation  of  an  adequately  tight  beam  for  laser 
propulsion  appears  feasible. 


(20)  Alan  Phillips,  S.A.I.,  private  communication 
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4.  PROPAGATION  THROUGH  HzO  PLUME 

Water  (or  ice)  appears  to  be  an  ideal  fuel  for  laser  propelled 
vehicles  because  a)  it  is  cheap,  b)  it  does  not  pollute,  and  c)  it  has  a low 
molecular  mass  in  the  vapor  phase  and  high  density  under  ambient  con- 
ditions (as  opposed  to  N2>  H^,  etc.  which  have  to  be  cooled). 

A drawback  to  the  use  of  water  as  opposed  to  liquid  for  example, 
is  the  condensation  that  will  occur  in  the  plume  and  the  deleterious  effect 
that  condensation  may  have  on  laser  propagation.  Also,  high  temperature 
H^O  is  an  absorber  of  10.  6 p m radiation.  We  will  examine  the  effect  of 
water  vapor  transmission  and  water  condensation  on  laser  propagation  in 
three  regimes. 

a)  High  altitude  plume  which  expands  into  vacuum. 

b)  Plume  at  20  km  altitude. 

c)  Low  altitude  plume  near  sea  level. 

In  Section  VL  4.  a.  we  summarize  the  data  on  absorption  by  water  vapor  in 
the  10.  6 pm  window.  In  Section  VI.  4.  b.  we  calculate  the  exit  conditions 
at  the  nozzle  for  a typical  vehicle  with  mass  8 tons,  acceleration  a g 
(a  3)  and  specific  impulse  800  sec.  The  high  altitude  plume  condensation 
is  calculated  in  Section  VI.  4.  c.  and  low  and  intermediate  altitude  plumes 
in  Section  VI.  4.  a.  The  absorption  of  the  laser  beam  in  these  plumes  is 
also  considered  in  these  sections.  We  present  our  conclusion  in  Section 
VL  4.  e. 

a.  Absorption  by  Water  Vapor 

At  low  temperatures  the  absorption  by  water  vapor  in  the  region 
near  10.  6 pm  is  small  and  is  of  the  continuum  type.  The  variation  of 
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absorption  with  wavelength  is  consistent  with  the  interpretation  that  the 
absorption  is  due  to  wings  of  lines  originating  from  the  pure  rotation  band 
of  H2O  beyond  15  (lm  and  of  lines  from  the  6.  3 (im  vib ration -rotation  band. 
A theoretical  estimate  of  absorption  made  by  Kaplan^1  ^ due  to  the  wings 
of  distant  lines. 


k (v)  = 


4.  1 X 10 

(v  - 200)' 


3.  7 X 10 
(v  - 1 550)' 


-1  2 
g cm  . 


(382) 


where  v is  the  wavenumber  in  cm  \ is  in  fair  agreement  with  measure- 

ments  of  continuum  absorption  in  the  atmosphere  made  by  Bignall  et  al,  ' 

(23) 

and  Roach  and  Goody.  ' 

Analyses^^’  of  more  recent  measurements  by  Burch^^V 
McCoy  et  al^^'  and  Arafev^^  indicate  that  continuum  absorption  is  due 
to  a water  vapor  dimer,  the  absorption  coefficient  being  given  by 


(for  X = 10.  6 (im) 


a = 4.3  x 10'6p(P+  139  p)  km-1 


(383) 


This  theoretical  estimate  is  quoted  in  the  experimental  paper  by 
Bignell  ct  al.  , next  reference. 

Bignell,  K.  , Saiedy,  F.  and  Sheppard,  P. , J.  Opt.  Soc.  Am.  53.  466 
(1963). 

Roach,  W.  , and  Goody,  R.  , Quart.  J.  Roy.  Meteorol.  Soc.  (London) 

84,  319  (1958). 

Roberts,  R. , Selby,  J.  , and  Biberman,  L.  , Appl.  Optics  ^5,  2085  (1976). 

Adiks,  T.,  Aref'ev,  V.,  and  Dianov-Klokov,  V.;  Sov.  J.  Quant.  Electron, 
5,  481  (1975). 

Burch,  D. , "Investigation  of  the  Absorption  of  Infrared  Radiation  by 
Atmospheric  Gases,  " Semi-Annual  Technical  Report,  Aeronutronic 
Division,  Philco  Ford  Corporation,  Report  U-4784  (January  1971). 

McCoy,  J.  , Rensch,  D.  , and  Long,  R.  , Appl.  Opt.  8,  1471  (1969). 
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where  p is  the  partial  pressure  of  water  vapor  (in  Torr)  and  P the  ambient 
pressure  (in  Torr),  or,  equivalently,  by 

a = 1.  76  X 10"3  p + 1.  42  X 10'6  p2  exp  (- AH/RT)  km'1  (384) 

3 

where  p is  the  water  vapor  concentration  in  g/m  and  AH  = -4,  500 
cal/mole.  The  second  term  in  Eq.  (384)  is  due  to  the  water  vapor  dimer 
and  dominates  over  the  first  term  for  p > 2.  3 g/m  at  sea  level  (12% 
relative  humidity  at  20°C).  The  term  A H in  the  exponent  correlates 
well  with  the  estimated  binding  energy  of  the  water  vapor  dimer.  If 
absorption  at  10.  6pm  is  due  entirely  to  a fundamental  band  of  the  water 
vapor  dimer  then  the  absorption  coefficient  at  higher  temperatures  should 
be  given  by  the  second  term  of  Eq,  (384).  A plot  of  the  absorption  length 
L = l/k  as  given  by  Eq.  (384)  as  a function  of  temperature  and  water 
vapor  pressure  (or  concentration)  is  shown  in  Figure  89. 

Experimental  measurements  of  water  vapor  absorption  at  temp- 
eratures T s;  2000°C  and  pressure  p ^ 1 atm  were  reported  by  Ludwig 
et  al.  ' ' At  high  temperatures,  vibration -rotation  lines,  due  to  high 

lying  levels,  appear  in  the  10  pm  window  and  are  the  main  cause  for 
absorption.  Unfortunately,  the  measurements  of  Ref.  28  were  taken 
with  a resolution  of  25  cm  ^ and  the  absorption  coefficient  at  the  location 
of  the  CC>2  laser  lines  is  unknown.  We  have  used  the  compilation  of 
spectral  information  on  made  by  McClatchey ^ 2<^'  to  estimate  the 

(28)  Ludwig,  C.  , Ferriso,  L.  , and  Abeyta,  C.,  J.  Quant.  Spectr.  and 
Rad.  Transfer  j>,  281  (1965). 

(29)  McClatchey,  R.  A.  , et  al. , "AFCRL  Atmospheric  Absorption  Line 
Parameters  Compilation,"  Air  Force  Cambridge  Research  Labora- 
tory, AFCRL-TR-73-0096,  January  1973. 


137 


absorption  coefficient  at  the  specific  wavelengths  of  the  CO^  laser.  The 
(29) 

McClatchey  tape  contains  information  on  lines  that  appear  at  room 

temperature  only,  and  represents  a lower  bound  on  H^O  absorption,  since 
other  lines  may  become  important  at  higher  temperatures.  Calculations 
of  the  absorption  coefficient  at  two  pressures  are  presented  in  Table  8. 
Absorption  is  the  largest  for  the  P (20)  line  of  CO^,  this  line  falling  only 
0.  18  cm  l away  from  a line  of  I-^O.  The  line  shape  dvie  to  collision 
broadening  is  given  by 


S(v)  = 


So<T>  Y 

* [ (v  - vo)2  + y2  1 


(385) 


where  vq  is  in  the  position  of  the  line  center  and  y,  the  line  width,  is 
proportional  to  pressure  and  varies  inversely  with  the  square  root  of 
temperature.  The  variation  of  line  intensity  per  atom  with  temperature 


is  given  by 


So  <T>  = 


exp  - 


(") 


(386) 


where  e is  the  energy  of  the  lower  state  of  the  vibration  rotation  transi- 
tion considered  and  Qy  and  are  the  vibrational  and  rotational  partition 
functions.  For  water  vapor  at  temperatures  below  3 000°K,  w T^2 

and  Qy  ~ (l  - exp  - 0/ T)  ' where  fl  = 2290°K  (corresponding  to  the 
6.  3 flm  band  of  H^O).  At  standard  temperature  and  pressure,  the  line 
width  of  an  H^O  line  is  estimated  to  be^^  ~ 0.  07  cm  *,  so  that  at 
pressures  less  than  atmospheric,  the  CO2  lines  fall  well  into  the  wings 
of  the  H2O  lines  compiled  by  McClatchey.  The  absorption  coefficient 

(30)  Goody,  R.  , Atmospheric  Radiation.  Vol.  1 Theoretical  Basis,  Oxford 
University  Press  (1964),  p.  115. 
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TABLE  8.  ABSORPTION  COEFFICIENT  OF  PURE  H,0  VAPOR  AT  v (cm 


a at  temperature  T and  pressure  P is  obtained  from  the  absorption  coef- 
ficient a at  a reference  temperature  T and  pressure  P by  combining 


Eqs . (385)  and  (386).  We  find 

3 


«o  * (pP0)2  ( T°)  exp['  • i)}\ 


exp  (-  0/T) 
exp  (-  e/T0) 


(387) 


where  one  factor  (P/P  ) (T  /T)  takes  into  account  the  varying  amount  of 

absorber  (assumed  pure)  in  the  path,  another  factor  (P/P  ) (T  /T)  takes 

into  account  the  pressure  broadening  of  the  lines  and  the  remaining  factors 

are  due  to  the  partition  functions.  Absorption  of  the  P(16)-P(24)  lines  of 

CC>2  by  the  BLjO  vapor  were  found  to  follow  Eq.  (387)  fairly  well  in  the  range 

300  to  2500°K  and  can  be  represented  with  an  accuracy  better  than  30%  by 

P2  , _i 

a = A — - exp  f-  (e/kT)]  [ 1 - exp  (-  2290/T)]  cm  (388) 

T"5 

where  P is  in  atmospheres  and  A and  e are  given  by  Table  9. 


TABLE  9.  PARAMETERS  ENTERING  IN  EQ.  (388) 


Line 

P(16) 

P(18) 

P(20) 

P(22) 

P(24) 

A(cm~  1 Atm'2  (°K)3) 

1.  53  x 106 

6 

4.  48  x 10 

5.  55  x 107 

4.  1 x 106 

5.  5 x 106 

£/k  (°K) 

1870 

2950 

3330 

3150 

3010 

A comparison  of  Eq.  (388)  and  Table  9 for  the  P(20)  line  of  CO^  is  shown  in 
Figure  90. 

b.  Exit  Conditions  for  the  Study  of  Plume  Behavior 

We  assume  for  generality  a pulsed  device  where  the  expansion  time 
of  the  vapor  in  the  nozzle  is  r (note  T » Tlaser>  and  the  repetition  rate  is 
f.  We  define: 
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r 


I = Specific  impulse  (800  sec)  = V£/g 
T = Temperature 
m = Mass  of  vehicle  (8  tons) 

-23 

ju  = Molecular  mass  of  fuel  (3  x 10  g) 

-2. 

g = Acceleration  of  gravity  (1000  cm  sec  ) 
z = Compressibility  factor 
u = Velocity  of  gas  with  respect  to  vehicle 
A = Nozzle  area 
p = Vapor  density 
rh  = Mass  flow  rate 

e = Subscript  referring  to  exhaust  conditions 
s = Subscript  referring  to  stagnation  conditions 
The  equation  of  motion  of  the  vehicle  averaged  over  many  pulses  is: 

m IF  = ™ Ue  = {Pe  Ue  V Tf  Ue 


We  require,  furthermore, 


dV 

m ~dt  = a m g 


Where  typically  a - 3. 

Solving  for  pg  we  get 


a m 

Pe  I2  g rf  A 


(389) 


We  supplement  Eq.  (389)  with  the  condition  that  the  efficiency  of  the  engine 
is  q . 

The  ratio  of  exit  enthalpy  hp  «o  stagnation  enthalpy  hfi  is  given  by 


tt-  = l - 9 


(390) 


l 
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t I 


tLl- 


also 


u 2 2 2 

h = -JS_  = 1 -B- 
e 2 2 


(391 


For  I and  n given,  h and  h are  uniquely  determined  by  solving  Eqs.  (390) 
and  (391).  We  use  the  Mollier  diagram  for  H^O  presented  in  Figure  91. 

We  assume  isentropic  expansion  in  the  nozzle.  This  means  that  the 
stagnation  and  exit  conditions  lie  on  a vertical  line  s = const,  in  the  Mollier 
diagram.  In  order  to  be  definite,  we  choose  the  reasonable  (but  somewhat 
arbitrary)  stagnation  conditions. 


p =80  atm 

T = 5700°K 
s 


Stagnation  Conditions 


= 5.  3 x 10  j/g 


We  then  calculate  that  for  a 50%  efficient  engine  (q  = 0.  5) 


0.  06  atm 


T = 3000  K 
e 

h = 2.65  x 1CT 
e 

z =1.6 
e 


Exit  Conditions 

(point  of  B of  Figure  91 ) 


I = 760 

The  above  conditions  seem  reasonable  for  operation  at  high  altitudes 
where  p < 0.06  atm,  but  would  give  a very  inefficient  engine  at  sea  level. 
One  would  get  flow  separation  in  the  nozzle  and  shock  formation  that  would 
bring  the  exhaust  pressure  up  to  atmospheric.  The  computation  of  the  exit 
conditions  for  an  over-expanded  jet  taking  into  account  all  these  effects  is  a 
major  task.  We  will  here  make  the  drastic  assumption  that  the  fluid  in  the 
nozzle  expands  up  to  the  point  p = 1 atm  and  then  forms  a parallel  jet  which 
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2000 


ENTROPY  (cal  /gram  °K) 


91  Molli rr  Chart  for  H?0 


eventually  mixes  with  the  atmosphere.  This  will  result  in  a jet  cross  sec- 
tion that  is  much  less  than  the  nozzle  cross  section  (such  a jet  can  be  realized 
by  using  a double  nozzle).  The  efficiency  of  the  engine  r)  and  thrust  I will 
be  somewhat  reduced.  The  Mollier  diagram  for  the  above  mentioned  stagna- 
tion conditions  yields  (point  C of  Figure  91). 
h = 3.  2 x 104  J/g 

e t 

T = 3600°K 


p =1  atm 
*e 

z =1.9 
e 

n = 0.4 


Exit  Conditions  at  Low  Altitudes 


1 = 650  sec 

In  order  to  determine  the  exit  area,  we  solve  Eq.  (389)  for  A using 
p as  given  by  the  Mollier  diagram.  The  smallest  area  corresponds  to  t f 
1,  i.e.,  cw  operation.  We  find  for  the  high  altitude  case  (with  rf  = 1 ) 

pe  = 1 * 35oo  x x 10  x 29  x 1 9 = 3.2  x 10  g/cm 


„ 3x8x10 

A = r 

(650r  X ( 1 0 0 0 ) x3 . 2 x 10 


= 1.7  x 10'^  cm^ 


The  exit  conditions  calculated  above  for  the  two  cases  are  listed  in  Table  10. 
We  look  at  condensation  in  the  plume  and  its  effect  on  laser  propagation  in 
the  next  sections. 

c.  Absorption  in  High  Altitude  Plumes 


For  the  high  altitude  case  we  assume  that  the  plume  is  expanding  into 
vacuum,  that  the  chemistry  is  frozen  and  make  a constant  y = 1.33  assump- 
tion for  the  gas  (frozen  vibration,  3 degrees  of  freedom  for  rotation  of  H-,0). 
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The  structure  of  highly  under  expanded  jets  has  been  investigated  by  many 

.(3 1) 


authors  and  we  summarize  the  results  relevant  to  this  research  below. 


The  jet  geometry  is  shown  in  Figure  92.  Sufficiently  far  away  from  the  noz- 


zle the  flow  is  practically  radial  and  the  density  varies  as  r~2.  Also  at  a 


given  radial  position  the  density  varies  as  [cos  )|  2/^Y  ' 1 * where  9 

L.  (l  ( 


is  the  Prandtl  Mayer  expansion  angle  and  9 the  polar  angle  with  respect  to 
the  axis  of  the  vehicle.  The  results  quoted  above  do  not  take  into  account 


the  effect  of  condensation.  Since-the  enthalpy  of  condensation  (11  kcal/mole 


3 

3 x 10  j/g)  is  ten  times  less  than  the  kinetic  energy  of  the  expanding  plume 


2 / 4 

(ue  / 2 - 3 x 10  J/g)  we  expect  that  condensation  will  affect  the  flow  in  a 


negligible  way. 


( 1 ) Study  of  Condensation  and  its  Effect  on  Laser  Propagation 


.-4 


The  characteristic  time  in  our  system  is  r = '\rK~/  u ^ 10'  sec 

e'  e 


and  condensation  will  start  occurring  when  the  condensation  rate  is  high 
enough  so  that  an  appreciable  fraction  of  the  flow  would  condense  over  the 
time  t. 

The  rate  of  formation  of  condensation  nuclei  of  critical  radius  r , 
following  condensation  theory,  is  given  by 

1 /zu«v/2 

2 


r 

P\\- 


^LUa^j 


X exp 


$ 2 ” 
4 n ( r ) ct 

3 kT 


i • -3  -1 

nuclei  cm  sec 


392) 


.vhere  p is  the  density  of  condensed  phase  (1  g/cc),  a is  the  surface  tension 


for  ice  (96  dyne/cm)  and  where  the  critical  radius  of  a condensation  nucleus  is 
related  to  the  super  saturation  ratio  as  follows 


(31)  H.  Ashkenas  and  F.  Sherman  in  Experimental  Methods  in  Rarefied 
Gasdynamics , J.  H.  de  Leeuw,  Academic  Press,  New  York  (1966) 
Vol.  2,  p.  84;  F.  Albini,  A1AA  J.  , 1535  (1965);  F.  Boynton,  A1AA 

J.  , 5,  1703  (1967). 
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* 2 P a 

A — rr-»  /I  __  / . 


P kT  Pn  (p/p  .) 


The  factor  I'  varies  over  several  orders  of  magnitude  depending  on  the  nuclea- 
tion  theory  considered.  Because  of  the  strong  temperature  dependence  of  the 
exponent,  however,  the  actual  value  of  F turns  out  to  be  not  too  important.  In 
our  calculations  we  take  r = 1 which  should  be  a good  approximation  for  water 

vapor.  We  obtain  the  condensation  rate  by  multiplying  FJq.  (392)  by  the  factor 

*3 

4ir/i  (pw/p)  r which  is  number  of  water  molecules  in  a critical  condensa- 
tion nucleus.  The  saturated  vapor  pressure  of  ice  is  approximately  given  by 


p : 4 x 10^  exp  (-6170/T)  atm 
sat 

Assuming  that  the  expansion  of  the  gas  is  adiabatic,  i.e.  , 


(394) 


= 1 = 


(395) 


we  can  relate  T and  p to  p.  Using  Eqs.  (392)  and  (393)  and  the  exit  condition 
of  Table  10  we  obtain  the  results  shown  in  Figure  93.  We  plotted  in  Figure  93 
three  curves  as  a function  of  Mach  number  in  the  flow.  The  first  two  curves 
are  temperature  and  saturation  ratio  and  show  that  the  flow  becomes  saturated 
when  T = 200°K.  The  third  curve  represents  the  ratio  of  mass  condensed  per 
unit  volume  [as  given  by  Eq.  (392)]  during  the  characteristic  time  t = \rA  / u 
to  the  mass  per  unit  volume  p.  Whenever  this  ratio  reaches  a number  of  order 
unity  we  expect  condensation  to  occur.  From  Figure  9 3 we  see  that  condensa- 
tion occurs  at  T = 150°K,  i.e.  , 50°K  below’  the  saturation  temperature.  We 
therefore  take  r = r as  the  radial  position  where  T = 150°K  and  pc  = 3.  5 x 10 
g/cm  \ In  order  to  make  calculations  of  transmission  we  assume  that  all  the 

More  prec isely , one  should  use  r = r /u^,  see  discussion  next  page. 
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Figure  93  Condensation  in  High  Altitude  Plume 

a)  Temperature  vs  Mach  Number 

b)  Supersaturation  Ratio  vs  Mach  Number 

c)  Normalized  Condensation  Rate  vs  Mach  Number 


vapor  is  condensed  for  r > r^.  This  should  yield  an  upper  bound  on  absorp- 
tion by  the  plume.  Absorption  will  be  maximum  for  a beam  traveling  along 
the  axis  of  the  vehicle.  The  integrated  water  content  per  unit  area  along  the 


axis  is  given  by 


or> 

-f 


P (r)  d r = pc  rc 


/(£)r-'c%  /(^)Vp. 


where  the  subscript  c refers  to  conditions  at  the  condensation  point.  The 
radial  distance  r^  is  determined  by  using  continuity  of  flow 


nuA  = n u nr  x 
fe  e e Kc  c 


/ [“*(*£)] 


2/(Y  - 1) 


sin  0 d 0 


For  a Mach  number  4.45  flow  with  y = 1.33,  we  have  0 = 69  and  the 

integral  is  0.08.  Letting  pc  = 3.5  x 10  ^ g/cm3,  u^  a<  u^  = 7.6  x 10^ 

cm/sec  and  using  the  exit  conditions  of  Table  8 we  find  = 140  m and 

- 6 Z 

w = 5 x 10  g/cm  . 

We  calculate  the  absorption  length  due  to  condensed  water  assuming 
that  all  droplets  are  less  than  10  pm  in  diameter  (Rayleigh  region).  The 
absorption  cross  section  per  droplet  of  radius  R is: 

a = n R2  Q , 


where  Q . at  10.  6 am  is  shown  in  Figure  94.  One  sees  from  Figure  94 
abs  ' 

2 

that  Q , =3x10  R where  R is  in  cm.  Adding  the  contribution  of  all 

dbs 

particle  j in  a unit  volume  we  obtain  the  absorption  coefficient  a 


a = 3 


x io2  ’rR 3 = t-  3 * 1-°-  y;  ^ r. 

L-J  J 4 w p L-*i  3 j 


= 70  p cm 


(396) 


where  p = density  of  water  ( = 1 g/cm  ) and  p is  the  mass  concentration 
of  water  in  g/cm3. 

* This  assumption  will  give  maximum  absorption. 
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The  transmissivity  of  the  condensed  part  of  the  plume  to  a 10.  6 pm 

laser  beam  traveling  along  the  axis  is  thus 

rp  -70  W , . c -4 

T = e = 1 - 3.5  x 10  =*  1 

Condensation  should  have  therefore  a negligible  effect  on  laser  propagation. 

As  a closing  remark  we  note  than  the  characteristic  time  for  conden- 

4 

_ „ ; ..  l . , 1 • • „ 10  „ . n“  2 


sation  along  the  axial  stream  line  is  r /u  =* — » 1 0~  sec  rather 

_4  C c 8 x 10' 

than  10  sec  or  stated  earlier.  This  does  not  effect  our  results  significantly, 


as  can  be  seen  from  Figure  93(c),  where  a shift  in  ordinate  by  a factor  of  10^ 
modifies  (and  thus  T^)  by  only  a small  amount. 

(2)  Absorption  by  Water  Vapor  Before  Condensation 


We  have  just  found  that  condensation  in  the  plume  at  high  altitudes  will 
affect  laser  propagation  in  a negligible  way.  We  now  study  the  effect  on  laser 
transmission  of  water  vapor  content  in  the  plume  between  the  exit  nozzle  and 
the  condensation  point.  We  consider,  as  we  did  before,  the  stream  line  along 
the  axis  of  the  vehicle  which  contains  the  largest  amount  of  vapor.  The  vapor 
starts  expanding  when  the  Mach  line  from  the  edge  of  the  nozzle  intersects  the 
axis  of  the  vehicle.  This  occurs  at  a distance. 

L = R M 

e 

We  thus  have  T = T^  and  p = pe  for  x < L.  For  x > L we  approximate 


the  flow  field  by 


P (*>  = Pe 2 

e (x  - L)  + 1/ 


(397) 


and  use  the  isentropic  relation  (394)  to  calculate  T and  p.  We  can  then  cal- 
culate absorption  due  to  high  temperature  lines  of  H.,0  by  use  of  Eq.  (397). 
We  have  calculated  the  transmission  through  the  plume  given  by 


r = 


exp 


J 


in  the  following  two  cases. 

(a)  Frozen  Vibration,  Frozen  Chemistry 
Let  the  absorption  coefficient  a = at  the  exit  nozzle.  Then  at  a 
location  x the  absorption  coefficient  is 

a 

3 - Y 


a (x) 


We  have 


L 


/V.  ^ /»L 

° d , = j oc  d X + J O d X 


= a L + 
e 


% l / 


X /L 
c' 


d (x/L) 


(x/L  - if  + 1 ] 


3 -V 


=*  2 a L 


where  in  the  last  step  we  made  x -►  <*  and  assumed  3 - y = 1.5. 


(b)  Equilibrium  Vibration,  Frozen  Chemistry 


/“ 


d x = a L + 


ne  L 

e + 2 | 1 - exp  (-  0/T  )] 

e T /T 

c e 


/ dt^~rri tl_v  exp[- 


k T 


- (t1" 


x [ 1 - exp  (-  0 t 1 ~ Y)] 

Numerical  integration  of  the  integral  for  « T^,  \ = 1.33  and  T^  = 
3000°K  yields  a value  of  0.  78. 


/■ 


n d x = n L <1  + — 


0.  78 


0 - exp  -0/T) 


= l . 88  o L 
e 


The  above  result  is  for  all  practical  purposes,  the  same  as  for  frozen 


vibration . 


i e c. 


We  evaluate  the  absorption  for  the  P(20)  line  of  CO.,.  We  have 


L = R M = 
e 


5 x 10 


« = x ,0.06,2xexp|- 

e I3000)3 


.-4 


( 333o\/, 
3000  A 


x 4.5  = 300  cm 


1 - exp  - 


2890 
3000  J 


^ = 1.5x10  ^ cm  ^ 


Thus  o L = 4.  5 x 10"  and  absorption  by  water  vapor  is  negligible, 
e 

d.  Absorption  and  Condensation  in  Lower  Altitude  Plumes 

At  low  altitudes  the  external  pressure  prevents  the  jet  of  vapor  from 
expanding  free’  ’ and  cooling  of  the  jet  arises  from  entrainment  of  and  mixing 
with  the  ambient  air.  The  cooling  of  the  jet  and  possible  condensation  de- 
pends on  a)  the  amount  of  air  entrained,  b)  the  ambient  temperature,  and 
c)  on  the  degree  of  recombination  of  the  partially  dissociated  water  vapor. 

We  have  used  a three  dimensional  axisymmetric  computer  code  which 
numerically  solves  mass,  momentum,  enthalpy  and  species  conservation 
equations  using  turbulent  diffusivity  and  finite  rate  chemistry.  The  results 
of  two  runs  at  20  km  altitude  (p  = 0.  06)  and  sea  level  (p  = 1 atm)  were  made 
and  will  be  discussed. 

( l ) Condensation  of  Water  Vapor  in  Plumes 

Before  presenting  results  of  the  computer  runs,  we  first  discuss  the 
possibility  of  condensation  in  the  plume  and  show  from  general  energy  consid- 
erations that  condensation  is  unlikely  to  occur  if  the  air  temperature  is  larger 
than  » 230°K. 

Consider  a mass  rn^  of  plume  material  that  entrains  a much  larger 
mass  ni  of  air.  Assume  that  the  exit  pressure  is  the  same  as  the  ambient 


pressure  p . Conservation  of  momentum  and  energy  yields 
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(398) 


in  u = (m  + m ) u 
p e a p 

m (h  + 4"  u *")  + m h (T  ) = m h (T)  + m h (T)  + \r  (m  + m ) u2. 
p'e  2 e a a o pp  a a ' 2 ' p a 


(399) 


We  assumed  in  the  above  that  the  air  was  at  rest  and  that  the  vehicle  velocity 
is  much  less  than  u^,  so  that  the  relative  velocity  of  the  plume  with  respect 
to  the  air  is  ug.  Tq  is  the  ambient  temperature  and  T is  the  temperature  of 
the  final  mixture.  Eliminating  u between  Eqs.  (398)  and  (399)  yields 


h +Iu  2 - h (T)  = (-^~A  — 

e 2 ' p V *P  / 


C AT  flu  2 
pa  2 e 


(400) 


4-  +(W  ) — 

p p 


where  AT  = T - Tq  is  the  temperature  rise  of  the  air,  C^a  is  the  specific 

heat  of  air,  4 is  the  molar  concentration  of  the  plume  constituents  in  the 

mixture,  and  p.a  and  fj  are  the  molecular  weights  of  air  and  plume  material 

respectively.  For  mixtures  in  which  T =*  we  can  neglect  h (T)  as  com- 

1 2 

pared  to  the  stagnation  enthalpy  h^  + ^- u^  , anc^<  as  we  shaH  see  later,  we 
can  neglect  the  second  term  on  the  right  hand  side  of  Eq.  (400).  We  therefore 


1 2 (l  - V _ . _ 

h = h + u = — C _ AT 

s e 2 e </p  Mp  Pa 


Assuming  that  recombination  is  complete  and  ij1  « 1,  we  have 


p h2o 


C AT 
na 


,lH  ,0 


(401) 


Condensation  will  occur  only  when 


p HzO  Z Psat 


(402) 


Where  p is  given  by  Eq.  (395) 
sat 


Psat  = B e 


B = 4 x 10  atm 


.0  = 6170  K 


(403) 


We  write  T = + AT  with  AT  « T^  and  expand  the  exponent  in  Eq.  (394). 

We  look  for  conditions  under  which  inequality  (402)  is  just  satisfied.  This  will 
yield  the  temperature  T and  concentration  t/-  at  which  condensation  may  (if 
the  time  scale  is  long  enough)  just  occur.  Using  Eqs.  (401)  and  (394),  Equality 


(4021  gives 


h B a -0/  T 
s 0 e ' o 


2°  ^ Cr 


2 P„ 


(404) 


where  we  defined  AT# /T ^ = y.  Now,  the  left  hand  side  of  Eq.  (398)  has  a 

maximum  at  y = 1 of  e ^ = 0.  37.  If  Ty  is  so  large  as  to  make  the  right 
hand  side  of  (404)  larger  than  0.  37  then  there  is  no  solution  and  condensation 
will  not  occur.  We  thus  write  the  condition  for  any  condensation  to  occur  as 


P'H-O  hs  B 


^a  pa  T 


0 e'^o 

r 2 


< 0.  37 


(405) 


We  discuss  the  range  of  validity  of  Ineq.  (405),  or  rather  Eq.  (404) 
from  which  it  was  derived.  We  assumed  that  hg  represents  the  stagnation 
enthalpy  of  the  plume.  If  the  chemistry  is  slow  (i.e.  , recombination  is  in- 
complete) and  radiation  losses  are  important,  one  should  subtract  from  h^ 
the  amount  of  energy  remaining  in  dissociation  and  the  energy  per  unit  mass 
lost  by  radiation.  In  order  to  take  into  account  the  dissociation  remaining, 
the  right  hand  side  of  Eq.  (404)  should  be  multiplied  by  1 /f  where  f is  the 
fraction  of  H^O  (f  = 1 if  recombination  is  complete).  The  computer  runs 


that  we  have  made  at  p = 0.  06  atm  indicate  that  less  than  25%  of  the  initial 
enthalpy  of  dissociation  remains  when  a significant  amount  of  air  has  been 
entrained.  Also,  under  the  times  of  interest  (t  < 0.  1 sec)  we  have  found 
that  losses  by  radiation  are  insignificant. 

The  air  temperature  for  which  inequality  (405)  is  satisfied  is  shown 
in  Figure  95.  In  order  to  plot  the  boundaries  in  Figure  95  we  used  = 
3000°K,  ug  = 7.  6 x 105  cm/sec,  hg  = 5.3  x I04  J/g,  Cpa  = 1 J/g°K, 
o = 29/l8  = Several  boundaries  are  indicated  in  Figure  95, 

corresponding  to  enthalpy  available  to  heat  the  gas  being  100%,  50%,  or 
20%  of  the  initial  stagnation  enthalpy. 

The  temperature  rise  of  the  air  when  condensation  occurs  is  AT  = 

2 

y T /O  where  y is  the  solution  of  Eq.  (404).  The  maximum  rise  in  tempera- 
ture corresponds  to  the  maximum  value  of  y,  namely  y = 1.  Therefore, 


(A 


TV  o (200)  _ z°K 

max  0 6170 


The  partial  pressure  of  H^O  at  condensation  is  also  a maximum  when  y = 1, 
From  Figure  95  one  would  expect  condensation  to  occur  in  a small  altitude 
range  around  12  km  where  pQ  = 0.2  atm,  T^  = 220  K.  The  concentration  of 
condensed  and  uncondensed  water  under  the  above  conditions  would  be 

P = (Pair>- 


STP 


||^  x B exp  (-  0/T) 


= 1.2  x 10"3  HI  x x 4 x 107  exp  - (6170/226) 

= 5 x 10^  g/cm3 . 

1 he  Absorption  coefficient  at  10.6  jjTr.  due  to  cor.der.6ed  water  is  give  t by 
t q.  v Wt>)  and  will  be  at  the  most  3.  f \ 10  ° cm  “.  We  can  r.eglect  the 
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Figure  95  Map  for  Plume  Condensation  vs  Altitude 
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effect  of  condensed  water  on  laser  propagation,  therefore,  unless  we  are 
dealing  with  path  lengths  through  condensed  zones  in  excess  of  1 km.  We 
can  now  show  that  the  neglect  of  the  second  term  on  the  right  hand  side  of 
Eq.  (400)  as  compared  to  the  first  one  was  justified.  Since  if  « the 

neglect  of  the  second  term  on  the  right  hand  side  of  Eq.  (400)  is  justified  if 

2 ? 

v C AT 

...Pa 


we  havt 


•-■•(£) -fell 

pV1  j ^ ^4  x 107  e^-mo/23°y  = 


(406) 


-7 


= 2 x 10 


t-hile 


(y  ^ 


AT 


= 2 x (1.6)  x 


10  x 6 


= 5 x 10 


-4 


e (8  x 10') 

Inequality  (406)  is  therefore  well  satisfied. 

(2)  Plume  at  20  km  Altitude 

The  results  of  the  three  dimensional  code  for  a plume  with  the  exit 

conditions  given  by  the  first  line  of  Table  10  are  shown  in  Figure  96.  The 

ambient  air  temperature  is  T . The  contours  of  constant  H,0  concentration 

are  found  to  follow  the  contours  of  constant  temperature  as  is  to  be  expected. 

The  computer  run  was  lengthy  and  the  results  of  5 hours  computation  on  a 

CDC  3174  led  to  a down  stream  distance  of  80  exit  diameters  (x  = 100  meters). 

Contours  beyond  this  distance  were  extrapolated. 

The  temperature  and  concentration  of  H^O  along  the  axis  of  the  plume 

are  shown  in  Figure  97.  Since  the  self  broadening  of  H^O  is  (roughly)  five 

(32) 

times  larger  than  broadening  due  to  collisions  with  air  molecules,  we 

compute  the  absorption  coefficient  as  follows: 

iU)  Burch,  D.  , Singleton,  E.,  and  Williams,  D.  , Appl.  Optics  J_,  359(1962). 
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° (po’  '•  ' H-,0}  " n (po’  1 ’ ^H,0  n ^ 1L,0  4 11,0  + 
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1 - 4 


(407) 


where  V ,,  is  the  concentration  of  H_0  and  n (p  , 1,  (/  ,,  1)  is  tabu- 

H-O  2 'ro  11,0 

C c 

lated  in  Table  2.  We  used  Table  9,  Figure  97  and  Eq.  (407)  to  evaluate  the 

absorption  coefficient  to  the  P(20)  line  of  CO,  along  the  axis  of  the  plume 

where  absorption  is  the  largest.  The  absorption  coefficient  is  shown  in 

Figure  9b  as  well  as  its  integral  up  to  ij  ^ q ~ 100  m.  The  absorption 

1 00  1 


f n dx 


4 x 10 


i.e,,  is  negligible.  One  should  keep  in  mind  that  we  used  the  McClatchey 
tape  to  compute  water  vapor  absorption  and  that  broad-band  measurements 
at  2000°K  indicate  absorption  100  times  larger^*”'  ' than  that  given  by  the 
McClatchey  tape.  Even  with  this  factor  of  100  taken  into  account,  absorption 
over  a length  of  100  meters  of  plume  is  small. 

In  spite  of  the  fact  that  we  have  calculated  a very  small  absorption 
coefficient,  absorption  may  still  be  significant  if  the  laser  is  directed  straight 
up  the  axis  of  the  vehicle  and  has  to  traverse  the  whole  length  of  the  plume 
which  may  be  quite  long.  One  could  eliminate  these  very  long  paths  by  re- 
quiring a small  tilt  ( ~ l-S°)  of  the  vehicle  with  respect  to  the  laser  beam. 

Also,  the  presence  of  a transverse  wind,  which  usually  exists  at  high  alti- 
tudes, will  cause  a bending  of  the  plume  and  shorten  the  path  length.  We 
consider  below  the  effect  of  such  a transverse  wind. 

Assume  that  the  axial  plume  velocity  is  U and  that  the  transverse 
wind  velocity  is  W.  Consider  a unit  axial  length  of  plume  material.  The 
equation  of  motion  for  the  transverse  velocity  V of  this  element  of  plume  is: 


G5«6  AXIAL  DISTANCE  (m) 


Figure  ^8 


Absorption  by  20  km  Altitude  Plume  for  Axial  Path 
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* Pp  T dF  = CPa  (V  - W)  T 

where  C is  a drag  coefficient  of  order  1.  Integration  from  0 to  t,  with  the 
boundary  condition  V (t  = 0)  = 0,  yields 


V = W 


1 


1 


1 + t J 
o 


where  tQ  = p^  r D/(2  C Wp^). 

Integrating  once  again,  we  get  the  transverse  y coordinate  of  the 


y = W 


*•  ln(‘  +t). 


(408) 


We  have  plotted  in  Figure  99  the  velocity  of  a fluid  element  on  axis  as  given 
by  the  computer  run,  versus  time.  A good  approximation  to  this  velocity  is 


u = u t < T 
e 


(/)• 


(409) 


t > r 


5 -4 

where  u = 7.6  x 10  cm  sec,  T = 6.5  x 10  sec,  a = 0.43. 
e 

Integration  over  time  yields  the  axial  position 

x = u t < T 
e 


M (?)'] 


(41  01 


t > T 


Using  the  above  expressions  for  x and  y we  have  plotted  in  Figure  100  the 

position  of  a laser  beam  with  respect  to  the  plume  when  the  laser  beam  is 

pointed  straight  up  the  axis  of  the  vehicle  and  the  plume  is  submitted  to  a 

cross-wind  of  velocity  W = 10  cm/sec.  We  chose  in  the  numerical  calcula- 

- S «,  A 

tions  D = 120  cm,  p^  =8x10  g/cc,  p^  ~ Pc  = 2.  8 x 10  g/cm  , 

166 


PH ""  ■ — 


C = 1,  resulting  in  tQ  = 6.  6 x 10  sec.  The  laser  beam  under  such  con- 
ditions traverses  only  a section  of  plume  300-400  meters  long.  The  absorp- 

_4 

tion  would  be  at  most  8x10  if  we  use  the  results  of  Figure  98. 

(3)  Plume  at  Sea  Level 

Finite  rate  chemistry  could  not  be  included  in  the  plume  entrainment 
program  at  low  altitudes  because  the  chemistry  is  too  fast  and  the  axial  step 
size  required  in  order  to  follow  accurately  the  chemistry  would  require  exces- 
sively long  computer  times.  We  therefore  ran  the  program  using  frozen 
chemistry.  Evaluation  of  the  reaction  rate  at  sea  level,  however,  shows  that 
the  chemistry  must  be  near  equilibrium.  We  therefore  used  the  temperature 
and  concentration  profiles  furnished  by  the  frozen  chemistry  run  and  then 
modified  these  profiles  by  calculating  the  equilibrium  concentrations  of  vari- 
ous constituents,  keeping  the  total  enthalpy  constant. 

The  exit  composition  of  the  plume  was  as  follows:  tp ^ = 0.25,  4 qj.j  = 

0.15,  4>h  = 0.15,  4 0 = 0.12,  V 0 = 0.06  and  Q = 0.23.  Other  pa - 

2 2 2 
rameters  can  be  found  in  the  second  line  of  Table  3.  The  enthalpy  of  recom- 
bination of  the  vapor  components  at  3 500°K  to  form  one  mole  of  11,0  was 
calculated  from  the  JANAF  tabulations'  to  be  74.  13  Kcal.  Similarly  at 
298°K,  the  enthalpy  of  recombination  was  calculated  to  be  73.  4 Kcal/ mole 
of  H^O.  It  is  reasonable  to  assume  that  the  enthalpy  of  reaction  should  be 
fairly  constant  over  the  whole  range  of  temperature  encountered.  For  final 
temperatures  below  2000°K  practically  all  of  the  plume  material  should  be 
recombined  into  H^O.  The  plume  program  that  we  used  allows  for  cooling 
and  heating  of  all  components  at  constant  total  enthalpy  but  with  frozen  chemi- 
stry we  took  the  temperature  given  by  the  plume  program  and  raised  it  by  an 
amount 

(33"j  JANAF  Thermomechanical  Tables,  Second  Edition,  U.S.  National  Bureau 
of  Standards  NBS-37,  Government  Printing  Office  (1971). 
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AT  = X Cp  (H20)  + (1  - x)  Cp  (air)  x 10  (41  1 > 

where  x is  the  fraction  of  recombined  in  the  final  mixture,  and  Cp  is 

the  specific  heat  in  cal/mole  °K  [ Cp  (H^O)  - 13  cal/mole  °K,  Cp  (air)  » 

7 cal/mole  °KJ.  Equation  (411)  should  be  accurate  for  T + AT  < 2000°K. 

At  higher  temperatures  we  used  air  tables  and  the  Mollier  diagram  (Figure 
91)  for  water  vapor  since  partial  dissociation  has  to  be  included. 

Temperature  and  water  vapor  concentration  profiles  of  the  low  alti- 
tude plume  are  shown  in  Figure  100.  The  temperature  and  water  vapor 
concentration  along  the  axis  is  shown  in  Figure  101.  Absorption  by  water 
vapor  for  an  axial  path  length  was  calculated  following  the  procedure  outlined 
in  the  previous  section  and  is  shown  in  Figure  102.  The  integrated  absorp- 
tion for  an  axial  path  length  extending  100  meters  from  the  exit  nozzle  of  the 
vehicle  is  found  to  be  4%.  Longer  path  lengths  do  not  lead  to  much  larger 
absorption  because  of  the  rapid  decrease  in  temperature.  The  effect  of  a 
10  m/sec  cross  wind  on  absorption  was  considered  and  is  shown  in  Figure 
100.  It  sould  reduce  absorption  only  slightly  since  most  of  this  absorption 
occurs  in  the  first  40  meters  of  plume  and  the  plume  is  appreciably  bent  by 
the  wind  only  beyond  40  m.  Basing  oneself  on  the  II.-,  O absorption  lines  tabu- 
lated  by  McClatchey  one  would  not  expect  the  absorption  by  the  plume  to 
be  significant  at  sea  level.  The  broadband  measurements  of  Ludwig  et  al. 
indicate  that  absorption  could  be  as  much  as  100  times  higher,  leading  to 
catastrophic  absorption  by  the  plume  at  sea  level.  High  resolution  measure- 
ments of  H.,0  absorption  in  the  region  of  the  CO^  lines  is  needed  before  any 
firm  conclusions  can  be  drawn.  If  the  laser  beam  is  incident  at  an  angle  of 
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Figure  101  Temperature  and  H^O  Concentration  Along  Axis  of 
Plume  at  Sea  Level 
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10  with  the  axis  the  absorption  can  be  reduced  by  a factor  of  5 to  10,  but 
the  absorption,  following  the  measurements  of  Ref.  28,  would  still  be  signi- 
ficant 50%). 

e . Summary  and  Conclusions 

We  have  looked  at  condensation  and  absorption  of  CO^  radiation  in 

water  vapor  plumes  at  different  altitudes.  In  the  high  altitude  case  the  plume 

was  found  to  exapnd  and  cool  leading  to  condensation  when  the  temperature 

was  1 50°K  (50°K  below  the  saturation  temperature).  Absorption  by  the  con- 

_4 

densed  I-^O  was  found  to  be  negligible  (=*  10  ) and  absorption  by  the  vapor 

-4 

was  also  very  small.  (5x10  on  axis).  For  intermediate  and  low  altitude 

plumes  where  the  plume  does  not  expand  freely  but  mixes  with  the  surrounding 

air  it  was  found  that  condensation  should  not  occur,  except  possibly  for  a 

small  altitude  region  near  10  km  altitude  if  the  ambient  temperature  is  below 
o 

- 220  K (the  exact  temperature  depends  on  the  amount  of  enthalpy  stored  in 
dissociated  plume  products  that  have  not  recombined).  Even  if  the  condesa- 
tion  were  to  occur,  the  condensed  water  content  would  be  s*  5 x 10  g/cm  , 
i.  e.  , too  small  to  produce  significant  absorption. 

A three  dimensional  axisymmetric  plume  computer  code  was  used  to 


compute  temperature  and  H^O  concentration  profiles  in  water  vapor  plumes 
at  20  km  altitude  and  at  sea  level.  Using  these  profiles  we  estimated  that  the 
effect  of  H^O  absorption  on  the  P(20)  CO^  laser  line  propagation.  Absorption 

_4 

was  found  to  be  negligible  at  20  km,  (4x10  ) and  small  at  sea  level  (4%),  if 

one  used  the  H-,0  lines  tabulated  by  McClatchey  in  Ref.  29.  Measurements  on 
emissivity  of  H^O  flame  products  indicated'  ' an  absorption  100  times  larger 
than  those  quoted  above,  leading  us  to  believe  there  may  well  be  catastrophic 
absorption  in  sea  level  plumes.  The  measurements  of  Ref.  28  are  broadband, 
however,  and  high  resolution  transmissivity  measurements  at  the  exact  location 


t 


of  the  CO^  lines  are  needed  before  definite  conclusions  may  be  drawn. 
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5.  IGNITION  OF  AIR  BREAKDOWN  WAVES 

Fnder  the  high  laser  fluxes  considered,  solid  aerosols  in  the  atmos- 
phere can  heat  up,  vaporize  and  act  as  initiation  centers  for  air  breakdown 

waves.  These  waves  are  of  tw'O  types,  subsonic  laser  supported  deflagra- 

(34  35) 

tion  waves  ’ (LSA  waves)  that  are  semi-transparent  to  the  laser  beam 

/ o ZL  \ 

and  laser  supported  detonation  waves'  (LSD  waves)  that  are  supersonic 
and  totally  absorbing  to  the  laser  beam.  The  second  type  of  waves  are  of 
particular  concern  since  they  spread  in  a plane  perpendicular  to  the  beam 
at  a velocity 


v 


1.1x10 


t1/3 


cm/  sec 


where  I is  in  W/cm  and  the  numerical  coefficient  obtains  at  sea  level. 
Letting  1=  10  W/cm  , one  finds  that,  after  10  /jsec,  for  each  particle 
having  ignited  an  LSD  wave,  a 17  cm^  section  of  the  beam  will  be  blocked 
out.  The  formation  of  air  breakdown  waves  will  require,  however,  that 
the  particle  igniting  the  wave  have  a size  larger  than  a threshold  size. 
There  will  also  be  a delay  time  before  ignition  occurs.  Ignition  of  air 


(34)  Bunkin,  F.  , Konov,  V.I.,  Prokhorov,  A.M.,  and  Fedorov,  V.B., 
Zherf  Pis.  Red  _9,  609  (1969)  (English  translation  JE  TP  Letters,  9 , 
371  (1969). 

(35)  Raizer,  P,,  Yu,  Zhetp  Pis.  Red  Id,  195  (1969)  (English  translation 
JETP  Letters  11,  120  (1970). 

(36)  Raizer,  P.  , Yu,  Zhetp  (USSR)_48,  1508  (1965)  (English  translation, 
Soviet  Physics  JETP  21_,  1009  (1965). 
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breakdown  waves  is  a complex  subject  and  is  the  subject  of  continuing  ex- 

137) 

perimental  investigation. 

We  present  in  Section  VI.  5.  a a brief  survey  of  high  altitude  aersol 
distributions  and  properties.  In  Section  VI.  5.  b we  look  at  the  microphysics 
of  an  aerosol  in  a laser  beam.  We  study  absorption  due  to  breakdown  of 
the  aerosol  vapor  and  perform  a scaling  with  altitude  of  this  absorption. 

The  effect  of  vapor  absorption  is  found  to  be  negligible.  We  next  look  at 
the  time  scales  required  for  vaporization,  breakdown  and  initiation  of 
laser  supported  absorption  waves.  Breakdown  times  scale  inversely  with 
pressure.  Basing  ourselves  on  breakdown  times  for  the  most  pessimistic 
assumption  where  all  the  laser  heating  goes  into  cascade  ionization  of  the 
electrons,  we  calculate  initiation  times  that  at  high  altitude  are  much 
longer  than  the  laser  pulse  times  considered  for  a pulsed  device. 

We  present  in  Section  V.  6 the  threshold  versus  altitude  for  main  - 
taining  LSD  waves,  once  ignited.  The  calculations  are  based  on  Raizer 
theory  and  closely  model  the  previous  work  done  by  Nielsen  and  Canavan. 
The  threshold  decreases  with  increasing  beam  size  and  decreases  with  in- 
creasing altitude  up  to  a certain  altitude  where  it  reaches  a minimum. 

5 2 

For  a beam  size  of  1 meter  this  minimum  is  less  than  10  W/cm  , and  is 
well  below  the  laser  fluxes  that  we  are  considering.  Thus  the  long  onset 
time  for  breakdown  off  particles  mentioned  above  is  essential  in  order  to 
obtain  pulsed  laser  propagation  at  the  intensities  required. 


(38) 


(37)  Junge,  C.,  Chagnon,  C.,  and  Manson,  J.,  of  Meteor,  1 8,  81,  (1901). 

(38)  Nielsen,  P.  , and  Canavan,  G.  , "Laser  Absorption  Waves  In  The 
Atmosphere"  in  Laser  Interaction  and  Related  Plasma  Phenomena, 

Vol.  1A.  , Schwarz,  H.,  and  Bora,  H.,  editors,  Plenum  Press  (1973). 


a . Aerosol  Concentration  In  The  Atmosphere 

A review  of  the  concentration,  size  distribution  and  composition 


of  aerosols  in  the  atmosphere  up  to  30  km  altitude  can  be  found  in  Junge'  s 


paper  (Ref.  37)  and  is  summarized  below.  Aerosols  below  the  tropopause 
(h  < 10  km)  have  a lifetime  of  a few  days  because  they  are  washed  down  by 
rain.  They  are  composed  of  insoluble  dust  and  various  salts.  Above  the 
tropopause  the  aerosols  are  either  of  atmospheric  origin  or  of  cosmic  origin. 
The  particles  larger  than  1 pm  are  entirely  of  cosmic  origin  because  of 
their  short  residence  time.  Chemical  analysis  show  that  stratospheric 
aerosols  contain  iron,  sulfur  and  silicon. 

The  fall  velocity  vs  altitude  for  particles  of  different  radius  r is 
shown  in  Figure  103,  and  the  residence  time  in  Figure  104.  The  size  dis- 
tribution is  shown  in  Figure  105. 

Figure  105  shows  that  from  r = 0.  1 to  10  pm,  one  has  --  ^ cc  r 

i.e.,  the  number  of  particles  within  a constant  relative  interval  A r/r  varies 

as  r . The  concentration  of  particles  larger  than  0.  1 pm  is  shown  in 

_ 3 

Figure  106  and  is  < 0.1  cm 

Basing  himself  on  the  fall  velocity  of  particles  in  the  stratosphere 

and  on  the  fact  that  the  tropopause  acts  as  a sink  for  particles,  Junge  showed 

that  most  of  the  aerosol  particle  above  the  tropopuase  must  be  of  atmospheric- 

origin.  These  particles,  that  are  now  believed  to  be  solutions  of  H in 

water,  present  no  problem  to  laser  propagation  because  they  will  vaporize 

at  very  low  temperatures  and  can  not  trigger  laser  supported  absorption 

7 2 

waves  at  the  fluxes  which  we  are  considering  (I  < 10  W/cm  ). 

The  particles  that  we  are  concerned  with  are  the  dust  particles  of 
extra-atmospheric  origin  of  size  10  pm  or  larger.  They  are  so  few  in 
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ALTITUOE  (KILOMETERS) 


PARTICLE  CONCENTRATION  AS  PER  CM 


C.  JUNGE  ETAL.  J.  OF  METEOROl 
18.  81  ( 1961 ) 


G57I5  PARTICLE  RADIUS  (/xm) 


Figure  105  Compilation  of  all  Size- Distribution  Data  Determined 
from  Samples  Collected  as  Indicated.  The  numbers 
given  after  the  altitude  refer  to  the  stage  of  the  Ca sella 
impactor  on  which  the  sample  was  collected. 


J UNGE, CHAGNON  AND  MASON 
J.  METEOROL.  18, 81. 1961 


Figure  106 


PARTICLE  CONCENTRATION  (CM3) 


Vertical  Profile  of  Particles  Collected  with  the  General 
Mills  Impactors  on  9 December  1959,  and  the  Available 
Temperature  Soundings  Closest  in  Time  and  Space.  The 
verticle  bars  indicate  the  altitude  interval  over  which  the 
sample  the  sample  was  collected. 


number  that  their  distribution  in  size  and  altitude  is  unknown.  Junge  pre- 
sents data  for  10  /j,m  particle  at  only  two  altitudes  18.4  and  24.5  km.  Since 
these  particles  are  presumed  of  extra  atmospheric  origin,  in  order  to  derive 
the  concentration  vs  altitude,  we  use  the  continuity  equation  V • n u = 0, 

This  tells  us  that  n a l/u  and  since  u is  constant  up  to  35  km  (see  Figure 
103),  we  must  take  n = const,  for  Z < 35  km.  Above  35  km  the  fall  veloc- 
ity increases  since  the  flow  surrounding  a falling  particle  of  radius  10  |jm 

Z /h 

or  less  becomes  free  molecular  and  we  take  n to  vary  as  e ' , where  h 

is  the  scale  height  (h  =*  8 km).  Assuming  that  the  distribution  in  size  of 
large  particles  follows  the  same  law  as  the  small  ones,  we  estimate 
n(r  > 10  fim)  =* 

(z  -35\ 

- 3 Z > 3 5 km 


n(r  > 10  (im)  - 


10-  5 


(Z  -35\ 

8 ). 


1 5 < Z < 35  km 

(412) 


cm 


Basing  ourselves  on  the  previous  considerations,  we  can  expect 

that  the  number  of  particles  greater  than  1 fj.m  in  a beam  of  length  10  km 

2 

between  15  and  25  km  and  cross  section  1 m will  be 

/ A 7 7 

x 10°  cm  x 10^  cm4=  10.  Similarly  for  r > 10  ^m  we 
to  be  10^. 


find  this  number 


Note  added  in  proof.  More  recent  measurements  (Science,  191 , 1270 
(1976))  have  just  come  to  our  attention  which  indicates  that  the  number  of 
particles  larger  than  10  jj.m  is  =*  10"  9 to  10"  10  cm“  i.  e.  , more  than 
four  orders  of  magnitude  below  those  given  by  Eq.  (412)  and  used  in  this 
report.  The  particles  ;n  the  range  1 to  10  (im  would  mainly  be  spherical 
AI2O3  from  rocket  exhaust  and  those  above  10  jim,  would  be  extra  ter- 
restrial origin. 


181 


b.  Breakdown  Induced  by  Aerosols 


( 1 ) Evaporation  and  Breakdown  of  an  Aerosol  Particle 

(39)  (40)  (41  42) 

Analysis  of  experimental  studies  of  Avco  and  elsewhere  ’ 

show  that  particles  that  initiate  air  breakdown  waves  do  so  by  vaporizing  with 

subsequent  breakdown  in  the  vapor  phase  (Al2C>3,C)  and  by  driving  strong 

shocks  in  air  (CaC03,  H2O).  These  conclusions  were  drawn  from  the  fact 

that  the  threshold  depends  on  particle  size,  and  that  particles  that  vaporize 

at  low  temperatures  have  very  high  thresholds  (H2O).  Small  particles  need 

high  fluxes  in  order  to  vaporize  since  heating  is  volumetric,  while  cooling 

by  heat  conduction  with  the  air  varies  as  r,  the  asymptotic  temperature 

2 

reached  varying  therefore  as  r . 

Basing  ourselves  on  sea  level  measurements  of  the  effect  of  aero- 
sols on  air  breakdown,  we  can  derive  some  simple  scaling  relations  with 
altitude. 

Once  a particle  has  vaporized,  the  volume  occupied  by  the  vapor  is 


V * 


4 irr ' 


4 n r 


3 


P sol 

Pair 


T . 
air 


v Q , n4  3-1 
a 8 x 10  r p 
o 


where  rQ  = initial  radius  of  solid  particle,  Tv  = vaporisation  temperature 
and  p = the  pressure  in  atmospheres.  In  the  last  step  we  assumed  Ty  = 
3000°K,  ps  j = 2 and  that  the  molecular  mass  of  the  vapor  was  the  same 
as  air.  The  l/p  dependence  comes  from  the  variation  of  air  density  with 
pressure . 

(39)  "Propagation  Through  Atmospheric  Aerosols,  " Semi-Annual  Technical 
Report,  AERL,  July,  197S. 

(40)  Schlier,  R.  , Pirri,  A.,  and  Reilly,  D.  , "Air  Breakdown  Studies, 

Avco  Report  AFWL-TR-72-74,  February,  1 Q 7 3 . 

(41)  Hall,  B. , Maher,  W.  , and  Wei,  P.  , "An  Investigation  of  Laser 
Supported  Detonation  Waves,"  Boeing  Report  A FW  L-  T R-  7 3-  28 , 

June,  1973. 

(42)  See  Lincoln  Laboratory  Optics  Research  Reports  for  the  period 
1970-1973. 


The  vapor  will  heat  up  due  to  inverse  brems strah lung  absorption, 

but  will  transfer  heat  to  the  surrounding  air  by  heat  conduction.  The  heat 

conductivity  of  air  varies  as  K (T)  =K  Ta  where  a =0.8.  The  heat  lost 

o 

by  a sphere  of  radius  r at  temperature  T due  to  conduction  to  the  surround- 
ing air  is 


4 7r  r KT  -1/3 

6 T = r~  a p ' r 

v L a + 1 K c 

The  absorption  of  laser  energy  on  the  other  hand  is: 


^abs  = a 'ne  nN  1 


4 it  3 


e nN  ^laser  x 3 r 


where  a' = inverse  Br ems strahlung  absorption  factor  and  n^  and  n^  are  the 
concentrations  of  electrons  and  neutrals  respectively.  Assume  T^  independ- 
ent of  altitude  (pressure),  then  nN  <x  p.  Initially,  we  assume  Te  = Tv  (this 
will  not  be  true  at  high  fluxes  or  when  most  of  the  gas  is  dissociated)  and 
that  the  electron  concentrated  is  in  equilibrium  with  T . We  will  have 

1 /?  “fr/^kT'  1 / y 

n « A n.  / e cc  p1'  0 where  e = ionization  potential,  and  A is  a 

e N “ 1 

factor  depending  weakly  on  T.  We  thus  get 


t 3/2-1  T 

*abs  * Ilaser  P P * 1 


1/2  3 


laser  P ro 


2 j / ^ l/2  5 / in  2 

Equating  0abs  with  <f>L>  we  see  that  Ilaser  « rQ"  p"  ' = P"  r0‘ 

l 

Figure  107  presents  curves  of  marginal  stability  for  an  Al^O^  aero- 
sol which  has  a boiling  temperature  of  4000°K  and  an  electron  concentra- 
tion in  the  vapor  phase,  forTv=  4000,K,  Psaj  = 1 atm,  of  8 xlO  part/cc. 
To  the  right  of  each  curve  the  vapor  tends  to  cool  down.  To  the  left  of  each 
curve  the  vapor  will  heat  up. 


PRESSURE  (atm) 


Figure  107  Stability  Boundaries  for  Aerosol  Vapor 
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ALT  I 


The  intersection  of  a constant  I line  with  the  dashed  boundary  at  the 
bottom  of  the  figure  corresponds  to  the  smallest  aerosol  particle  that  will 


vaporize  under  the  given  flux. 

We  can  make  an  estimate  on  the  effect  of  aerosols  on  laser  propaga- 
tion under  the  following  assumptions: 

a)  At  a given  I all  particles  with  r > r given  by  Figure  107 
vaporize  and  the  vapor  expands  and  heats  to  17,000°K  tempera- 

(43) 

ture  of  maximum  absorption  where  the  vapor  is  practically 
totally  ionized. 

b)  We  assume  that  no  detonation  wave  forms  and  that  no  deflagra- 
tion wave  propagates  through  air. 

c)  We  use  the  distribution  of  particle  sizes  of  Figure  105  for 

h > 10  km  with  concentration  given  by  Figure  106  up  to  20  km 
and  by  Eq.  (412)  above.  For  h > 10  km,  we  use  the  distribu- 
tions and  concentrations  under  the  above  assumption  tabulated 
by  McClatchy^^  corresponding  to  a 5 km  haze.  We  find  an 
absorption  coefficient  for  the  vapor 

2 - 1 

k = 0 . 25  x p cm 

and  an  absorption  cross  section  per  particle 

- 7 3 2 

o (r  ) = 6 x 10  r p cm 
o o 

where  r is  expressed  in  nm.  Integrating  over  the  particle 
distribution  of  Figure  106  from  r^  to  10  pm,  we  calculate  the 
absorption  length  shown  in  Figure  108.  This  figure  shows  that 

(4  1)  K a i /. e r , Yu..  Soviet  Physics  ,T?lTP  31,  l 1 48  (1970). 
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ALTITUDE  (Km) 


absorption  by  aerosols  is  only  important  at  sea  level.  If  tin- 
laser  is  situated  at  a mountain  top  one  can  completely  neglect 
the  attenuation,  due  to  the  aerosol  vapor  itself. 

( 2)  Time  Scales  Involved 

The  time  it  takes  an  Al00..,  aerosol  to  heat  up,  melt  and  vaporize 
is  shown  in  Figure  109.  The  relevant  parameter,  I x t,  which  is  propor- 
tional to  the  energy  deposited  in  the  particle,  is  tire  abscissa  and  the  radius 
of  the  particle  is  the  ordinate.  Under  a given  flux  I,  only  particles  above  a 
given  radius  will  either  heat  up  to  T,  melt  or  vaporize  because  of  heat 
losses  to  the  surrounding  air.  This  radius  is  given  by  the  intersect  of  the 

horizontal  dashed  lines  with  the  appropriate  full  curve.  Figure  109  shows 
7 2 

that  for  1=10  W/crn*'  the  minimum  time  for  vaporization  to  occur  is 
0.8  p sec  (for  r = 0.  3 pm).  This  time  becomes  8 psec,  and  80  p sec  for 
1 = 1 o'  and  10^  W/cm°,  respectively.  The  times  for  complete  vaproization 
to  occur  are  twenty  times  longer. 

We  now  look  at  breakdown  in  the  vapor  bubble.  The  electron  distri- 
bution function  (which  determines  T and  ionization  rate)  is  determined  by 
a balance  between  Joule  Heating  and  elastic  and  inelastic  collisional  losses. 

Since  the  elastic  collision  frequency  is  less  than  the  laser  frequency,  both 

(44) 

heating  and  cooling  rates  are  linear  with  pressure  and  the  electron  dis- 
tribution function  will  therefore  depend  only  on  I.  . We  expect  that  the 
vapor  will  break  down  at  a threshold  ltp  independ ent  of  altitude.  The  time 
to  breakdown  should  however  scale  as  l/p. 


(44)  This  statement  is  only  correct  if  we  neglect  the  effect  of  pressure 
on  the  degree  of  dissociation. 


RADIUS  (/im) 


We  derive  an  estimate  of  the  breakdown  time  by  assuming  this  time 

to  be  ten  avalanche  breakdown  times  (i.e.,  enhancement  of  the  electron 
4 

density  by  1 0 ) and  assuming  that  all  the  Joule  heating  goes  into  ionization. 

- 1 5 2 

Using  an  elastic  cross  section  of  10  cm  and  an  ionization  potential  of 
l 3eV  and,  assuming  the  vapor  at  temperature  T’v(0K)  to  be  at  the  ambient 
pressure  p,  we  find 

6 (T  / 4000) 

T „ , , = rT5 — usee  (413) 

Breakdown  n/in's  ^ 1/2  F 

(1/10  ) Te  p 

2 

where  T is  the  electron  temperature  in  eV , I is  in  W/cm  and  p is  in  at- 
mospheres. The  breakdown  time  is  plotted  in  Figure  110.  For  1 = 10^  W/cm 
p = 1,  T ^ = 1 eV,  we  find  that  = 60  psec  and  will  very  probably  be 


longer  than  the  laser  pulse  duration.  Note  that  Eq.  (413)  gives  a lower 
bound  on  the  breakdown  time  since  any  inelastic  process  will  result  in  gas 
heating  and  lengthen  the  breakdown  time. 


PRESSURE 


ALTITUDE  (KM) 


6.  MAINTENANCE  OF  LASER  SUPPORTED  ABSORPTION  WAVES 


Once  the  aerosol  vapor  has  broken  down,  the  breakdown  front  will 
move  into  the  surrounding  air,  initiating  an  air  breakdown  wave.  A good 
understanding  of  the  initiation  mechanism  and  the  dependence  of  threshold 
flux  on  pressure  is  at  present  lacking. 

Subsonic  laser  supported  absorption  waves  ignited  off  solid  targets 

were  maintained  at  fluxes  as  low  as  4 x 10^  W/cm*'  with  a beam  area  of 

2 (45)  -1  6 

0.  5 cm  and  the  threshold  was  found  to  vary  as  r . Studies  of 

ignition  of  laser  supported  detonation  waves  off  aerosols  were  made  at 
AERL,  and  Lincoln  Laboratories^*^  and  give  a threshold  at  sea  level 

of  c*  1.5  x 10^  W/cm*L  Computer  simulations  at  Livermore , ) on  the  other 

hand,  give  a threshold  of  3 x 10°  W/cm  . The  detonation  wave  behavior 
down  to  the  fluxes  quoted  above  were  found  to  agree  well  with  the  Raizer 
detonation  wave  theory.  The  Raizer  theory  assumes  a)  complete  ab- 
sorption of  the  laser  radiation  by  the  air,  b)  that  the  air  can  be  treated  as 
a fluid  with  constant  y (y  - 1.17  for  10,  000  < T < 20,  000°K)  and  c)  that  the 
wave  obeys  the  Chapman  Jouguet  conditions  under  equilibrium  conditions. 

Calculations  of  the  absorption  length  at  the  Chapman  Jouguet  point  show  that 

7 2 

at  sea  level,  for  fluxes  below  10  W/cm  , the  absorption  length  is  larger 

than  1 cm,  i.e.,  much  larger  than  the  beam  diameter  in  all  the  experiments 

^This  low  threshold  may  be  explained  by  the  fact  that  inelastic  collisions 
between  electrons  and  molecules  were  neglected,  resulting  in  electron 
temperature  much  higher  than  heavy  particle  temperature.  At  low  fluxes 
(3x10^  W/cm^)  where  the  gas  is  not  completely  dissociated  departure 
from  LTE  is  not  expected  to  occur. 

(45)  Smith,  D.  C.  and  Fowler,  M.  , Appl.  Phys.  Lett.  22,  500  (1973). 

(46)  Edwards,  A.,  Ferriter,  N.  , Fleck,  J.  and  Winslow,  A.,  "A  Theoretical 
Description  of  the  Interaction  of  a Pulsed  Laser  and  a Target  in  a Laser 
Environment,"  Lawrence  Livermore  Laboratory,  UCRL-51489, 

November  1973. 
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I 

quoted  above.  This  indicated  that  the  threshold  for  maintaining  detonation 

waves  is  determined  by  conditions  in  an  absorption  zone  that  may  be  far 

from  equilibrium.  In  order  to  explain  the  data  and  correctly  predict  the 

threshold  for  maintaining  LSD  waves  one  must  look  at  the  structure  of  the 

detonation  waves.  An  analytic  approach  was  attempted  at  Lawrence  Liver- 

more  Laboratory  by  Ferriter  and  Winslow  and  a computer  model  was 

developed  at  the  Air  Force  Weapons  Laboratory  by  Stamm  and  Nielsen^4^ 

in  order  to  study  the  threshold,  but  none  of  their  results  should  be  given,  in 

our  opinion,  much  credibility,  since  their  treatments  are  too  simplistic. 

The  initiation  threshold  of  detonation  waves  at  10.  6 jim  versus  pres- 

(4  \ ) 

sure  was  measured  by  Hall,  et  al. , who  ignited  the  waves  off  various 

surfaces  in  a chamber,  the  inside  pressure  of  which  could  be  varied.  They 

found  that  the  threshold  either  increased  or  remained  constant  as  one  de-  J 

-2 

creased  the  pressure  from  1 to  10  atm.  In  no  case  were  they  able  to 

6 2 ] 
ignite  a detonation  wave  at  fluxes  below  6x10  W/cm  . At  low  pressures 

the  sharpness  of  the  front  was  lost  and  it  was  difficult  for  the  experimentators  > 

I 

i J 

to  ascertain  whether  they  had  truly  ignited  a detonation  wave.  Their  laser 
pulse  lasted  only  10  /jsec  and  it  is  not  clear  whether  longer  pulses  would 

r 

have  resulted  in  a lower  threshold. 

The  altitude  dependence  of  threshold  for  maintaining  detonation  waves 

(38 ) 

was  studied  by  Nielsen  and  Canavan.  7 They  based  their  calculations  on 

(47)  Ferriter,  N.  and  Winslow,  A.,  "Calculated  Intensity  Threshold  for  the 
Maintenance  of  Laser  Supported  Detonation  Waves  with  Various  Elec- 
tron Densities,"  Lawrence  Livermore  Laboratory,  UCRL-51606. 

(48)  Stamm,  M.  and  Nielsen,  P.  , A1AA  J.  1 3,  205  (1975). 
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Raizer's  Theory  which  takes  into  account  the  lateral  expansion  of  the  gas 
behind  the  front.  If  f is  the  absorption  length  at  the  Chapman  Jouguet  point 
and  R is  the  radius  of  the  beam,  then  the  flux  necessary  to  maintain  the 
wave  following  Raizer  is  increased  by  a factor  (1  + f/R)  over  the  flux  re- 
quired to  maintain  the  wave  if  R were  infinite.  The  threshold  flux  required 
to  maintain  a detonation  wave  at  a given  altitude  is  the  minimum  of  value 
of  I’  defined  by 

I'  = I (1  + f/R) 

where  C is  a function  of  I.  We  extended  the  calculations  of  Nielsen  and 
Canavan,  which  were  done  for  R ^ 1 cm,  to  larger  beam  radii.  The  results 
are  shown  in  Figure  111.  A discussion  of  the  equations  used  to  derive  the 
threshold  flux  is  given  in  the  appendix.  At  low  altitudes,  in  the  linear  por- 
tion of  the  graph,  our  R = 1 cm  curve  is  a factor  of  1.  5 below  the  same 
curve  given  by  Nielsen  and  Canavan,  a consequence  of  the  different  choice 
of  numerical  coefficient  used.  The  important  point  to  be  made  here  is  that 
the  threshold  for  maintaining  detonation  waves  at  high  altitudes  and  large 
beam  radii  is  very  low.  For  example,  with  a beam  radius  of  10  m one  would 
be  able  to  maintain  detonation  waves  at  an  altitude  of  15  km  with  a laser  flux 
as  low  as  4 x 104  w/cm\  Whether  such  detonation  waves  can  be  ignited, 
however,  is  open  to  question.  The  curves  drawn  in  Figure  11  1 were  ter- 
minated at  the  low  density  side,  wherever  f became  equal  to  R. 

The  dashed  lines  in  Figure  111  represent  the  beam  intensity  re- 
quired to  drive  a Chapman  Jouguet  detonation  wave  with  a given  absorption 
length  f . The  derivation  of  these  curves  is  explained  in  the  Appendix. 


One  would  think,  using  Nielsen  and  Canavan's  criterion,  that  by 
going  to  larger  and  larger  beams  one  could  maintain  detonation  waves  at 
arbitrarily  low  fluxes.  Radiation  losses  from  the  laser  heated  region  will 
prevent  this  from  occurring.  If  one  defines  T ^ = f/u  as  the  heating  time 
(where  u =:  D/2  is  the  velocity  of  the  fluid  with  respect  to  the  shock),  one 
must  require  that  the  radiative  decay  time  at  the  Chapman  Jouguet  point  be 
much  longer  than  the  heating  time. 

(4  9) 

We  used  the  values  of  radiation  from  hot  air  calculated  by  Kivel 
to  compare  the  radiative  losses  during  the  heating  time  to  the  energy  con- 
tent of  the  fluid  at  the  Chapman  Jouguet  condition.  Details  of  this  study  are 
given  in  the  Appendix  and  the  results  are  presented  in  Figure  112.  In  this 
figure  we  plotted  the  intensity  required  to  maintain  a detonation  wave  of  a 
given  absorption  length  as  a function  of  density  (altitude).  Radiative  losses 
corresponding  to  20%  and  100%  of  the  energy  absorbed  by  the  gas  are  shown 
(r)  = 0.  2 and  1 curves).  Since  the  fluid  mechanical  losses  at  threshold, 
following  the  criterion  of  Nielson  and  Canavan,  are  roughly  20%  of  the 
absorbed  energy,  radiative  losses  should  dominate  fluid  dynamic  losses  in 
the  region  to  the  left  and  below  the  q = 0.  2 curve.  Radiative  losses  were 
calculated  assuming  the  plasma  is  optically  thin  to  its  own  radiation.  This 
assumption  is  only  true  on  the  portion  of  the  curve  to  the  left  of  the  minimum, 
i.  e.  , at  an  altitude  h > 22  km  for  the  q = 0.  2 curve.  We  expect  radiative 
losses  to  be  important  at  fluxes  below  3 x 10^  W/cm“  and  altitudes  above 
22  km.  A proper  account  of  these  losses,  however,  can  only  be  made  by 
looking  at  the  detailed  structure  of  the  detonation  wave,  which  is  beyond  the 
scope  of  the  present  study. 

(49)  Kivel,  R.  , "Radiation  from  Hot  Air,  " AERL  Research  Report  RR-79, 
October  1939. 
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ALTITUDE  ( Km) 


G7056  DENSITY  P0  (g/cm3) 


Figure  112  Radiative  Losses  in  Detonation  Waves 
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a.  Conclusion 


We  have  found  that  aerosols  in  the  atmosphere  should  have  a negli- 
gible effect  on  laser  beam  propagation  at  fluxes  < 10  ‘ W/cm^,  if  a)  one 
operates  a few  kilometers  above  ground  level  (i.  e.  , on  a mountain)  so  as 
to  avoid  hazy  conditions  and  b)  aerosols  do  not  drive  laser  supported  ab- 
sorption waves  in  the  atmosphere.  Basing  ourselves  on  a very  simple 
model  of  breakdown,  we  find  a delay  for  breakdown  given  by  Eq.  (412)  that 
is  longer  than  the  pulse  durations  that  are  presently  being  considered. 


* 


7.  PROPAGATION  THROUGH  CLOURS 

l'he  prt'St'nct'  of  condensed  \v.\lcr  in  the  path  of  the  laser  will  cause 
an  attenuation  of  the  laser  beanvbv  scattering  and  absorption.  We  estimate 
these  effects  in  this  section.  A review  of  particle  concentrations  and  size 
distributions  in  various  type  clouds  and  toys  was  performed  in  Ref.  ,P>,  and 
we  summarize  the  result  below.  The  size  distribution  is  presented  in 
Figure  11  f.  The  propagation  of  the  laser  beam  is  governed  by  the  equation 


T U . (r.)| 

*— • ,1  ext  i 1 


(4)4) 


where  the  sum  extends  over  all  particles  i of  radius  r.  in  a unit  volume  and 

where  the  extinction  coefficient  at  10.  0 pm  for  water  (sum  of  scattering  4 

absorption)  is  given  by  Figure  114.  One  sees  from  this  figure  that  Q a.  r 

t*  xt 

for  r < 10  pm.  Figure  11  t shows  that  most  particles  in  clouds  and  hazes 
have  a size  v 10  pm.  The  summation  in  Fq.  (4  14)  can  then  be  carried  out 
leading  to  the  simple  relation  for  the  sum  of  the  absorption  ♦ scattering 


coefficient 


a = 70  p 


(4  1*4) 


where  p is  the  mass  density  in  g/cm  of  condensed  water,  (equation  (4  IS) 
is  valid  for  water  at  X 10.  6 pm.  Extinction  at  other  wavelengths  can  be 
obtained  using  Mie  Theory^*'"  and  the  complex  index  of  refraction  of  water*  '' 
or  ice.  ' *"  The  complex  index  of  water  and  ice  over  the  wavelengths  range 
1 to  20  pm  is  shown  in  Figures  1 1 and  116. 


(MH  Van  dor  Hulst,  Tight  Scattering  by  Small  Particles.  John  Wiley  and 
Sons  Publishing.  N.Y.  (19c'7). 

(^  1 ) Centano,  M.  , J.  Opt.  Soc.  Amer.  d 1 . 244  (1°41). 

( “4* J > Oikm.in,  N.  . Advances  in  Physics  7,  P>o  (p)Sg) 


log 


VOLUME  CONCEr 


63931 


Figure  113  MiMsurt'd  Water  Aerosol  Distribution  Functions 
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We  reproduce  in  Table  11  the  data  on  sizes  and  water  content  of 


( 

I 

| 

i 

clouds  and  fogs  that  was  compiled  in  Ref.  39. 

The  meaning  of  R j , R^  and  N in  Table  11  becomes  apparent  when 

. one  refers  to  Figure  117. 

If  one  excludes  thick  fog  and  cumulus  clouds  one  estimates  a water 
T 3 

content  of  10  g/cm  or  less,  which  means  that  a laser  beam  at  10.6  pm 
will  propagate  a distance  of  ^ 1 km  through  clouds  and  fogs  before  suffer- 
ing appreciable  attenuation. 

Cloud  structure  having  a vertical  extent  of  1 km  or  more  will  cause 
significant  attenuation  of  the  laser  beam  and  this  precludes  the  use  of  laser 
propulsion  on  most  cloudy  days.  Attenuation  of  the  beam  will  be  lower  if 
one  operated  in  the  vicinity  of  9 and  4 pm  due  to  the  smaller  absorptance  of 
water  and  ice.  There  is  no  gain  however,  in  operating  much  below  4 pni 
because  of  the  increased  importance  of  scattering  laser  which  scale  as  A”* 
and  will  dominate  over  absorption  losses. 
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TABLE  II.  WATER  AEROSOL  DATA 


Aerosol  Class  N 

Thick  Fog  2 

Thin  Fog  2 

Coast  Haze  2 

Stratus  3 

Strato  Cumulus  3 

Mist  4 

Strato  Cumulus  4 

Fair  Weather  Cumulus  >6 

Cloud  Burst  4 

Cumulus  Congestus  5 

Moderate  Rain  4 


R 1 R 2 P 


A 

(Bm) 

L* 

(Bm) 

(g/cm3) 

10"1 

10 

4 x 10'7 

2 x 10'2 

1 

6 x 10_1  1 

2 x 10_1 

2 

2 x 10'8 

10_1 

20 

4 x 10'7 

5 

20 

io-7 

100 

800 

1 . 6 x 1 0'9 

7 

20 

1 o' 7 

7 

15 

10'7 

800 

3000 

2.  3 x 10'6 

17 

50 

2.  5 x 10' 6 

300 

2 x 103 

6 x 10“8 
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VII.  ECONOMICS 


The  advantage  of  laser  propulsion  as  a means  of  obtaining  direct 
launch  to  geosynchronous  orbit  (GSO)  lies  in  two  facts:  high  volume  through- 
put and  cost  savings.  The  first  is  implicit  in  the  use  of  a high  power  laser 
for  a launch  whose  duration  is  ~ 200  seconds.  A preliminary  outline  of  the 
second  will  be  given  in  this  section. 

Pricing  will  be  in  1976  dollars.  We  note  that  each  item  addressed 
here  will  require  detailed  analysis  in  future  work. 

The  system  cost  required  is  broken  into  the  following  main  components: 


A.  CAPITAL  COST: 

Laser  (closed  cycle,  air  breathing) 
Power  Plant  (peaking  6 GW) 

Beam  control  system 

B.  OPERATING  COST: 


Fuel  cost 
Maintenance 
Propellant 
Vehicle  cost 


1.  CAPITAL  COST 
a.  Laser 

Since  the  major  capital  component  is  likely  to  be  the  laser  array, 
we  will  consider  the  carbon  dioxide  laser  using  the  pulsed  electric  closed- 
cycle  air  breathing  system  to  maximize  efficiency  (over  20%  efficiency  of 
conversion  is  obtained). 
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i 


PKECSDIUO  FiOE  fcu-AAf 


(53  54) 

Systems  studies  * of  laser  scaling  at  AERL  result  in  an  effective 

capital  cost  of  60c  per  watt,  for  pulsed  CO^  lasers.  In  the  present  system  we 

9 

require  a 10  watt  final  beam  aiming  at  the  vehicle.  Allowing  for  atmospheric 
transmission  loss  (20%)  and  pointing  and  tracking  system  loss  (9%),  we  re- 
quire a beam  at  the  laser  exit  of  1374  MW.  Using  a three  stage  MOPA 
system  with  10:1  branching  ratio,  this  power  is  reached  using  1111  1.  4 MW 

units,  which  will  be  stacked  in  a phased  array,  shown  schematically  in 

' 

Figure  118.  The  laser  capital  cost  is  estimated  at  $M933. 

b.  Power  plant 

Estimating  the  electrical  cost  for  a system  of  this  type  may  be  done 
either  by  taking  the  capital  cost  for  the  power  plant  required  and  including 
it  in  the  capital  amortization,  or  by  adding  the  purchased  cost  of  electricity 
(at  4c  per  kwh  to  the  operating  cost.  If  nuclear  or  hydroelectric  power  is 
used  then  amortization  of  the  capital  cost  is  the  major  cost  of  the  power. 

We  therefore  treat  the  power  plants  as  a capital  investment.  Taking  the  f 

capital  cost  of  generating  equipment  at  \l{  per  watt,  a standard  power 

I 

plant^^  capable  of  peaking  at  6 GW  is  estimated  at  $M1020. 

• } 

c.  Beam  control  system  , 

t 

In  the  present  design  a phased  array  of  lasers  is  used  to  provide  a 
correctible  beam.  Conceptually,  phase  corrections  can  be  applied  to  this 

’ ! 

beam  either  before  or  after  the  final  amplifier  stage.  Two  large  mirrors  ! 

are  required:  a chandel,  at  which  the  laser  beams  are  combined  into  a 
single  coherent  phase-corrected  beam,  and  a steerable  beam  director, 
which  aims  the  beam  at  the  target  rocket  engine  (Figure  119).  A pointing 

(53)  Neacy,  W.  , 1976,  private  communication. 

(54)  Feinberg,  R.  and  Roos,  J.  , 1976,  private  communication. 

(55)  Boston  Ediston  Inc.  , 1975,  private  communication. 
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- CHANOEL 
BEAM  COMBINER 
C PHASE  CORRECTOR 


BEAM  DIRECTOR 


Figure  119  Laser  Launching  System 


r 
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and  tracking  system  controls  the  beam  direction.  This  system  is  estimated 
at  $M200. 

2.  OPERATING  COST 

The  costs  will  be  normalized  to  that  required  for  launch  of  1 ton  to 

C.SO. 

a.  Fuel 

The  laser  system  operates  on  the  fuel  JP-4,  the  combustion  products 
of  which  provide  the  lasing  gas  mixture.  For  the  launch  duration  of  256  sec 
required  for  the  standard  case  of  insertion  of  1 ton  into  geo-synchronous 
transfer  ellipse,  we  require  40  tons  of  JP-4,  costing  $4000. 

b.  Propellant 

Possible  propellents  and  costs^^'  are  listed  in  Table  12. 

TABLE  12.  PROPELLANT  COST 


Propellant  (liquid) 

Cost  (1976) 

$ per  ton  liquid 

Ammonia 

198 

Argon 

876 

f 

1 

Hydrogen 

2750 

f 

Methane 

4940 

Nitrogen 

70 

~ 1 

Oxygen 

62 

) 

W ater 

0 

t 

In  the  present  case  we  take  liquid  nitrogen  as  propellant,  and  obtain  the 
cost  for  launch  of  1 ton  to  GSO  as  $150. 


(■>6)  Air  Products  Inc.  , 1975,  private  communication. 
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Maintenance 


The  major  consumable  in  the  capital  equipment  is  the  capacitors 

used  in  the  power  conditioning  circuit  for  the  laser  supply.  These  are 

4 

estimated  to  be  good  for  4x10  launches,  yielding  a cost  of  $2000  per 
launch.  Other  maintenance  is  allowed  at  $500  per  launch, 
d.  Vehicle  Cost 

The  vehicle  is  to  be  kept  as  simple  as  possible,  and  in  one  elementary 
form  consists  of  a reinforced  block  of  ice  surmounted  by  the  payload, 

(Figure  120).  Guidance  is  then  performed  by  the  laser  beam,  the  flux 
density  of  which  is  varied  across  the  rocket  base  to  provide  torques  for 
guidance.  In  another  form  liquid  fuel  is  injected  into  the  engine  cavity 
through  a porous  plate,  and  there  evaporated  and  subsequently  heated.  In 
this  case  also  the  concept  of  steering  the  rocket  engine  by  varying  the  laser 
beam  on  the  ground  may  apply. 

We  note  that  omitting  on-board  guidance  systems  would  render  ex- 
tremely low  rocket  cost  feasible.  However,  a considerable  systems  analysis 
will  be  required  before  a reliable  cost  figure  can  be  estimated  for  the  engine. 
The  uncertainties  of  vehicle  price  suggest  that  it  should  be  left  as  a parameter 
in  the  final  estimate  of  cost  versus  usage. 

3.  SUMMARY 

LASER  PROPULSION  COST  ESTIMATE 
CAPITAL  COST 

1111  1.4  MW  lasers  at  60  c/w  $ 9 32M 

Electron  beam 
Pulse  forming  network 
Cavity  optics 

Gas  generation  system  (closed) 

PM  loss  9% 

Atmospheric  loss  20% 
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1 TON 
PAYLOAD 


5.4  TON 

SOLID  PROPELLANT 


SHOCK  ABSORPTION 
(OPTIONAL) 


NONSYMETRICAL 
OFF  AXIS  LASER 
BEAM  GENERATES 
CONTROL  TORQUE 


THRUST  N 

NORMAL  TO 
BASE  / 

ft 

V 


10  fisec.  2 x 107  WATTS/  cm2 
6 M JOULE,  160 pps 


G 6275 


-ABLATIVE  CASE 
& 1-D  SKIRT  LINER 
0.5  TON  STRUCTURE 


SLIDING  1-D  SKIRT 


Figure  120  LSD  Wave  Rocket  Engine 
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1 PEAKING  POWER  PLANT,  6 GW 


$ 1 02  OM 

1 REAM  CONTROL  SYSTEM  $ 200M 

Phase  corrector 
Beam  combiner 
Beam  director 
Pointing  ki  tracking 

TOTAI^  $1Mi 

OPERATING  COST  f 1 TON  TO  GSO] 


Fuel  cost  2K 

Maintenance  2.  5K 

Propellant  0-  4K 

TOTAL  8 IK 


4.  UTILIZATION  CURVE 

Using  8%  annual  interest  and  20  year  amortization,  the  cost  per 
kilogram  launched  to  GSO  versus  launch  frequency  is  shown  in  Figure  121, 
where  vehicle  cost  is  given  as  parameter.  These  preliminary  estimates 
will  be  refined  when  the  complete  system  is  designed,  and  a systems 
analysis  of  the  laser  propulsion  engine  permits  proper  estimates  of 
vehicle  costs.  At  present  we  note  that  a vehicle  cost  of  $25000  results 
in  net  launch  costs  of  $40/kg.  at  a rate  of  10.  000  metric  tons  per  year. 


214 


Figure  121 


Transportation  Cost  for  1 Ton  for  Geosynchronous 


VIII.  CONCLUSION 


An  investigation  of  the  potential  and  possible  methods  of  laser 
propulsion  has  been  completed.  We  summarize  the  results  obtained  to 


Trajectory  analysis  identifies  the  optimal  mission  as  insertion  of 

a 1 ton  (brennschluss  weight)  vehicle  into  a GEO  transfer  trajectory.  The 

3 

first  10  km  would  be  powered  by  the  ground  laser  station.  Volume  on 

the  order  of  10^  ton/yr  is  feasible. 

The  laser  required,  allowing  for  propagation  losses,  would  be 
o 

-1.4x10  watt.  We  envisage  a stacked  array  of  phase-locked,  repetitively 
pulsed  air-breathing  CO^  EDL's.  These  can  be  arranged  either  as  a phase- 
locked  oscillator  or  an  oscillator  amplifier.  Net  efficiency  anticipated  is 


A beam  combiner  (chandel)  reflects  the  beam  onto  the  beam  director, 

which  is  a 20m  diameter  steerable  mirror.  Laser  location  on  a mountain 

top  (~  10,  000  feet  altitude)  virtually  eliminates  atmospheric  H^O  absorption. 

CC>2  net  absorption  is  reduced  to  - 10%  by  the  bleaching  effect,  and  blooming 

correspondingly  reduced.  Use  of  phase  correcting  mirrors  will  give  a 
3 

2 meter  spot  at  10  km,  as  required  for  the  rocket  engine. 

Examination  of  various  possible  laser  propulsion  rocket  engines 
leads  us  to  select  the  double -pulse  laser  supported  detonation  wave  (LSD) 
engine;  an  initial  laser  pulse  evaporates  a selected  distribution  of  propellant, 
and  a second  pulse  produces  an  LSD  wave  which  heats  the  propellant  uniformly 


PRECEDING  PAGE  ULA1X 


to  a temperature  on  the  order  of  1 eV.  The  propellant  then  expands  out  of 

the  engine.  Engine  efficiencies  near  50%  are  predicted. 

4 

Preliminary  analysis  shows  that  launch  of  ~ 10  tons  per  year  to 
geo-synchronous  (maximum),  at  a cost  of  $20/lb,  should  be  feasible. 

We  conclude  that  laser  propulsion  offers  a promising  alternative 
method  for  economically  conveying  large  amounts  of  material  in  units  of 
-2/3  ton  (allowing  for  chemical  kick-stage)  to  GEO.  Construction  and 
proof-testing  of  the  laser  propulsion  engine  on  the  largest  pulsed  lasers 
available  is  required  as  the  next  step. 
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APPENDIX  A 
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ASTRODYNAMIC  FUNDAMENTALS 


This  Appendix  provides  a very  brief  review  of  elementary  astro- 


dynamics,  together  with  a derivation  of  the  expressions  for  mission 


velocity,  etc.  , which  are  used  in  Section  II  of  this  report.  It  is  based 


primarily  on  class  notes  from  Course  16.41,  a graduate  course  in  the 


Department  of  Aeronautics  at  M.I.T.,  taught  by  one  of  the  authors  of 


this  report  (P.K.  Chapman).  To  facilitate  comparison  with  the  literature, 


conventional  units  are  used  in  this  Appendix,  except  that  specific  results 


used  in  the  next  are  also  given  in  the  units  of  Table  1 . 


Law  of  Newtonian  Gravitation:  Two  bodies  of  masses  m^  , 


located  exert  forces  on  each  other  given  by 


F j = Gm1  £ R =rnll?l 


(A- la) 


— 2 = 1 = m2-2 


(A-lb) 


where  G = 6.67  x 10*  1 ^N-m^/^R^ 


* = *2'  *1 


Then 


* = *2^1 


. Gm]  R'3R  - Gm?R"3R 


(A-  2) 


= ft  R^JL 


wh  ere  p = 


G(m1  + m^l 


The  equation  of  relative  motion  of  the  two  bodies  is  thus 


R + pR'  \il  = 0 (A-  3) 

In  what  follows,  unless  otherwise  stated  one  of  the  bodies  will  be 
the  Earth  and  the  second  body  will  be  a satellite,  of  so  much  smaller 
mass  that  we  may  take 

,,  = Cm  q 

and  regard  R as  the  inertial  rather  than  relative  acceleration  of  the 
satellite  (sine  the  gravitational  attraction  of  the  satellite  does  not  move 
the  Earth  much).  In  this  case,  R is  the  geocentric  position  vector  of  the 
satellite. 

Circular  Orbits:  One  possible  solution  of  the  equation  of  motion 
(A- 3),  and  the  simplest  possible  case,  is  a circular  orbit  about  the  Earth, 
in  which  the  centripetal  acceleration  of  the  satellite  equals  the  gravitational 
acceleration: 

v2/R  = «2R  = p/R2  (A-4) 

where  v is  the  orbital  velocity,  R is  the  (constant)  orbital  radius,  12 
is  the  orbital  angular  velocity.  Thus 

v2  = p/R  (A-  5) 

122  - p/R3  (A- 6) 


From  the  last  equation,  the  periods  of  two  circular  orbits  of 
radii  R and  R,  are  related  by 

1 A* 

Vt2  ■ [Vh] 3/2 

For  very  low  orbit,  with  R = R^  = 6 178  km,  the  radius  of  the  Earth, 


(A-  7) 
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the  acceleration  is  equal  to  the  surface  gravitational  acceleration, 

2 

= 9.81  m/sec  , so  that  Eq.  (A-4)  gives 

v = v (g  R ) = 7.91  km/sec  (A-8) 

o 1 o o 

S = -y  (g^/R^)  = Schuler  frequency  (A-9) 

The  period  of  a very  low  orbit  is  thus  the  Schuler  period,  84.4  min. 

In  general,  the  values  of  parameters  for  other  orbits  may  be  expressed  in 
terms  of  the  low-orbit  values;  e.g.,  the  radius  of  a geosynchronous 
(24-hour)  orbit  is  given  by  Eq.  (A-7)  as 

R = T /T  Z^R  = 24 7°-  2/3R  = 6.63R  = 42270  km  (A-10) 

a so  o 84,4  o o 


Radial  Escape: 
a radius  Rj  to  infinity, 
(A-  3)  gives 


or 


Another  very  simple  case  is  that  of  escape  from 
along  a geocentric  radius.  The  equation  of  motion 

R = inr  (tv2)=  -^/rZ 

(A- 1 1 ) 


2 2,11 

V1  " v2  = Z[i  Rj  ' Rz 


The  minimum  velocity  at  Rj  which  is  required  to  escape  is  given 
by  this  with  v^  = 0,  R^  = 

v2  = 2fi/R, 

= 2v2 


(A- 12) 


from  Eq.  (A- 5).  The  velocity  required  to  escape  from  any  point  is  thus 
\fZ  times  the  circular  velocity  at  that  point.  In  particular,  the  velocity 
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A 


required  to  escape  from  the  surface  of  the  Earth  (neglecting  aerodynamic 
drag,  etc.  ) is 


v - \IZ  v = \l( 2g  R~  11.18  km/ sec 
e o o o 

In  the  units  of  Table  1 , Eq.  (A-  1 2)  reads 

2 _ _2 2 

Vl  rl  1 + *1 

More  General  Orbits:  If  we  take  the  vector  product  of  R with  the 
equation  of  motion  (A- 3)  we  find 

R x (R  + p R”  3 R)  = R x R_  = — (R  x R) 


(A-  I 3) 


(A-  1 2a) 


(A-  14) 


since  R x R = 0 and  R x R = 0.  Thus  the  vector 


h = R x R = constant 


(A- 15) 


When  multiplied  by  the  mass  of  the  satellite,  this  vector  gives  the 
orbital  angular  momentum  of  the  satellite:  Eq.  (A- 15)  is  a statement  of 
angular  momentum  conservation.  Since 

R_  • h = R_.  (R  x R)  = (R  x R)  • R = 0 (A- 16] 

the  motion  takes  place  in  the  plane  normal  to  h. 

If  now  we  take  the  vector  product  of  h with  Eq.  (A- 3),  we  obtain 


R x h = — — (R  x h)  (since  h constant) 


- mr"3  (!i  xR) 

- pR'3|(R  x R)  1 

pR-3  ( R2R  - (R.R)R] 


(A-  17) 
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Now 


A.  __L  _ jL  ,r  j/2 

dt  R ' dt 


= - (R.R)"3//2R.R 


= - r"3r.r 


(A-  1 8) 


so  that  Eq.  (A-17)  may  be  written 


dt  (R  x h)  = R £ + ~ dt  (r)] 


Integration  yields 


= P W <R/R> 


R x h = fj  [ R/R  + ej 


R = 


JlfZlL 


(A-  1 9) 


(A- 20) 


where  £ is  a vector  constant  of  integration.  Taking  the  scalar  product 
of  Ft  with  this  equation  gives 

R.  (R  x h)  = (R  x R)  . h = h2 
= p[R+R.el 
= p R(  1 + e cos  f) 

where  f is  the  angle  between  R and  £.  The  equation  of  the  orbit  is 
thus 


(A- 21 ) 


1 + e cos  f 

This  is  the  polar  equation  of  a conic  section,  with  origin  at  the 
focus.  The  type  of  conic  section  depends  on  the  value  of  the  constant  e, 
as  follows: 


(A-  22) 
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Circle:  e = 0 


Parabola:  e = 1 


Ellipse:  0 < e > 1 Hyperbola:  e > 1 

Properties  of  Elliptic  Orbit:  From  the  definition  (A- 15) 


h = RR  sin  y = R Rx  (A-23) 

where  y is  the  angle  between  the  velocity  and  the  radius  vector  and  R^ 
is  the  component  of  velocity  perpendlar  to  the  radius  vector;  but,  as  the 
radius  vector  sweeps  through  the  angle  f,  the  perpendicular  component 


of  velocity  of  the  tip  is 


Rx  = R f 


(A- 24) 


so  that 


h = Rl  = 2 A 


(A-  25) 


where  A is  the  rate  at  which  the  radius  vector  sweeps  out  area  in  the 
ellipse.  The  area  rate  is  thus  a constant:  this  is  Kepler's  second  law 
of  planetary  motion. 

It  is  clear  from  Eq.  (A- 22)  that  R is  a minimum  when  f = 0,  so  that 
f is  the  geocentric  angle  away  from  the  direction  of  the  point  of  closest 
approach  to  Earth  in  the  orbit,  which  is  called  perigee.  The  point  of  maxi- 
mum distance  from  Earth,  called  apogee,  occurs  when  f = 180°.  From 
Figure  A-l  and  Eq.  (A-22),*  the  sum  of  perigee  and  apogee  radii  ie 


R + R = 2a  = 
p a 


Jl2  /_J_  + _J_\ 

(U  \1+®  1 “ e) 


(A- 26) 


fid  - * ) 
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F2 


C ^F1 


Thus 


h = p a(  I - e ) 


(A- 27) 


R = = a(l  - e) 

p 1 + e 


(A-  28) 


R = h = a(  1+  e) 
a 1 - e 


(A- 29) 


From  the  figure,  the  distance  from  the  center  of  the  ellipse  to 


a focus  is 


c = a - R = ea 
P 


(A- 30) 


When  the  satellite  is  at  the  tip  of  a minor  axis  of  the  ellipse,  the 


geocentric  angle  is  given  by 


cos  f,  = - c/r. 
b b 


from  the  figure,  where  r,  is  the  radius  vector  at  that  time.  Also, 
from  Eq.  (A- 22), 

r = *■%- 

b 1 - ec/r. 


(A- 31) 


2 ..2/  .2  , ..  2.  2 ,2  2 

rb  = (h  /p  + ec ) = [ a ( 1 - e ) + e a ) = a 

2 .2  2 2.2 

= c+b  = ea+b 


(A-  32) 


Thus  the  eccentricity  e is  given  in  terms  of  the  semi-major  and 
semi-minor  axes  of  the  ellipse  by 


= 1 - b"/* 


(A-  33) 


The  area  of  an  ellipse  is  7rab.  If  the  period  of  the  orbit  is  T,  the 


rate  at  which  area  is  swept  out  is 


from  (A- 25).  The  period  of  the  orbit  is  thus  given  by 


(A- 34) 


T = = ZrraZ^l-e2±  = 

' N/l|ia(l-e  ) 


and  depends  only  on  the  semi-major  axis.  Thus  all  orbits  having  the 
same  maximum  diameter  have  the  same  period. 

To  obtain  an  expression  for  velocity  in  the  orbit,  return  to  the 
equation  of  motion  (A- 3)  and  take  the  scalar  product  with  R: 


R-4  = t -ir  <*•*>  = -jr  (tv2) 

= - (iR'3R.R 

= ^ ft  (if) 

the  last  step  being  from  Eq.  (A- 18).  Integration  gives 


where  the  integration  constant  C may  be  obtained  from  conditions  at 
perigee,  where  the  velocity  is  v^  = R^f  (being  perpendicular  to  the 
radius  there).  Thus, 


i % - ^rp 
i >y>2  - n/% 


(A- 3 5) 


(A- 36) 


(A- 37) 
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= fh/Rp)2  - m/Rp  (from  (A- 25) 


. __i± — - (from  A- 27),  (A- 28) 

(1  - e) 


(A- 38) 


_ 2.  .2  a 

2a  (1  - e) 


_ _ _LL 

2a 


so  that,  finally,  Eq.  (A- 37)  reads 


v = 4 


(+-*) 


(A-  39) 


This  very  useful  expression  is  sometimes  known  as  the  vis-viva 
integral.  Note  that  velocity  in  an  orbit  at  a given  radius  depends  only  on 
the  semi-major  axis  of  the  orbit.  In  particular,  if  the  semi-major  axis 
is  infinite,  so  that  the  vehicle  is  on  an  escape  orbit  Eq.  (A- 39)  reduces  to 


v2  = (A-40) 

R 

which  is  the  same  form  found  for  radial  escape  (Eq.  (A- 12).  Thus  escape 
velocity  from  any  point  is  independent  of  the  direction  in  which  the  vehicle 
is  projected. 

Transfer  Ellipse  to  Geosynchronous  Orbit 

If  a transfer  ellipse  with  apogee  at  geosynchronous  altitude  (A- 10' 
is  required,  the  vis  viva  integral  (A- 39)  can  be  used  directly  to  give  the 
required  injection  velocity  only  if  the  perigee  altitude  (and  hence  the  emi- 
major  axis)  is  specified.  However,  from  Eq.  (A- 15),  the  magnitude  of  the 
orbital  angular  momentum  h at  any  point  is  equal  to  the  orbit  rad'us  times 
the  component  of  orbital  velocity  perpendicular  to  the  orbit  radius.  From 
die  geometry  of  Figure  A- 2,  just  at  burnout  we  thus  have 


h = VjRj  sin  ( y - ) 


(A- 4 1 ) 


where  Yj  the  flight  path  angle  at  burnout,  relative  to  the  local  vertical 
at  the  laser  station,  and  is  the  geocentric  angle  traversed  during 
boost.  Clearly 


tan  5j  = pj  sin  Oj/(Ro  pj  cos  Oj ) 


where  R is  the  radius  of  the  Earth.  In  the  units  used  in  Section  II, 
o 


this  gives  Eq.  (108). 
From  Eq.  (A- 29) 


Ra 

e - - 1 

a 


and,  if  we  substitute  this  in  Eq.  (A- 27),  we  find 


(A-  42) 


(A- 4 3 ) 


h = p 2Ra  - 


(A-  44) 


The  semi-major  axis  may  now  be  eliminated  by  means  of  the  vis  viva 
integral.  With  Eq.  (A-41),  we  have 

,2  2 2 2. 

h = v R sin  (y  - 0 ) 


= 2PR,->‘Ra 

2mr»2  (tt  - 


(A-  4 5) 


whence,  finally 


v:  = 2p 


(ri  ■ *. 


1 2 

1 - R-  sin 

a 


,-8.) 


which  gives  the  mission  velocity  Eq.  (107)  for  a transfer  ellipse  to  GSO. 
Circularization  AV 

Finally,  we  calculate  the  AV  required  to  circularize  a coplanar 
transfer  ellipse  to  GSO  at  geosynchronous  altitude.  The  expression  (110) 


for  the  velocity  at  apogee  in 


the  transfer  orbit  is  obtained  directly  by 


noting  that  a - 2(  Rp  RJ.  so  that  Eq.  (A-  39)  gives 


f—  - — 1 

\Ra  % 


2|lR, 


R 


R (R  Ra 

a p a 


) 


(A-  46) 


For  the  applications  in  Section  II,  where  the  perigee  radius  is  not 
directly  known,  a more  useful  expression  is  found  by  equating  the  orbital 
angular  momentum  to  that  at  burnout,  given  by  Eq.  (A-41) 

h = v R = v R.  sin  (y  . - 6)  (A-47) 

a a 11  1 


so 


R, 


v = v, 
a 


“1  , a 1 (A- 48) 

'1  — sin  S - V 

The  circular  velocity  at  synchronous  altitude  is  given  by  (A- 5)  as 


2_ 

R 


so  that  the  AV  required  of  the  circularization  kick  stage  is 


AV, 


■k  - 3i"  (V| ' 

which  is  equivalent  to  Eq.  (172) 

Notice  that,  if  Eq.  (A-48)  is  used  in  Eq.  (A-45),  one  obtains 


i-  v2 

2 1 


_JL.  Iv2 


R 


2 a 


-Ji- 

lt 

a 


(A-  49) 


(A-  SO) 


(A-Sl) 


wMch  states  that  the  sum  of  kinetic  an, I Bravitational  potential  enerpy  of  the 
vehicle  at  injection  to  the  tcanstec  ellipse  is  equal  to  that  at  apopoc.  Thus. 
Kq.  (A-451  could  have  been  obtained  simply  from  enerpy  conservation. 
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APPENDIX  B 


PRELIMINARY  SYSTEM  CONCEPT 


1.  AEROSOL  INJECTION 

The  detonation -wave  laser  propulsion  (LP)  engine  requires  a uniform 
distribution  of  propellant  gas  in  the  chamber  before  the  laser  supported  de- 
tonation wave  travels  through  it.  This  distribution  may  be  prepared  by  ex- 
pansion of  gas  evaporated  from  a condensed  phase  using  a preliminary  laser 
pulse;  it  may  also  be  prepared  by  electric-arc  induced  heating  of  a condensed 
propellant  layer.  In  this  Appendix  we  discuss  a third  method,  the  generation 
of  an  aerosol  whose  mean  density  equals  that  required  in  the  propellant. 

If  an  aerosol  is  generated  and  ejected  into  the  vacuum  from  the  en- 
gine feedplate,  then  provided  the  particles  do  not  evaporate  too  rapidly  in 
the  interpulse  period,  when  the  main  laser  pulse  arrives  the  aerosol  will  be 
evaporated  extremely  rapidly.  The  released  gas,  expanding  in  all  directions 
from  each  droplet,  will  provide  a uniform  gaseous  propellant.  The  laser 
pulse,  having  propagated  (by  evaporating  the  droplets)  through  the  aerosol 
to  the  end  plate,  will  then  strike  a detonation  wave  which  will  travel  back 
through  the  propellant  gas,  thus  heating  it  to  the  required  levels.  We  will 
examine  aerosol  evaporation,  generation  and  distribution. 

2.  AEROSOL  EVAPORATION 

The  behaviour  of  a droplet  in  a vacuum  is  predicted  from  first  prin- 
ciples by  calculating  the  energy  loss  rate  due  to  evaporation  and  allowing 
for  the  effect  of  the  consequent  temperature  decrease  on  the  evaporation 
rate,  using  the  Clapeyron  equation.  The  rate  of  change  of  radius  r will  be 


PHKCKD1N0  FAC* 


dr 

dt 


(B-l ) 


n c M 
4 Pn  o 

where  c is  local  soundspeed,  M molecular  weight,  p droplet  density  and 
n^  Avogadro's  number.  The  local  gas  density  n near  the  droplet  may  be 
taken  as  that  causing  the  vapour  pressure  (obtained  from  Clapeyron's  equa- 
tion), so  that: 


-LM/RT 

nTT 


(B-2) 


w 


o Q 

here  n'  is  Loschmidt's  number,  T'  = 273  K,  pQ  = 5.67  x 10  for  H^O. 


and  L is  the  droplet  vaporization  latent  heat.  Since  the  heat  loss  rate  will 
be 


do  4 3 _ dT  . . 2 dr 
dt  = 3 *r  P cn  dT  = L*4*p  r hT 


p dt 


dt 


( B -2a) 


where  is  the  droplet  specific  heat,  we  obtain 


T-  & « 

P 


(B-3) 


where  and  r^  are  the  initial  values  of  temperature  and  radius,  com- 
bining Eqs.  (B-l),  (B-2)  and  (B-2a)  we  obtain  the  differential  equation  for 
the  droplet  temperature  as: 


dT 

dt 


3AL 
r C T 
1 P 


C 

(Tj  - T)  - LM/RT 


(B-4) 


9 o 

where  A is  a numerical  constant;  A = 1 . 25  x 10  cm  K/sec  for  H^O.  Sup- 
posing that  T is  limited  to  (relatively)  small  excursions,  Eq.  (B-4)  can  then 
be  integrated  and  yields  an  expression  for  droplet  radius: 
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(B-5) 


The  fractional  mass  loss  is  then 


- fc)’  ' 


(B-5a) 


We  note  that  if  the  droplet  undergoes  a phase  change  and  freezes,  which  will 
happen  if  a suitable  crystallizing  nucleus  is  present,  since  the  droplet  will 
become  supercooled  very  soon  after  its  exposure  to  the  vacuum,  then  the 
fractional  radius  change  required  to  induce  freezing  will  be: 


M"1# )' 


(B-6) 


where  L'  is  the  latent  heat  and  p'  the  density  of  ice.  Consequently  ~ 5%  of 
the  radius  is  lost  in  the  phase  change  for  the  freezing  of  11.,  O droplets. 

Combining  this  with  Eqs.  (B-5)  and  (B-5a)  gives  the  fractional  mass 
loss  as  function  of  initial  radius  for  various  times  after  exposure  to  the 
vacuum.  This  is  shown  in  Figure  B-l,  where  the  values  of  Am/m  are 
shown  for  both  frozen  and  supercooled  droplets. 

For  a typical  pulsed  laser  engine  design,  the  interpulse  period  would 
be  10  msec.  It  is  probable  that  droplets  remain  supercooled  for  this  period. 
We  observe  that  for  droplets  larger  than  20  pm  in  diameter,  less  than  1 0°'<» 
of  the  mass  would  be  lost  to  evaporation  in  this  period.  Thus  if  an  aerosol 
of  this  type  could  be  created  in  vacuo,  it  would  exist  for  long  enough  to  be 
vi  s cd . 
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3.  AEROSOL  PRODUCTION  AND  DISTRIBUTION 

Whereas  the  stability  of  an  aerosol  in  vacuo  appears  to  be  satisfac- 
tory for  our  purposes,  its  production  and  distribution  is  more  problemati- 
cal. The  main  difficulty  is  that  although  aerosols  can  be  produced  in  air  by 
inducing  a Rayleigh  instability  at  the  interface  of  a liquid  stream  or  film  and 
a gas,  in  a vacuum  this  is  impossible.  We  are  therefore  forced  to  use  the 

Rayleigh  jet  instability.  We  note  that  alternate  methods,  such  as  genera- 

(1  2) 

tion  of  droplets  by  high  electric  fields,  while  theoretically  possible,'  ’ 

do  not  permit  a large  decrease  in  droplet  radius  over  that  produced  by  the 

jet  instability.  A useful  rule  of  thumb  concerning  the  latter  is  that  the  mean 

( 3) 

droplet  diameter  is  roughly  twice  the  jet  diameter.  Since  we  require 
droplets  with  diameters  near  the  laser  wavelength  of  10  pm,  at  which  Mie 
absorption  (which  corresponds  to  the  maximum  absorption  efficiency  per 
unit  mass)  will  become  important,  this  implies  that  we  require  an  aerosol 
plate  with  a very  large  number  of  very  small  holes,  through  which  the  liquid 
is  forced  by  ahigh  pressure  difference.  The  velocity  of  the  emergent  drop- 
lets will  then  be  u = \l  p/ p , where  p is  the  stagnation  pressure.  For  a 
continuously  operating  droplet  source,  we  require  that  the  droplets  travel 
a distance  during  the  interpulse  phase  equal  to  that  traversed  by  the  deton- 
ation wave  during  the  main  laser  pulse.  Thus  the  detonation  chamber  length 


d = u r.  = V T 
i P 

( 1)  Peskin,  R.  L.  and  Raco,  R.  J.  , Drop  Size  from  a Liquid  Jet  in  a Longitudinal 
Electric  Field,  Technical  notes,  AIAA  Journal,  781,  (1964). 

(2^  Lapple,  C.  E.  , Electrostatic  Phenomena  with  Particulates,  Adv.  in  Chem. 
Eng.  8,  1 (1970),  (A.  P.  New  York). 

( 3)  Hidy,  G.  M.  and  Brock,  J.  R.  , The  Dynamics  of  Ae rocolloidal  Systems, 
Pergamon,  1970. 
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where  j.  and  j are  the  interpulse  and  laser  pulse  durations,  and  V is  the 
1 p 

detonation  wave  velocity.  The  latter  is  given  by 


- [tv 


where 


1 


for  H^O  vapor,  pc  is  the  aerosol  mass  density,  and 
Now  the  aerosol  density  must  be 


Pc 


rh  t. 
1 


u Ti  A 


c 


is  laser  flux  density. 


where  the  mass  flow  rh  = A.up;  A^  is  chamber  area  and  Aj  is  combined 
droplet  jet  area.  Consequently 

Pc  = ^ • P ■ (B-7) 

The  density  needed  in  the  detonation  chamber  must  be  at  least  enough  to 
provide  a Mie  absorption  length  less  than  the  chamber  length.  The  Mie 
absorption  length^  for  ICO  at  10.6  pm  may  be  written 


2 

L 


4.44  x 10 


Mie 


cm. 


Since  we  require 


Lw.  < d,  then 
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4 . 4 x 10 
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(4)  Reilly,  J.  P.  et  al.  , Multiple  Pul se  Propagation  in  Fog,  Rain  and 
Dust  at  10.  6 pm.  Interim  Report,  30  Sept.  197S  to  1 April  1976. 
Contract  No.  N00 1 73-76-C-0059  (ONR). 
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S 2.  36  x 10  ^ g cm  \ for  r ? w 20  psec  and  $ ~ 10^  W/cm^ 


Then  the  open  fraction  of  the  aerosol  plate  must  be 


A. 


> 2.  36  x 1 0 


-4 


Thus  considering  10  pm  diameter  jet  holes,  which  will  produce  drop- 
lets at  the  upper  limit  of  acceptable  diameter,  we  note  using  Eq.  (B-7)  that 
we  will  require  more  than  300  such  holes  per  square  centimeter  of  detona- 
tion chamber  aerosol  plate.  This  is  a considerable  engineering  problem. 

Lot  us  ignore  for  a moment  the  exact  mechanism  of  droplet  forma- 
tion. It  should  be  noted  that  large  droplets  will  be  shattered  very  rapidly 

(41 

by  the  laser  pulse,  thus  producing  smaller  particles  whose  Mie  absorp- 
tion efficiency  will  be  much  higher.  For  fluxes  </>  > 1 0 , explosive  vapor- 
ization and  shattering  will  occur  in  particles  with  radius  r > 10  pm,  and 

(4) 

in  a time  scale  on  the  order  of  the  sound  travel  time  across  the  droplet 
(about  1 p sec  for  r > 10  pm  droplets).  It  is  thus  imaginable  that  we  distrib- 
ute larger  aerosol  particles  and  allow  the  natural  shattering  process  to  pro- 
vide us  with  an  efficiently  absorbing  aerosol  of  Mie  droplets.  Then  the  number 
of  aerosol-generating  holes  is  relaxed.  We  take  the  above  numbers  for  the 

case  of  a large  pulsed  laser,  and  obtain  r = 20  psec,  r*  = 4 x 10  ' sec. 

7 2 

Then  taking  the  flux  as  10  W/cm  , lower  limit  of  surface  breakdown,  we 

1/3  -4 

have  d = Vr^,  V = Vo  (<,3/pc)  ' ; taking  p . - 2.  36  x 10  , as  above,  then 
V = 9.  31  x 10  % IQ6  cm/sec,  and  d = 18.6  cm.  Then  the  velocity  u is 


4 3 


u = -£  v 

T. 

1 


5 x 10' 


which  implies  a stagnation  pressure  P = Pu  = 25  atm,  or  355  psi.  The  mass 

flow  m = A.  up  is  then  1.  18  g/sec.  For  a 500  kW  laser  the  peak  power  i6 

8 2 
10  W.  Area  of  the  aerosol  plate  should  therefore  be  ~10  crn  in  order  to 

get  10 ' W/cm^,  so  the  total  mass  flow  will  be  11.  8 g/ sec,  and  the  enthalpy 

4 

per  gram  will  be  4.  24  x 10  J. 

The  stagnation  pressure  (355  psi)  is  lower  by  a factor  of  5 than  the 
pressure  usually  used  for  producing  unstable  jets  that  generate  fine  aerosols. 
We  note  that  a further  possibility  would  be  to  produce  a sufficiently  fine  aero- 
sol with  the  conventional  stagnation  pressure  (~  1500  psi).  The  droplets  would, 
in  this  case,  be  traveling  away  from  the  plate  at  more  than  twice  the  required 
velocity.  They  could  be  decelerated  by  changing  the  droplets  and  using  prop- 
erly arranged  electric  fields.  This  arrangement,  however,  would  greatly 
improve  the  complexity  of  the  laser  feed.  It  appears  preferable  to  find  ways 
of  producing  aerosols  (e.g.  , by  colliding  jets)  using  the  stagnation  pressure 
above  (355  psi).  Experimental  work  is  needed  to  determine: 

1)  Whether  adequately  fine  aerosols  can  be  produced  at  this  stag- 
nation pressure. 

2)  Whether  the  laser  beam  will  efficiently  shatter  oversized  drop- 
lets, thus  producing  a medium  which  will  absorb  it  efficiently. 

While  it  is  possible  that  electrostatic  aerosol  production  or  ultra- 
sonic jet  disintegration  could  be  used  to  produce  the  small  droplets  needed, 
it  appears  at  present  as  though  alternate  methods  of  producing  the  required 
propellant  distribution  will  be  easier  to  handle  and  cheaper  to  produce. 
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APPENDIX  C 


THERMODYNAMIC  PROPERTIES  OF  SELECTED  PROPELLANTS 


The  available  thermal  energy  for  conversion  to  kinetic  energy  and 
ultimately  to  specific  impulse  is  most  conveniently  observed  for  a specific 
gas  on  a Mollier  Chart.  This  chart  provides  the  functional  form  of  h (p,  T, 

<£  ) and  s (p,T,$)  where  h is  the  enthalpy,  s is  the  entropy,  p is  the  pres- 
sure, T is  the  temperature,  and  <£  is  a measure  of  the  "chemical"  state 
of  the  gas,  i.  e.  , whether  it  is  dissociated  or  ionized,  etc.  For  a simple 
system  thermodynamic  variables  (h,  s)  are  a function  of  only  two  variables 
p and  T.  At  elevated  temperatures  (j>  may  be  considered  as  a compressi- 
bility factor.  Steady  propulsion  system  are  most  conveniently  described  by 
pressure  and  temperature.  Unsteady  (or  pulsed)  systems  naturally  favor  the 
use  of  two  different  thermodynamic  variables  - the  specific  internal  energy 
and  the  density. 


1 ' H2° 


Water  is  one  of  our  most  common  chemical  compounds.  We  recog- 
nize water  in  its  solid  state  as  ice,  in  its  liquid  state,  and  of  course  in  its 
gaseous  state  as  steam.  Much  is  known  about  water  at  these  low  tempera- 
tures (<  1 000°K)  since  it  has  played  such  a big  part  in  steam  generators  as 
well  as  many  other  technological  applications.  For  our  laser  propulsion 
applications  we  require  high-temperature  properties  of  H^O  at  sub-  and 
super-atmospheric  pressure  levels. 

In  order  to  underline  the  fact  that  H^O  exhibits  rather  unique  behavior 
we  first  describe  the  variation  of  the  ratio  of  specific  heats  y with  T (°K)  and 
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p (atm).  At  very  low  temperatures  y = 4/3  which  is  typical  for  a polyatomic 
gas,  however  as  the  temperature  increases  there  is  very  little  range  where 
y may  be  approximated  to  be  a constant. 

Figure  C-l  shows  y = y (p,T)  for  both  equilibrium  and  frozen  chem- 
istry as  computed  at  AERL  using  an  equilibrium  chemistry  code.  The  in- 
crease of  y between  20Q0°K  and  3000°K  is  a result  of  the  dissociation  of 
H-,0  since  dissociation  essentially  makes  the  gas  more  diatomic  - or  at  very 
high  temperatures  more  monatomic.  This  behavior  can  be  seen  in  Figures 
C-2  to  C-4  w'here  we  display  the  composition  of  H^O  at  1,  10,  and  50  atm. 

The  Mollier  Chart  for  H^O  is  given  in  Figure  C-5  where  we  have  the 
pressure  and  temperature  on  the  enthalpy-entropy  plane.  Note  that  the 
dashed  lines  correspond  to  the  compressibility  factor 

Z = 0.2196  (C-l) 

p (gm/cm'  ) T (°K) 

from  Z,  p,  and  T we  can  compute  p for  any  position  on  the  chart. 

2.  OXYGEN,  02 

In  Figure  C-6  we  show  the  Mollier  Chart  for  oxygen.  We  note  immed- 
iately that  in  the  temperature  and  pressure  range  of  the  chart  (same  as  that 
for  H^O)  there  is  less  available  enthalpy.  Furthermore,  at  first  glance  the 
efficiencies  of  constant  entropy  expansion  processes  appear  to  be  inferior 
to  H^O.  In  general,  the  shape  of  the  02  chart  is  the  same  as  that  of  H^O. 

This  statement  can  not  be  made  for  N^. 

3.  NITROGEN,  N2 

In  Figure  C-7  we  display  the  nitrogen  Mollier  Chart.  The  general 
shape  of  this  chart  is  different  from  that  of  H20  and  02 . Note  the  region 
at  high  pressures  which  is  not  very  influenced  by  dissociation  effects.  This 
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Figure  C-3  H.,0  Composition  at  10  atm 
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Figure  C-5  Mollier  Chart  for  11,0  (p,  T,  7) 
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Figure'  C-6  Mollier  Chart  for  O. 
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Figure  C-7  Mollier  Chart  for 


feature  implies  that  one  can  recover  this  energy  via  constant  entropy  ex- 
pansion and  obtain  efficient  cycles.  Extrapolation  of  conditions  seems  to 
indicate  that  operation  at  8000  - 9000°K  under  100  atm  would  be  rather 
favorable  - leading  to  a specific  impulse  of  about  750  seconds. 

4.  ARGON,  Ar 

In  Figure  C-8  we  see  a Mollier  Chart  for  Argon.  This  chart  has 
a much  simplier  structure  than  the  previous  three  charts.  This  is  pri- 
marily due  to  the  lack  of  real  gas  phenomena  - i.e.  , ionization  of  Ar.  The 
chart  also  implies  that  we  would  require  ve ry  high  temperature  levels  to 
obtain  laser-propulsion  enthalpy  increases. 

5.  HYDROGEN,  H2 

In  Figure  C-9  we  see  the  Mollier  Chart  for  hydrogen.  This  curve 
was  generated  by  NASA  (SP-3002,  1964).  We  have  outlined  the  regions  of 
calculation  for  H-,0,  O^,  N^,  Ar  which  were  generated  at  AERL  using  avail- 
able JANAF  tapes  to  temperatures  of  6000°K.  In  order  to  provide  a point 
of  comparison  for  available  enthalpy,  we  cite  the  Ah  differential  between 
100G°K  about  6000°K.  This  enthalpy  difference  is  1.  5 x 10^  j/gm  which 
exceeds  the  differences  over  the  same  temperature  range  for  H^O,  02>  N 
and  Ar. 
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Figure  C-8  Mollier  Chart  for  Ar 
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APPENDIX  D 


CALCULATION  OF  THE  MAINTENANCE  THRESHOLD  OF 
DETONATION  WAVES  IN  AIR 


We  start  from  the  following  equations: 

1 A pZxZ 

1 " -pTT 


B 


t3/2  e-0/T 


1 - x 


P 


I = Pp  e3/2  (1  + C/R) 


e = C 


(vi  c / v 0.  12 

■ T.  V 

i o, ooo j yp  f 


(D-l) 

(D-2) 
(D-3 ) 

( D-4 ) 


where  f = absorption  length  in  the  plasma  (cm),  R = radius  of  the  beam,  p, 
T and  x are  the  density  (g/cm\  temperature  (°K)  energy  per  unit  mass 
(erg/g),  and  fractional  ionization  at  the  Chapman  Jouguet  point.  Equation 
(D-l)  gives  the  absorption  length  due  to  electron  inverse  bremmstrahlung. 
At  A = 10.  6 pm  the  coefficient  A is: 

3/2 


A = 2 x 1 


(D-5) 


.-7 


Equation  (D-2)  is  an  approximate  Saha  equation  for  air  with  B = 10  and 
P = 1 60,  000° K (13.6  eV).  Equation  (D-3)  can  be  derived  from  Raizer'  s 
theory*1'  and  relates  conditions  at  the  Chapman  Jouguet  point  to  the  laser 
flux.  The  coefficient  (3  ~ (y  - 1 ) ^2  (y  + 1 )/(2y ' ^2 ) ^ 0.  4,  where  y is  the 


(It  Raizer,  Yu.,  Soviet  Physics  JETP  21 , 1009  (1965). 
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effective  coefficient  of  the  adiabat  (y  — 1.  17).  Equation  (D-4)  is  the  energy 
equation  of  state  for  air'*-'  with  c = 2.  8 x 10**  in  cgs  units  and  with  p g the 
standard  density  (1.  3x10'  g/cm  ).  Note  that  Eq.  (D-4)  implies  y = 1.  24, 
so  that  our  choice  of  y = 1.  17  in  the  evaluation  of  (3  is  not  completely  con- 
sistent with  Eq.  (D-4). 

System  (D-l)  through  (D-4)  represents  4 equations  with  6 unknowns, 

p,  p,  T,  X,  e and  I.  If  we  fix  p we  can  eliminate  x,  T and  e and  evaluate  I 

as  a function  of  f.  We  then  seek  the  minimum  of  I when  p is  varied.  This 

(3) 

is  essentially  the  procedure  used  by  Nielson  and  Canavan  in  order  to  de- 
termine the  maintenance  threshold.  The  resulting  threshold  is  plotted  in 
Eigure  111.  At  low  altitude  and  long  absorption  lengths  inverse  bremm- 
strahlung  absorption  due  to  electron  collisions  with  neutrals  rather  than  ions 
is  predominant  and  we  took  these  into  account  by  modifying  the  coefficient  A. 

'VZ  * 0.571  + ‘-2*  ‘°3 * *  T2-8 


= 2 x 10lfc  2 t 0.  57_ 


The  last  term  is  due  to  electron  collisions  with  O and  N atoms.  For  T = 

10,  000°K  both  terms  are  of  same  magnitude  when  the  degree  of  ionization 
is  1%.  The  inclusion  of  electron  neutral  collisions  has  the  effect  of  reducing 
the  threshold  flux  by  a factor  of  1 . 5 at  sea  level  for  R = 1,  000  cm.  The 
effect  at  smaller  beam  radii  and  higher  altitudes  was  found  to  be  neglibible. 

Simple  analytic  results  can  be  derived  when  the  degree  of  ionization 
is  small,  i.  e.  , x«l.  Eliminating  x between  Eqs.  (D- 1 ) and  (D-2),  we  get 

d 


T = 


Pn  (AB  p p) 


(D-6) 


(2)  Zel'dovich,  Ya.  and  Raizer,  Yu.,  Vol.  I of  Physics  of  Shock  Waves  and 
High  Temperature  Hydrodynamic  Phenomena,  Academic  Press,  N.  Y. 

(1«66),  1,  209. 

(3)  Nielson,  P.  and  Canavan,  G.  , "Laser  Absorption  Waves  in  the  Atmosphere' 

in  Laser  Interaction  and  Related  Plasma  Phenomena,  Vol.  3A  (H.  Hora  and 

H.  Schwarz.,  Editors)  Plenum  Press.  N.  Y.  (1973),  p.  177. 
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. 


Inserting  this  value  for  T into  Eqs.  ( D — 3 ) and  (D-4)  and  eliminating  e be- 


tween the  two  we  get: 


1 = p e3^2  p 


o-  82( « V/4  (,  + i) 

\fn  AB  pP  / V R/ 


( D-7 ) 


Taking  the  derivative  of  (D-7)  with  respect  to  P,  neglecting  the  variation  of 
A with  T (or  P),  we  get  the  implicit  equation  from  which  p can  be  calculated 
in  order  to  find  the  threshold  flux. 


e - — ^ / D 

* " 4 Pn  (AB  pP  ) - 9/4 

The  solution  to  Eq.  (D-8)  is  shown  in  Figure  D-l  as  a function  of  p for 
various  values  of  R.  One  sees  from  this  figure  that  p/R  varies  in  the 
narrow  range  0.  13  to  0.  36  when  p varies  over  4 orders  of  magnitude  and 
L varies  for  1 cm  to  100  cm.  The  straight  portion  of  the  threshold  curves 
in  Figure  1 1 1 is  therefore  well  approximated  by  setting  p/R  =■  0.  2 in  Eq. 
(D-3)  and  the  slope  of  the  curves  I(p)  is  0.82.  The  separation  between 
the  curves  corresponding  to  different  beam  radii  is  due  to  the  factor 


(D-8) 


[ Pn  (AB  pP)) 


"9//4  ~ f pn 


[(n  (0.2ABpR)| 


It  is  of  interest  to  compute  the  laser  flux  at  a given  altitude  required 
to  drive  a detonation  wave  with  a given  absorption  length  p.  In  order  to  do 
this,  one  solves  the  set  of  Eqs.  (D-l)  through  (D-4)  with  Eq.  (D-3)  replaced 


I = Pp  e' 


(D-9) 


Curves  l(p)  corresponding  to  constant  values  of  f are  plotted  in  the  overlay 
to  Figure  111. 

An  analytic  representation  for  the  straight  portion  of  the  constant  f 
curves,  when  x « 1,  is  given  by: 
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I = 1.41x10 


<P<A.) 


1 + 6.5x10  (n  f 


where  p^  standard  density  and  p^  = upstream  density  (=  y p). 

The  constant  f curves  have  a minimum,  and  terminate  on  the  low 
density  side  with  a vertical  slope.  The  reason  the  curves  terminate  is  that 
the  solution  to  Eqs.  (D-l)  and  (D-2)  is  double  valued  in  T,  the  lower  value 
of  T being  the  stable  solution.  The  low  altitude  limit  corresponds  to  the 
convergence  of  the  two  solutions  given  by  the  condition  (dp/d  T) ^ = 0 or 
(dp/d  l)  = 0,  i.e.  , the  curves  terminate  with  vertical  slope  in  the  overlay. 
The  minimum  density  or  maximum  altitude  at  which  a detonation  wave  with 
a given  absorption  length  can  be  maintained  is  obtained  by  eliminating  x be- 
tween Eqs.  (D-l)  and  (D-2),  which  leads  to 

c9/T  rj,  3/4 

" ' ^ ~ jn 

and  writing  that  (dp/d  T)^  = 0.  This  leads  to  the  additional  equation 


6/7  2 t3/4 

e “4 


(D-l 0) 


Combining  Eq.  (D-10)  with  Eqs.  (D-l)  and  (D-2)  one  can  solve  from  T,  x 
and  p at  the  maximum  altitude  as  a function  of  (. 

The  maintenance  threshold  that  we  derived  and  that  is  shown  in 
Figure  l 11  does  not  include  the  effect  of  radiative  losses.  We  now  study 
these  losses  in  order  to  ascertain  how  important  these  may  be  in  deter- 
mining the  threshold.  Figure  D-2,  taken  from  Ref.  4 shows  the  energy 
radiated,  u),  per  unit  volume  of  hot  air  at  density  p and  temperature  T.  In 


(41  Kivel,  B.  , "Radiation  from  Hot  Air  and  Stagnation  Heating,"  RR-79, 
Avco  Everett  Research  Laboratory,  Inc.,  October  H>59. 


ENERGY  RADIATED  PER  VOLUME  (watts/cm 


H4362 


AVCO  SHOCK  TUBE 
MEASUREMENTS  - 


P/PQ  = 10. 


/ i 

r Jr- 


I d_5^ — 


//  / 
//  y 
/ / 
/ / 


P • 1.23  x 10  3gm/cm3 


2 4 6 8 10  12  14  16 

TEMPERATURE  (T/I000°K) 


Figure  D-2  Total  Emitted  Radiation  Energy  per  Unit  Volume  of  High 
Temperature  Air  in  Full  Equilibrium  as  a Function  of 
Temperature  for  Constant  Values  of  the  Density 
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the  temperature  range  8-12,  000°K  and  density  range  (p/p  ) = 1 to  10"^  we 


approximate  the  curves  by 


W = 10' 


(pj)(io.T°oo) 


erg/(cm  sec) 


(D-ll) 


where  p^  is  the  standard  density. 

We  require,  for  radiative  losses  to  be  negligible,  that  the  amount 
radiated  per  unit  mass  during  the  residence  time  of  the  flux  in  the  absorption 
zone  be  much  smaller  than  the  energy  per  unit  mass  stared  in  kinetic,  disso- 
ciation and  ionization  energy.  This  requirement  can  be  written  as 


= T)  « 1 


(D-12) 


where  e is  given  by  Eq.  (D-4),  (f/u)  is  the  residence  time,  and  the  velocity 
u of  the  fluid  with  respect  to  the  shock  is  given  by,  ^ ' 


(0-13) 


with  use  of  Eqs.  (D-4),  (D-6),  (D-7),  (D-ll)  and  ( D — 1 3 ) , in  Eq.  (D-12) 
becomes 


7.  75 

((cm)  « 52  j^l  + 6.  5 x 10"2  (n^ j p 


(D-14) 


We  have  plotted  in  Figure  D-3  the  absorption  length  f vs  normalized  density 
p /p  for  which  a constant  fraction  q of  the  absorbed  energy  is  radiated 
during  the  time  p/u.  The  results  can  be  plotted  in  a form  similar  to  that 
used  by  Nielson  and  Canavan  by  expressing  ( as  a function  of  1 and  p through 
Eq.  (D-7).  This  is  shown  in  Figure  112  where,  as  dashed  curves  corre- 
sponding to  q = 0.2  and  q = 1. 
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Figure  D-3  Radiative  Losses  in  Detonation  Wave 
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